Search NASA⌕ Search

SEARCH · Search NASA

Results for “entrainment processes”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Aerosol-boundary layer interaction modulated entrainment process

Abstract Aerosol-boundary layer interactions play an important role in affecting atmospheric thermodynamics and air pollution. As a key factor in dictating the development of the boundary layer, the entrainment process in the context of aerosol-boundary layer interactions is still poorly understood. Using comprehensive field observations made at a superstation in Beijing, we gain insight into the response of the entrainment process to aerosols. We found that high aerosol loading can significantly suppress the entrainment rate, breaking the conventional linear relationship between sensible heat fluxes and entrainment fluxes. Related to aerosol vertical distributions, aerosol heating effects can alter vertical heat fluxes, leading to a strong interaction between aerosols and the entrainment process in the upper boundary layer. Such aerosol-entrainment coupling can inhibit boundary layer development and explains the great sensitivity of observed entrainment rates to aerosols than can traditional calculations. The notable impact of aerosols on the entrainment process raises holistic thinking about the dynamic framework of the boundary layer in a polluted atmosphere, which may have a significant bearing on the dispersion of air pollutants and the land-atmosphere coupling.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1982

Clouds and entrainment

The concept of entrainment is discussed, and the knowledge that has been acquired about entraining turbulence in clear air is surveyed. Attention is then given to observations and models of stratocumulus entrainment, which show some of the ways in which cloudiness modifies the entrainment process. Entraining cumuli are also considered, with some of the differences between cumulus entrainment and stratocumulus entrainment noted. Emphasis is placed on the dynamical aspects of entrainment; its effects on droplet spectra and other microphysical aspects are mentioned only briefly.

Randall, D. A.↗
DOE OSTI2021

Vertical variation of turbulent entrainment mixing processes in marine stratocumulus clouds using high-resolution digital holography

Marine stratocumulus clouds contribute signi?cantly to the Earth’s radiation budget due to their extensive coverage and high albedo. Yet, subgrid variability in cloud properties such as aerosol concentration, droplet number and precipitation rates lead to considerable errors in global climate models. While these clouds usually have small vertical ex-tent, turbulent entrainment-mixing and precipitation can generate signi?cant variations in droplet number, size and relative dispersion with altitude. In this paper, we analyze turbulent entrainment-mixing processes and the variability in cloud microphysical properties as a function of height within a warm marine stratocumulus cloud layer over the Eastern North Atlantic. We use high resolution airborne holographic measurements and compare them with local turbulence measurements. We ?nd that entrainment-mixing is primarily inhomogeneous near cloud top leading to larger droplet sizes and homogeneous near cloud base leading to smaller droplet sizes. Further analysis of Damk¨ohler number measurements are able to explain the mixing mechanisms at di?erent cloud heights, reinforcing the importance of turbulent mixing and microphysical time scales in deter-mining cloud microphysics.

Desai, Neel↗
DOE OSTI2024

Using Machine Learning to Predict Cloud Turbulent Entrainment–Mixing Processes

Different turbulent entrainment–mixing mechanisms between clouds and environment are essential to cloud–related processes; however, accurate representation of entrainment–mixing in weather/climate models still poses a challenge. This study exploits the use of machine learning (ML) to address this challenge. Four ML (Light Gradient Boosting Machine [LGB], eXtreme Gradient Boosting, Random Forest, and Support Vector Regression) are examined and compared. It is found that LGB performs best, and thus is selected to understand the impact of entrainment–mixing on microphysics using simulation data from Explicit Mixing Parcel Model. Compared with traditional parameterizations, the trained LGB provides more accurate microphysical properties (number concentration and cloud droplet spectral dispersion). The partial dependences of predicted microphysics on features exhibit a strong alignment with physical mechanisms and expectations, as determined by the interpreting method, thus overcoming the limitations of the “black box” scheme. The underlying mechanisms are that the smaller number concentration and larger spectral dispersion correspond to more inhomogeneous entrainment–mixing. Specifically, number concentration after entrainment–mixing is positively correlated with adiabatic number concentration and liquid water content affected by entrainment–mixing, and inversely correlated with adiabatic volume mean radius. Spectral dispersion after entrainment–mixing is negatively correlated with liquid water content affected by entrainment–mixing, turbulent dissipation rate and relative humidity of entrained air. Sensitivity analysis further suggests that number concentration is mainly determined by cloud microphysical properties whereas spectral dispersion is influenced by both cloud microphysical properties and environmental variables. The results indicate that the LGB scheme has the potential to enhance the representation of entrainment–mixing in weather/climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Quantifying the Turbulent Entrainment‐Mixing Processes Based on Z ‐LWC Relationships of Cloud Droplets

