Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Surface fluxes”

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 541 records · Page 30

Radiometric Assessment of the First Three Years of the NOAA-20 VIIRS Reflective Solar Bands Calibration

NASA’s Clouds and the Earth’s Radiant Energy System (CERES) SYN1deg Ed4.1 product utilizes geostationary (GEO) satellite measured radiances and retrieved cloud properties to account for the regional diurnal fluctuations in the Earth’s radiant broadband fluxes for times between the CERES measurements gathered from the Aqua (1:30 PM) and Terra (10:30 AM) sun-synchronous satellites. In CERES Edition 4 products, a global uniformity in the cloud properties and computed surface fluxes across the GEO satellite domains is maintained by scaling the radiance observations from more than twenty GEO visible imagers in the CERES record to a common radiometric reference scale, i.e., Aqua-MODIS. With the new-generation GEO (Himawari-8/9 and GOES-16/17) imagers having multiple reflective solar bands (RSB) that are spectrally similar to those of VIIRS, the CERES Imager and Geostationary Calibration Group (IGCG) is preparing to use NOAA-20 VIIRS as the reference imager for maintaining the radiometric uniformity across the GEO imager constellation. Given the recent Aqua satellite anomaly, this transition may happen sooner than the projected date of Aqua’s de-orbit. This paper presents an independent performance evaluation of the first three years of the NOAA-20 VIIRS RSB calibration in the NASA VIIRS Land Science Investigator-led Processing System (Land SIPS) L1b Collection 2 dataset. The temporal radiometric stability is assessed using multiple invariant Earth targets, including tropical deep convective clouds and the Saharan desert. The invariant target anisotropic reflectance at the top of the atmosphere was modeled using five years of stable satellite observations acquired from the previous VIIRS instrument onboard the Suomi National Polar-orbiting Partnership (S-NPP) satellite. The anisotropic corrections are essential for detecting temporal trends with high statistical confidence. In addition, the radiometric consistency between the RSB of the two VIIRS instruments will be evaluated.

Rajendra Bhatt↗

On the Potential Impact of Daytime Surface Sensible Heat Flux on the Dissipation of Martian Cold Air Outbreaks

The Martian daytime soil surface temperature is governed primarily by the net irradiance balance and surface soil heat flux. Thus the outbreak of a cold air mass generates increased sensible heat flux that is conducive to daytime dissipation of the cold air mass thermal characteristics. Conceptual and scaling evaluations of this dissipation are provided while comparison is made with similar situations on Earth. It is estimated that sensible heat flux contribution to the dissipation of the original thermal structure of the cold air could be three times larger than the corresponding situation on Earth. Illustrative numerical model simulations provide scaling of the potential impact on the dissipation of cold air masses for various combinations of background wind speed and latitudes.

Segal, M.↗

On the coronae of rapidly rotating stars. IV - Coronal activity in F dwarfs and implications for the onset of the dynamo

X-ray observations of 14 early F dwarfs are reported and these stars are used, together with a complete sample from the literature, to examine how the characteristics of coronal X-ray emission vary from dwarfs of spectral type A through G. Evidence for a rotation-activity relation in stars redder than B - V = 0.45 is found. Stellar duplicity and age, except insofar as they influence the rotation rate, do not appear to be important in determining the coronal X-ray flux level in F dwarfs. It is suggested that the appearance of a relation between rotation and activity at B - V = 0.45 indicates the turn-on of a solar-like dynamo. The high X-ray surface fluxes and small variance thereof for dwarfs with B - V = 0.3-0.45 are also discussed.

Walter, F. M.↗

High resolution absolute flux profiles of the MC 2 h and k lines in evolved F8 to M5 stars

The central results of a survey of the Mg II resonance line emission in a sample of over 50 evolved late type stars, including spectral-luminosity type F8 to M5 and La to IV are presented. Observed and surface fluxes are derived and correlations noted. The major findings include: (1) Mg II k emission core asymmetry transition near K1 III, analogous to that known for Ca II K; (2) a small gravity and temperature dependence of the Mg II chromospheric radiative loss rate.

Stencel, R. E.↗

Impact of Assimilating GOES-Derived Land Surface Variables into the PSU/NCAR MM5