Abstract Turbulent entrainment‐mixing processes profoundly influence the relationship between radar reflectivity factor and liquid water content ( Z ‐LWC) of cloud droplets. However, quantification of the entrainment‐mixing mechanisms based on the Z ‐LWC relationship is still lacking. To address this gap, 12,218 entrainment‐mixing cases are simulated using the Explicit Mixing Parcel Model. We examine the variations of the parameters in the power‐law relationship Z = a LWC b , and the relationship between parameter b and homogeneous mixing degree ( ψ ), a measure quantifying entrainment‐mixing processes. The results indicate that parameter b distributes within the range of 1–2, with a positive correlation between parameter b and ψ . The b ‐ ψ relationship is fitted, which connects the Z ‐LWC relationship for various entrainment‐mixing types. The results suggest the potential for employing a remote sensing approach to investigate the entrainment‐mixing mechanisms of non‐precipitating small cumulus/stratocumulus clouds, thereby overcoming the limitations of traditional observational studies that rely solely on aircraft observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Dust-planetary boundary layer interactions amplified by entrainment and advections

Mineral dust contributes to more than half of the global aerosol loading. However, the radiative impacts of dust aerosols on planetary boundary layer (PBL) structure have not been explored sufficiently. During a typical dust storm event over Tarim Basin, dust aerosols exhibit a well-mixed distribution during the daytime in spite of a shallow layer of dust particles accumulated at higher altitudes. By contrast, nocturnal dust plumes are located near the surface due to stable stratification. We demonstrate that these differentiated vertical distributions determine the spatial heterogeneity of dust loading, radiative fluxes and PBL height variations. Dust aerosols cause daytime PBL suppression and nighttime PBL promotion through modulating surface and atmospheric radiative budgets. Specifically, dust-induced cooling effect within PBL directly inhibits the daytime PBL development. PBL suppression effect is then amplified by entrainment processes resulting in excessively low PBL height, especially for dust particles below but near the PBL top. Dust plumes weaken updrafts from PBL and downdrafts of the free atmosphere, which further reduce the entrainment mixing through attenuating horizontal and vertical advection, and eventually amplify PBL suppression. At night, near-surface dust aerosols stimulate a warm and unstable lower atmosphere, generate warm advection heating and promote the PBL development. Our study highlights the importance of specifying entrainment parameters and nighttime advection activities in quantifying the dust-PBL interactions.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1990

Laboratory experiments on stability and entrainment of oceanic stratocumulus. Part 1: Instability experiment

The existence and persistence of marine stratocumulus play a significant role in the overall energy budget of the earth. Their stability and entrainment process are important in global climate studies, as well as for local weather forecasting. The purposes of the experimental simulations are to study this process and to address this paradox. The effects of buoyancy reversal is investigated, followed by two types of experiments. An instability experiment involves the behavior of a fully turbulent wake near the inversion generated by a sliding plate. Due to buoyancy reversal, the heavy, mixed fluid starts to sink, turning the potential energy created by the mixing process into kinetic energy, thereby increasing the entrainment rate. An entrainment experiment, using a vertically oscillating grid driven by a controllable speed motor, produces many eddy-induced entrainments at a surface region on scales much less than the depth of the layer.

Shy, Shenqyang S.↗
NASA NTRS2000

The Turbulent/Non-Turbulent Interface Bounding a Far-Wake