Land surface processes are known to have a profound impact on the overlying atmosphere over a wide range of spatial and temporal scales. Many atmospheric numerical models include special parameterizations to improve the specification and partitioning of surface fluxes which are critical to the accurate prediction of warm season boundary layer behavior, organized mesoscale circulations, and convective precipitation. However, the added degrees of freedom resulting from the inclusion of vegetation and soil schemes require the specification of additional surface parameters such as vegetative resistances, green vegetation fraction, leaf area index, soil physical and hydraulic characteristics, and the vertical distribution of soil moisture. As satellite data have become more readily available in recent years, many investigations have attempted to use these new measurements to infer missing components of the surface energy budget. Sensitivity studies have shown land-skin temperature (LST) tendencies during the mid-morning hours are strongly sensitive to the surface moisture availability (a function of soil wetness and vegetation) and less sensitive to other parameters such as surface roughness. Based upon results from these and other studies, developed a simple technique that dynamically assimilates Geostationary Operational Environmental Satellite (GOES) derived land-surface products into the surface energy budget of a mesoscale model. The purpose of this paper is to demonstrate that assimilating the GOES satellite data has the potential to improve the representation of land surface characteristics within the model without prior knowledge of the land surface characteristics. The assimilation technique is presented in Section 2 and the numerical experiments are detailed in Section 3. Preliminary results and conclusions are presented in Sections 4 and 5, respectively.

Lapenta, William M.↗

Measurements of plasma parameters in the divertor island of Wendelstein 7-X through line-ratio spectroscopy on helium

We report electron temperature T e and density n e have been measured in the divertor island of Wendelstein 7-X by means of a new thermal He-beam system utilizing line-ratio spectroscopy. This system is equipped with a poloidal arrangement of five gas valves that allows to infer 2D plasma profiles T e ( R,Z ) and n e ( R,Z ) across the magnetic island forming the island divertor. A local maximum in T e ( R,Z ) was measured along the flux surfaces located between the island center and the outer separatrix. In the island center, characterized by closed field lines, a clear local minimum in T e was measured, indicating that no direct heat deposition occurs in this region and that the perpendicular heat flux from the core plasma into this domain is small. Increasing the island size leads to a reduction of the connection length in the scrape-off layer (SOL) and to a shift of the island center towards the divertor target. In this scenario, the T e ( R,Z ) are reduced in the entire island domain. Density profiles are found flat across the divertor island for both island size scenarios. These findings are supported by similar observations made with a reciprocating probe plunging through the same SOL island at the mid-plane of W7-X representing the upstream position in this 3D divertor geometry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Overview of the TJ-II stellarator research programme towards model validation in fusion plasmas

TJ-II stellarator results on modelling and validation of plasma flow asymmetries due to on-surface potential variations, plasma fuelling physics, Alfvén eigenmodes (AEs) control and stability, the interplay between turbulence and neoclassical (NC) mechanisms and liquid metals are reported. Regarding the validation of the neoclassically predicted potential asymmetries, its impact on the radial electric field along the flux surface has been successfully validated against Doppler reflectometry measurements. Research on the physics and modelling of plasma core fuelling with pellets and tracer encapsulated solid pellet injection has shown that, although post-injection particle radial redistributions can be understood qualitatively from NC mechanisms, turbulence and fluctuations are strongly affected during the ablation process. Advanced analysis tools based on transfer entropy have shown that radial electric fields do not only affect the radial turbulence correlation length but are also capable of reducing the propagation of turbulence from the edge into the scrape-off layer. Direct experimental observation of long range correlated structures show that zonal flow structures are ubiquitous in the whole plasma cross-section in the TJ-II stellarator. Alfvénic activity control strategies using ECRH and ECCD as well as the relation between zonal structures and AEs are reported. Finally, the behaviour of liquid metals exposed to hot and cold plasmas in a capillary porous system container was investigated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Annual Cycle of Planetary Boundary Layer Heights estimated from NOAA Profiler Network Data

The Planetary Boundary Layer (PBL) height over the Unites States Great and Central Plains during 1992-2012 is examined here using a combination of Wind Profiler-derived (WP) PBL height estimates and reanalysis fields from the Modern Era Retrospective Reanalysis Version 2 (MERRA-2). The combined analysis allows process study of the reasons behind the monthly mean behavior of the observed PBL heights under clear-sky conditions as well as the PBL height variability. WP PBL height monthly mean annual cycles were grouped into general categories of behavior, each analyzed using MERRA-2 fields of sensible and latent heat flux, surface temperature, net radiation and soil moisture. In the ’canonical’ category the latent heat plays little role in the determination of the monthly mean PBL height, and it follows the annual cycle of the net radiation. In the other categories, precipitation and latent heat flux had more influence in setting the annual cycle. An analysis of variance revealed that the role of latent heat in determining the PBL height variations is large (explaining up to 40% of PBL height variability) even in the ’canonical’ category for which latent heat played no role in setting the monthly mean. In other categories the latent heat explained up to 80% of the PBL height variations. The amount of that influence is shown to be related to the variability of column soil moisture.

MERRA-2↗

Influence of Idealized Heterogeneity on Wet and Dry Planetary Boundary Layers Coupled to the Land Surface

We examine the influence of surface heterogeneity on boundary layers using a large-eddy simulation coupled to a land-surface model. Heterogeneity, imposed in strips varying from 2-30 km (1 less than lambda/z(sub i) less than 18), is found to dramatically alter the structure of the free convective boundary layer by inducing significant organized circulations. A conditional sampling technique, based on the scale of the surface heterogeneity (phase averaging), is used to identify and quantify the organized surface fluxes and motions in the atmospheric boundary layer. The impact of the organized motions on turbulent transport depends critically on the scale of the heterogeneity lambda, the boundary layer height zi and the initial moisture state of the boundary layer. Dynamical and scalar fields respond differently as the scale of the heterogeneity varies. Surface heterogeneity of scale 4 less than lamba/z(sub i) less than 9 induces the strongest organized flow fields (up, wp) while heterogeneity with smaller or larger lambda/z(sub i) induces little organized motion. However, the organized components of the scalar fields (virtual potential temperature and mixing ratio) grow continuously in magnitude and horizontal scale, as lambda/z(sub i) increases. For some cases, the organized motions can contribute nearly 100% of the total vertical moisture flux. Patch-induced fluxes are shown to dramatically impact point measurements that assume the time-average vertical velocity to be zero. The magnitude and sign of this impact depends on the location of the measurement within the region of heterogeneity.

Houser, Paul↗

Suppressing trapped-electron-mode-driven turbulence via optimization of three-dimensional shaping

Abstract Turbulent transport driven by trapped electron modes (TEMs) is believed to drive significant heat and particle transport in quasihelically symmetric stellarators. Two three-dimensionally-shaped magnetic configurations with suppressed TEM-driven turbulence were generated through optimization that targeted quasihelical symmetry and the available energy of trapped electrons. Initial equilibria have flux surface shapes with a helically rotating negative triangularity (NT) and positive triangularity (PT). In gyrokinetic simulations, TEMs are suppressed in the reduced-TEM NT and PT configurations, showing that NT does not have the same beneficial turbulence properties over PT as seen in tokamaks. Heat fluxes from TEMs are also suppressed. Without temperature gradients and with a strong density gradient, the most unstable modes at low k y were consistent with toroidal universal instabilities (UIs) in the NT case and slab UIs in the PT case. Nonlinear simulations show that UIs drive substantial heat flux in both the NT and PT configurations. A moderate increase in β halves the heat flux in the NT configuration, while suppressing the heat flux in the PT geometry. Based on the present work, future optimizations aimed at reducing electrostatic drift wave-driven turbulent transport will need to consider UIs if β is sufficiently small.

Duff, J. M. (ORCID:0000000153086694)↗

GEWEX Water and Energy Budget Study