The velocity fields of a turbulent wake behind a flat plate obtained from the direct numerical simulations of Moser et al. are used to study the structure of the flow in the intermittent zone where there are, alternately, regions of fully turbulent flow and non-turbulent velocity fluctuations either side of a thin randomly moving interface. Comparisons are made with a wake that is 'forced' by amplifying initial velocity fluctuations. There is also a random temperature field T in the flow; T varies between constant values of 0.0 and 1.0 on the sides of the wake. The value of the Reynolds number based on the centreplane mean velocity defect and halfwidth b of the wake is Re approx. = 2000. It is found that the thickness of the continuous interface is about equal to 0.07b, whereas the amplitude of fluctuations of the instantaneous interface displacement y(sub I)(t) is an order of magnitude larger, being about 0.5b. This explains why the mean statistics of vorticity in the intermittent zone can be calculated in terms of the probability distribution of y(sub I) and the instantaneous discontinuity in vorticity across the interface. When plotted as functions of y - y(sub I), the conditional mean velocity (U) and temperature (T) profiles show sharp jumps Delta(U) and Delta(T) at the interface adjacent to a thick zone where (U) and (T) vary much more slowly. Statistics for the vorticity and velocity variances, available in such detail only from DNS data, show how streamwise and spanwise components of vorticity are generated by vortex stretching in the bulges of the interface. Flow fields around the interface, analyzed in terms of the local streamline pattern, confirm previous results that the advancement of the vortical interface into the irrotational flow is driven by large-scale eddy motion. It is argued that because this is an inviscid mechanism the entrainment process is not sensitive to the value of Re, and that small-scale nibbling only plays a subsidiary role. While mean Reynolds stresses decrease gradually in the intermittent zone, conditional stresses are found to decrease sharply towards zero at the interface. Using one-point turbulence models applied to either unconditional or conditional statistics for the turbulent region and then averaged, the entrainment rate E(sub b) would, if calculated exactly, be zero. But if computed with standard computational methods, E(sub b) would be non-zero because of numerical diffusion. It is concluded that the current practice in statistical models of approximating entrainment by a diffusion process is computationally arbitrary and physically incorrect. An analysis shows how E(sub b) is related to Delta(U) and the jump in shear stress at the interface, and correspondingly to Delta(T) and the heat flux.

Bisset, David K.↗
NASA NTRS2005

Numerical Analysis of Flow Evolution in a Helium Jet Injected into Ambient Air

A computational model to study the stability characteristics of an evolving buoyant helium gas jet in ambient air environment is presented. Numerical formulation incorporates a segregated approach to solve for the transport equations of helium mass fraction coupled with the conservation equations of mixture mass and momentum using a staggered grid method. The operating parameters correspond to the Reynolds number varying from 30 to 300 to demarcate the flow dynamics in oscillating and non-oscillating regimes. Computed velocity and concentration fields were used to analyze the flow structure in the evolving jet. For Re=300 case, results showed that an instability mode that sets in during the evolution process in Earth gravity is absent in zero gravity, signifying the importance of buoyancy. Though buoyancy initiates the instability, below a certain jet exit velocity, diffusion dominates the entrainment process to make the jet non-oscillatory as observed for the Re=30 case. Initiation of the instability was found to be dependent on the interaction of buoyancy and momentum forces along the jet shear layer.

Satti, Rajani P.↗
DOE OSTI2026

A Model Intercomparison Study to Investigate Mixing Characteristics in Non‐Precipitating Stratocumulus Clouds

Recent aircraft observations of marine stratocumulus clouds consistently showed that cloud microphysical relationships vary with altitude, indicating inhomogeneous mixing characteristics near cloud top and homogeneous mixing characteristics in mid-levels of clouds. Here, we conduct model intercomparison of an idealized, non-precipitating stratocumulus cloud to evaluate model consistency and examine whether simulations can reproduce the observed mixing characteristics. The results show that eleven large-eddy simulations with various dynamics and microphysics schemes show good agreement on the thermodynamical, microphysical, and dynamical properties of the stratocumulus-topped boundary layer in a steady state. The inter-model spread in steady-state liquid water path is significantly reduced compared to previous model intercomparison studies. This improvement might be due to better models and more consistent initial conditions than those used decades ago. In addition, most simulations, including a low-dimensional simulation, capture inhomogeneous mixing characteristics near the cloud top and homogeneous mixing characteristics inside the cloud. Moreover, simulations using Lagrangian microphysics schemes agree better with the observed mixing characteristics compared with those using the bin microphysics schemes. Since most simulations do not fully resolve the entrainment process, the apparent mixing characteristics arise from the variations in the resolved cloud properties. Our results support the vertical circulation mixing hypothesis, which suggests that homogeneous mixing characteristics in mid-levels of clouds are due to the vertical circulation of entrainment-affected and diluted parcels from the cloud top moved to lower levels.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A deep insight on the coal ash-to-slag transformation behavior during the entrained flow gasification process