Closing the global water and energy budgets has been an elusive Global Energy and Water-cycle Experiment (GEWEX) goal. It has been difficult to gather many of the needed global water and energy variables and processes, although, because of GEWEX, we now have globally gridded observational estimates for precipitation and radiation and many other relevant variables such as clouds and aerosols. Still, constrained models are required to fill in many of the process and variable gaps. At least there are now several atmospheric reanalyses ranging from the early National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR) and NCEP/Department of Energy (DOE) reanalyses to the more recent ERA40 and JRA-25 reanalyses. Atmospheric constraints include requirements that the models state variables remain close to in situ observations or observed satellite radiances. This is usually done by making short-term forecasts from an analyzed initial state; these short-term forecasts provide the next guess, which is corrected by comparison to available observations. While this analysis procedure is likely to result in useful global descriptions of atmospheric temperature, wind and humidity, there is no guarantee that relevant hydroclimate processes like precipitation, which we can observe and evaluate, and evaporation over land, which we cannot, have similar verisimilitude. Alternatively, the Global Land Data Assimilation System (GLDAS), drives uncoupled land surface models with precipitation, surface solar radiation, and surface meteorology (from bias-corrected reanalyses during the study period) to simulate terrestrial states and surface fluxes. Further constraints are made when a tuned water balance model is used to characterize the global runoff observational estimates. We use this disparate mix of observational estimates, reanalyses, GLDAS and calibrated water balance simulations to try to characterize and close global and terrestrial atmospheric and surface water and energy budgets to within 10-20% for long term (1986-1995), large-scale global to regional annual means.

Roads, J.↗

Carbon Monitoring System Flux Estimation and Attribution: Impact of ACOS-GOSAT X(CO2) Sampling on the Inference of Terrestrial Biospheric Sources and Sinks

Using an Observing System Simulation Experiment (OSSE), we investigate the impact of JAXA Greenhouse gases Observing SATellite 'IBUKI' (GOSAT) sampling on the estimation of terrestrial biospheric flux with the NASA Carbon Monitoring System Flux (CMS-Flux) estimation and attribution strategy. The simulated observations in the OSSE use the actual column carbon dioxide (X(CO2)) b2.9 retrieval sensitivity and quality control for the year 2010 processed through the Atmospheric CO2 Observations from Space algorithm. CMS-Flux is a variational inversion system that uses the GEOS-Chem forward and adjoint model forced by a suite of observationally constrained fluxes from ocean, land and anthropogenic models. We investigate the impact of GOSAT sampling on flux estimation in two aspects: 1) random error uncertainty reduction and 2) the global and regional bias in posterior flux resulted from the spatiotemporally biased GOSAT sampling. Based on Monte Carlo calculations, we find that global average flux uncertainty reduction ranges from 25% in September to 60% in July. When aggregated to the 11 land regions designated by the phase 3 of the Atmospheric Tracer Transport Model Intercomparison Project, the annual mean uncertainty reduction ranges from 10% over North American boreal to 38% over South American temperate, which is driven by observational coverage and the magnitude of prior flux uncertainty. The uncertainty reduction over the South American tropical region is 30%, even with sparse observation coverage. We show that this reduction results from the large prior flux uncertainty and the impact of non-local observations. Given the assumed prior error statistics, the degree of freedom for signal is approx.1132 for 1-yr of the 74 055 GOSAT X(CO2) observations, which indicates that GOSAT provides approx.1132 independent pieces of information about surface fluxes. We quantify the impact of GOSAT's spatiotemporally sampling on the posterior flux, and find that a 0.7 gigatons of carbon bias in the global annual posterior flux resulted from the seasonally and diurnally biased sampling when using a diagonal prior flux error covariance.

biased sampling↗

Experimental demonstration of transient CHI start-up using a floating biased electrode configuration

Transient coaxial helicity injection (CHI) current start-up in a spherical tokamak (ST) has been demonstrated for the first time using a reactor-relevant floating single biased (FSB) electrode configuration on QUEST. Unlike previous applications of CHI on HIT-II and on NSTX in which the entire ST vessel was split in two using two toroidal ceramic rings as vacuum breaks, in the FSB electrode configuration on QUEST, only the lower outer divertor plate is electrically insulated from the rest of the vessel. This is the CHI injector region. A CHI discharge is generated by driving current, using an external power supply, along magnetic field lines that connect the inner and outer divertor plates in the CHI injector region. The CHI plasma evolution and closed flux formation on QUEST are very similar to those on NSTX and have generated toroidal currents of 150 kA with 100 kA of current flowing on closed flux surfaces. An alternate, force-free model approach has been used to derive the CHI scaling relations. The resulting CHI scaling relations from the model are consistent with the conventional CHI scaling relations, which state that the CHI generated toroidal current is proportional to the CHI injector flux and the CHI injector current is inversely proportional to the toroidal field. These results from QUEST bode well for the application of CHI in a new generation of compact high field STs and low aspect ratio tokamaks in which the space for the central solenoid is very restricted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Remote sensing of atmosphere and oceans; Proceedings of Symposium 1 and of the Topical Meeting of the 27th COSPAR Plenary Meeting, Espoo, Finland, July 18-29, 1988

Papers are presented on the contribution of space remote sensing observations to the World Climate Research Program and the Global Change Program, covering topics such as space observations for global environmental monitoring, experiments related to land surface fluxes, studies of atmospheric composition, structure, motions, and precipitation, and remote sensing for oceanography, observational studies of the atmosphere, clouds, and the earth radiation budget. Also, papers are given on results from space observations for meteorology, oceanography, and mesoscale atmospheric and ocean processes. The topics include vertical atmospheric soundings, surface water temperature determination, sea level variability, data on the prehurricane atmosphere, linear and circular mesoscale convective systems, Karman vortex clouds, and temporal patterns of phytoplankton abundance.

Raschke, E.↗

Possible Influences of Air Pollution, Dust and Sandstorms on the Indian Monsoon

In Asian monsoon countries, such as China and India, human health and safety problems caused by air pollution are becoming increasingly serious, due to the increased loading of atmospheric pollutants from waste gas emissions and from rising energy demand associated with the rapid pace of industrialization and modernization. Meanwhile, uneven distribution of monsoon rain associated with flash floods or prolonged drought, has caused major loss of human life and damage to crops and.property with devastating societal impacts. Historically, air-pollution and monsoons research are treated as separate problems. However recent studies have suggested that the two problems may be intrinsically linked and need to be studied jointly. Fundamentally, aerosols can affect precipitation through radiative effects cif suspended particles in the atmosphere (direct effect) and/or by interfering and changing: the cloud and precipitation formation processes (indirect effect). Based on their optical properties, aerosols can be classified into two types.: those that absorb solar radiation, and those that do not. Both types of aerosols scatter sunlight and reduce the amount of solar radiation from reaching the Earth's surface, causing it to cool. The surface cooling increases atmospheric stability and reduces convection potential, Absorbing aerosols, however, in addition to cooling the surface, can heat the atmosphere. The heating of the atmosphere may reduce the amount of low clouds by increased evaporation in cloud drops. The heating, however, may induce rising motion, enhance low-level moisture, convergence and, hence, increases rainfall, The latent heating from enhanced rainfall may excite feedback processes in the large-scale circulation, further amplify.the initial response to aerosol heating and producing more rain. Additionally, aerosols can increase the concentration of cloud condensation nuclei (CCN), increase cloud amount and decrease coalescence and collision rates, leading to reduced precipitation. However, in the presence of increasing moist and warm air, the reduced coalescence/collision may lead to supercooled drops at higher altitudes where ice precipitation falls and melts. The latent heat release from freezing aloft and melting below implies greater upward heat transport in polluted clouds and invigorate deep convection. In this way, aerosols may lead to increased local convection. Hence, depending on the ambient large-scale conditions and dynamical feedback processes, aerosols' effect on precipitation can be positive, negative or mixed. In the Asian monsoon and adjacent regions, the aerosol forcing and responses of the water cycle are even more complex, Both direct and indirect effects may take place locally and simultaneously, interacting with each other. in addition to local effects, monsoon rainfall may be affected by aerosols transported from other regions and intensified through large-scale circulation and moisture feedback. Thus, dust transported by the large-scale circulation from the adjacent deserts to northern India may affect rainfall over the Bay of Bengal; sulphate and black carbon front industrial pollution in central, southern China and northern India may affect the rainfall regime over the Korean peninsula and Japan; organic and black carbon front biomass burning from Indo-China may modulate the pre-monsoon rainfall regime over southern China and coastal regions, contributing to variability in differential heating and cooling of the atmosphere and to the land-sea thermal contrast. During the pre-monsoon season and monsoon breaks, it has been suggested that radiative forcing by absorbing aerosols have nearly the same order of magnitude as the forcing due to latent heating from convection and surface fluxes. The magnitude of the total aerosol radiative cooling due to sulphates and soot is of the order of 20-40 W/m2 over the Asian monsoon land region in the pre-monsoon season, compared to about 1-2 W/m2 for global warng. However, the combined forcing at the surface and in the atmosphere, including all species. if aerosols, and details of aerosol mixing, and impacts on the energy and water cycles in the monsoon land regions, are not well known.