Recent research provided deep insight on the coal ash-to-slag transformation characterization during the entrained flow gasification process, with experimentation on a 40 kg/hr (dry basis) coal-fed opposed multi-burner (OMB) entrained flow gasifier and simulation via FactSage™ software. A general mechanism is presented to relate the gasifier design temperature, ash fluid temperature, and operating temperature with the degree of the slag polymerization. The change of the high temperature zone, the corresponding particle residence time in the high temperature zone, and syngas composition have obvious effects on the slag mineral transformation behavior. Mineral types formed on the wall of the gasifier chamber were mainly anorthite (CaAl 2 Si 2 O 8 ), aluminum oxide (Al 2 O 3 ), and calcium sulfide (CaS). These minerals transformed to anorthite and diopside (CaMgSi 2 O 6 ) at the slag hole zone, while the minerals at the lock hopper were anorthite, orthoclase (KAlSi 3 O 8 ), quartz (SiO 2 ), gypsum (CaSO 4 ), calcite (CaCO 3 ), and halite (NaCl). FactSage™ predicted minerals as anorthite, diopside, orthoclase, and albite (NaAlSi 3 O 8 ), etc., where the slag temperature was below the ash fluid temperature and when the ratios of CO/CO 2 and (CO + H 2 )/CO 2 were lower than 1.0 and 2.0, respectively. By simulation, residual carbon was found to be the dominant factor over syngas composition to cause mineral transformation, and this was verified experimentally. The Ca-based crystals, typically anorthite, was shifted to diopside, near the slag hole zone, and a linear relationship was found between the content ratios of diopside/(anorthite + diopside), CaO/SiO 2 , and (CaO + MgO)/SiO 2 . A dimensionless number, θ, was defined to characterize the changing chemical composition and the degree of slag polymerization, with temperature deviation from the design condition. Three zones of θ were identified and related to the deviation between the actual gasification condition from the design condition. Finally, a low slag polymerization degree corresponded with a higher temperature deviation between the actual condition and design condition, and this proved that increased residual carbon content and changing iron valence state increased the mineral types when the slag temperature was below the ash fluid temperature.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2025

Regime-based aerosol–cloud interactions from CALIPSO-MODIS and the Energy Exascale Earth System Model version 2 (E3SMv2) over the Eastern North Atlantic

This study investigates aerosol-cloud interactions in marine boundary layer (MBL) clouds using an advanced deep-learning-driven synoptic-regime-based framework, combining satellite data (CALIPSO vertically resolved aerosol extinction and MODIS cloud properties) with 1° nudged Energy Exascale Earth System Model version 2 (E3SMv2) simulation over the Eastern North Atlantic (ENA; ∼10°×10°, 2006–2014). The E3SMv2 captures observed seasonal variations in cloud droplet number concentrations (N d ) and liquid water path (LWP), though it systematically underestimates N d . We then partition ENA meteorology into four synoptic regimes (Pre-Trough, Post-Trough, Ridge, Trough) via a deep-learning clustering of ERA5 reanalysis fields, enabling regime-dependent aerosol-cloud interactions analyses. Both satellite and E3SMv2 exhibit an inverted-V LWP-N d relationship. In Post-Trough and Ridge regimes, the satellite shows stronger negative LWP-N d sensitivities than in Pre-Trough regime. The Trough regime displays a muted satellite LWP response. In comparison, the model predicts more exaggerated LWP responses across regimes, with LWP increasing too quickly at low N d and decreasing more sharply at high N d , especially in Pre-Trough and Trough regimes. These exaggerated model LWP sensitivities may stem from uncertainties in representing drizzle processes, entrainment, and turbulent mixing. As for N d susceptibility to aerosols, N d increases with MBL aerosol extinction in both datasets, but the simulated aerosol-cloud interactions appear oversensitive to meteorological conditions. Overall, E3SMv2 better captures aerosol effects under regimes that favor stratiform clouds (Post-Trough, Ridge), but performance deteriorates for regimes with deeper, dynamically complex clouds (Trough), highlighting the need for improved representations of those cloud processes in climate models.

Environmental sciences↗
NASA NTRS1982

Entrainment and detrainment in a simple cumulus cloud model

A cumulus cloud model, analogous to the mixed-layer models of the planetary boundary layer and the upper ocean, is developed using a single, unitary entrainment process in which the motion of the cloud boundary relative to the mean flow is permitted, produced, and controlled by turbulent processes. An alternate theory to the mixing-length theory of Asai and Kasahara (1967) is proposed which completely removes the strong scale-dependence of the Asai-Kasahara model. The model reintroduces scale-dependence by introducing including the pe5turbation pressure term of the equation of vertical motion. It is shown that the model predicts deeper clouds than the Asai-Kasahara model for a given sounding, due to the entrainment assumption and the effects of the perturbation pressure. Lateral entrainment dominates cloud-top entrainment, although finite-difference errors increase the cloud-top entrainment rate from zero to a positive value in actual situations. The fractional entrainment rate for updrafts is determined to vary only slightly with height and to decrease only slowly as the cloud radius increases, while the fractional detrainment rate for updrafts increases with height.

Randall, D. A.↗
NASA NTRS1990

Winds and turbulence above San Nicolas Island during FIRE

Winds and turbulence above San Nicolas Island (SNI) were continuously monitored during the First ISCCP Regional Experiment (FIRE) with a Doppler sodar (1600 Hz) and an ultrahigh frequencies (UHF) (400 MHz) wind profiler. Individual profiles were logged from 1 through 19 July 1987 at 30 sec to 2 min intervals; preliminary processing of hourly averaged winds from 100 to 2300 m is complete. These wind profiles have been combined with intermittent radiosonde profiles to construct time-height sections of parameters such as the Richardson Number. Analysis now in progress includes interpretation of the sodar (C 2/T) and profiler (C 2/N, C 2/q) reflectivity profiles. These profiles are being combined with radiometric and ceilometer data for investigations of vertical transport and cloud base/top entrainment processes. Finally, by using the SNI data in conjunction with soundings made from the research vessel Pt. Sur, the relationship between changes in the SNI winds and mesoscale and synoptic-scale weather features is being investigated.

White, Allen B.↗
NASA NTRS1986

Biological consequences of a recurrent eddy off Point Conception, California

The biological effects on three different time scales (100-day mesoscale, annual, and several-year) of a mesoscale anticyclonic eddy consistently found in shipboard surveys and satellite-sensed data several hundred kilometers southwest of Point Conception, CA, are described. A detailed shipboard study of the eddy in January 1981 found a complex system of fronts in surface chlorophyll at the northern edge of the eddy; microplankton and zooplankton distributions were strongly affected by entrainment processes at the surface and, apparently, at depth. Concurrent satellite coastal zone color scanner ocean color images show agreement with the general surface characteristics of the eddy chlorophyll field but do not reflect features deeper than about 25 m, including the contribution of the deep chlorophyll maximum to the integrated chlorophyll values. Satellite data for the period October 1980 through October 1981 and shipboard data from California Cooperative Oceanic Fisheries Investigations (CalCOFI) for December 1980 to July 1981 show the continued presence of the eddy in the sea surface temperature and color field and in the distributions of surface chlorophyll and zooplankton displacement volume. A review of the CalCOFI survey results from 1949 to the present time demonstrates the recurrent nature of the eddy system on a year-to-year basis. The eddy system appears to have a significant effect on the distribution of both oceanic and nearshore organisms. Offshore transport of coastal species occurs in the form of large entrained plumes or filaments.

Haury, Loren R.↗
DOE OSTI2023

Evaluation of liquid cloud albedo susceptibility in E3SM using coupled eastern North Atlantic surface and satellite retrievals