Lau, William K. M.↗

Resilient stellarator divertor characteristics in the Helically Symmetric eXperiment

Resilient divertor features connected to open chaotic edge structures in the Helically Symmetric eXperiment are investigated. For the first time, an expanded vessel wall was considered that would give space for implementation of a physical divertor target structure. The analysis was done for four different magnetic configurations with very different chaotic plasma edges. A resilient plasma wall interaction pattern was identified across all configurations. This manifests as qualitatively very similar footprint behavior across the different plasma equilibria. Overall, the resilient field lines of interest with high connection length L C lie within a helical band along the wall for all configurations. This resiliency can be used to identify the best location of a divertor. The details of the magnetic footprint’s resilient helical band is subject to specific field line structures which are linked to the penetration depth of field lines into the plasma and directly influence the heat and particle flux patterns. The differences arising from these details are characterized by introducing a new metric, the minimum radial connection min $(\delta_N)$ of a field line from the last closed flux surface. The relationship, namely the deviation from a scaling law, between min $(\delta_N)$ and L C of the field lines in the plasma edge field line behavior suggests that the field lines are associated with structures such as resonant islands, cantori, and turnstiles. This helps determine the relevant magnetic flux channels based on the radial location of these chaotic edge structures and the divertor target footprint. These details will need to be taken into account for resilient divertor design.

divertor↗

Impurity temperature screening in stellarators close to quasisymmetry

Impurity temperature screening is a favourable neoclassical phenomenon involving an outward radial flux of impurity ions from the core of fusion devices. Quasisymmetric magnetic fields lead to intrinsically ambipolar neoclassical fluxes that give rise to temperature screening for low enough $\unicode[STIX]{x1D702}^{-1}\equiv d\ln n/d\ln T$ . In contrast, neoclassical fluxes in generic stellarators will depend on the radial electric field, which is predicted to be inward for ion-root plasmas, potentially leading to impurity accumulation. Here, we examine the impurity particle flux in a number of approximately quasisymmetric stellarator configurations and parameter regimes while varying the amount of symmetry breaking in the magnetic field. For the majority of this work, neoclassical fluxes have been obtained using the SFINCS drift-kinetic equation solver with electrostatic potential $\unicode[STIX]{x1D6F7}=\unicode[STIX]{x1D6F7}(r)$ , where $r$ is a flux-surface label. Results indicate that achieving temperature screening is possible, but unlikely, at low-collisionality reactor-relevant conditions in the core. Thus, the small departures from symmetry in nominally quasisymmetric stellarators are large enough to significantly alter the neoclassical impurity transport. A further look at the magnitude of these fluxes when compared to a gyro-Bohm turbulence estimate suggests that neoclassical fluxes are small in configurations optimized for quasisymmetry when compared to turbulent fluxes.

Physics↗

Cloud‐Surface Coupling Alters the Morning Transition From Stable to Unstable Boundary Layer

Abstract Due to surface heating, the morning boundary layer transits from stable to neutral or convective conditions, exerting critical influences on low tropospheric thermodynamics. Low clouds closely interact with the boundary layer development, yet their interactions bear considerable uncertainties. Our study reveals that cloud‐surface coupling alters the morning transition from stable to unstable boundary layer and thus notably affects the diurnal variation of the boundary layer. Specifically, due to the reduction in surface fluxes, decoupled clouds can delay the process of eroding nocturnal inversion by 0.8‐hr and even prevent the transition of the boundary layer from happening for 12% of decoupled cases, keeping the boundary layer in a stable state during the noontime. On the other hand, when clouds are coupled with the surface, cloud‐top radiative cooling can directly cool the upper boundary layer to facilitate sub‐cloud convection, leading to an unstable boundary layer in the earlier morning.

54 ENVIRONMENTAL SCIENCES↗