The impact of aerosol number concentration on cloud albedo is a persistent source of spread in global climate predictions due to multi-scale, interactive atmospheric processes that remain difficult to quantify. We use 5 years of geostationary satellite and surface retrievals at the US Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) eastern North Atlantic (ENA) site in the Azores to evaluate the representation of liquid cloud albedo susceptibility for overcast cloud scenes in the DOE Energy Exascale Earth System Model version 1 (E3SMv1) and provide possible reasons for model–observation discrepancies. The overall distribution of surface 0.2 % CCN concentration values is reasonably simulated, but simulated liquid water path (LWP) is lower than observed and layer mean droplet concentration (N d ) comparisons are highly variable depending on the N d retrieval technique. E3SMv1's cloud albedo is greater than observed for given LWP and N d values due to a lower cloud effective radius than observed. However, the simulated albedo response to N d is suppressed due to a correlation between the solar zenith angle (SZA) and N d created by the seasonal cycle that is not observed. Controlling for this effect by examining the cloud optical depth (COD) shows that E3SMv1's COD response to CCN concentration is greater than observed. For surface-based retrievals, this is only true after controlling for cloud adiabaticity because E3SMv1's adiabaticities are much lower than observed. Assuming a constant adiabaticity in surface retrievals as done in top-of-atmosphere (TOA) retrievals narrows the retrieved ln N d distribution, which increases the cloud albedo sensitivity to ln N d to match the TOA sensitivity. The greater sensitivity of COD to CCN is caused by a greater Twomey effect in which the sensitivity of N d to CCN is greater than observed for TOA-retrieved N d , and once model–observation cloud adiabaticity differences are removed, this is also true for surface-retrieved N d . The LWP response to N d in E3SMv1 is overall negative as observed. Despite reproducing the observed LWP–N d relationship, observed clouds become much more adiabatic as N d increases, while E3SMv1 clouds do not, associated with more heavily precipitating clouds that are partially but not completely caused by deeper clouds and weaker inversions in E3SMv1. These cloud property differences indicate that the negative LWP–N d relationship is likely not caused by the same mechanisms in E3SMv1 and observations. The negative simulated LWP response also fails to mute the excessively strong Twomey effect, highlighting potentially important confounding factor effects that likely render the LWP–N d relationship non-causal. N d retrieval scales and assumptions, particularly related to cloud adiabaticity, contribute to substantial spreads in the model–observation comparisons, though enough consistency exists to suggest that aerosol activation, drizzle, and entrainment processes are critical areas to focus E3SMv1 development for improving the fidelity of aerosol–cloud interactions in E3SM.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Evaluating Sea Breezes and Associated Convective Cloud Evolution in the Model Gray Zone

We characterize convective clouds associated with sea‐breeze circulations (SBC) using multi‐agency observations and multi‐case ensemble model simulations. The focus is on assessing convective cloud lifecycle properties and their merging behavior, as well as the environmental conditions they are embedded in, particularly SBC features. In total, 46 SBC days over the Houston‐Galveston region are selected and simulated using the Weather Research and Forecasting (WRF) model at a gray zone scale with a forecast‐like parameterization setup. Advanced techniques, including change‐point detection, a Lagrangian cloud tracking method, and a newly developed cell merging and splitting detection algorithm, are applied and/or developed for this study. Our findings indicate that the WRF model at 1 km grid spacing well represents the thermodynamic conditions over the region, as well as SBC timing and intensity. However, for the associated convective cells, WRF overestimates the 30‐dBZ echo top height, cell area, and maximum radar reflectivity compared to radar observations. This overestimation is potentially due to under‐resolved entrainment processes, an overestimated merging frequency, and the overestimation of updraft intensity. Furthermore, the model exhibits a deficiency in simulating congestus clouds, showing a more rapid transition from shallow to deep convection compared to observed behavior. Moreover, observations indicate stronger, deeper, and wider clouds when merging happens. Conversely, in simulations, the merging process does not necessarily lead to higher or longer‐lived cells, as many cases experience rapid and frequent merging and splitting which may result in more variance in convective updraft velocity during the convection lifetime.

54 ENVIRONMENTAL SCIENCES↗