Search NASA⌕ Search

SEARCH · Search NASA

Results for “troposphere”

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

Effects of Wind Shear and Aerosol Conditions on the Organization of Precipitating Marine Stratocumulus Clouds

Here this study examines how wind shear affects precipitating marine stratocumulus clouds under different cloud droplet number concentrations (N d ). We performed a series of large eddy simulations (LES) of nocturnal marine stratocumulus clouds using Cloud Model 1 (CM1). The simulations show that N d is the dominant factor for cloud cellular organization in this cloud regime rather than wind shear. Low N d characterizes the open cellular structure with a high in-cloud liquid water path (LWP). When wind shear is increased, the cloud fraction tends to decrease along with LWP, suggesting the cloud top region is significantly influenced by the entrainment and mixing of dry air from the free troposphere. We also examine cold pools in open and closed cellular clouds. Open-cell clouds produce larger and deeper cold pools compared to closed-cell clouds. Interestingly, cold pools can exist without surface precipitation and are produced by evaporation of light precipitation (drizzle) below the cloud base with weak downdrafts. The evaporation of raindrops and drizzle play an important role in initiating new convection, particularly where colliding outflows occur downstream of the cloud. This secondary convection contributes to the development and maintenance of the cloud cellular organization and formation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Variability of Eastern North Atlantic Summertime Marine Boundary Layer Clouds and Aerosols Across Different Synoptic Regimes Identified With Multiple Conditions

Abstract This study estimates the meteorological covariations of aerosol and marine boundary layer (MBL) cloud properties in the eastern North Atlantic (ENA) region, characterized by diverse synoptic conditions. Using a deep‐learning‐based clustering model with mid‐level and surface daily meteorological data, we identify seven distinct synoptic regimes during the summer from 2016 to 2021. Our analysis, incorporating reanalysis data and satellite retrievals, shows that surface aerosols and MBL clouds exhibit clear regime‐dependent characteristics, whereas lower tropospheric aerosols do not. This discrepancy likely arises from synoptic regimes determined by daily large‐scale conditions, which may overlook air mass histories that predominantly dictate lower tropospheric aerosol conditions. Focusing on three regimes dominated by northerly winds, we analyze the Atmospheric Radiation Measurement Program (ARM) ENA observations on Graciosa Island in the Azores. In the subtropical anticyclone regime, fewer cumulus clouds and more single‐layer stratocumulus clouds with light drizzle are observed, along with the highest cloud droplet number concentration (Nd), surface cloud condensation nuclei (CCN) and surface aerosol levels. The post‐trough regime features more broken or multi‐layer stratocumulus clouds with slightly higher surface rain rate, and lower Nd and surface CCN levels. The weak trough regime is characterized by the deepest MBL clouds, primarily cumulus and broken stratocumulus clouds, with the strongest surface rain rate and the lowest Nd, surface CCN and surface aerosol levels, indicating strong wet scavenging. These findings highlight the importance of considering the covariation of cloud and aerosol properties driven by large‐scale regimes when assessing aerosol indirect effects using observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A Novel Thermodynamical Predictor of Tropical High‐Cloud Area Coverage: Estimated Anvil‐Outflow Stability

The stability at the lapse-rate tropopause (LRT) was previously defined as the upper-tropospheric stability (UTS and SUT) but might underestimate the stability control on high-level clouds, since there is little direct connection between convective processes and the thermal stratification at the LRT. Here, a novel estimated anvil-outflow stability (EAS) based on the minimum stability in the upper troposphere is proposed. At the Manus site, the results show that small and large values of the LRT-based UTS and SUT both correspond to strong divergence and frequent occurrence of high-level ice clouds. In contrast, smaller EAS corresponds to stronger convective outflows to produce likely more high-level ice clouds, with a strong negative correlation. In tropics, EAS better explains the geographic distribution of high-level cloud coverage (HCC) and its temporal variations than UTS and SUT. With a strong linear correlation, EAS is likely a simple useful predictor of HCC.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Aerosol Influences on Cloud Water: Insights From ARM EPCAPE Observations With Explainable Machine Learning

This study employs an explainable machine learning (ML) framework (XGBoost-SHapley Additive exPlanations analysis) to investigate controlling factors on cloud liquid water path (LWP) using EPCAPE observations near the California coast. Aerosols are found to be the dominant factor explaining LWP variability, surpassing meteorological factors (MFs). By isolating aerosol effects from meteorological influences, the ML reveals a negative linear relationship between LWP and cloud droplet number concentration (Nd) in log space, likely driven by entrainment drying via evaporation-entrainment feedback. This aligns with the negative regime of the inverted-V relationship reported in previous studies, while no positive LWP responses are found due to a limited number of precipitating cases in EPCAPE. Furthermore, the sensitivity of LWP to Nd shows a non-linear dependence on MFs like moisture contrast between surface and free troposphere and lower-tropospheric stability. This occurs due to the interplay between the MFs' direct effects on entrainment drying and indirect effects through LWP adjustments.

EPCAPE observations↗
DOE OSTI2025

The Efficiency of Water Vapor on Top‐of‐Atmosphere Radiation

Earth's climate sensitivity is greatly affected by the compensation between temperature feedback and water vapor (WV) feedback. Using abrupt 4xCO 2 experiments, we show that the global-mean WV feedback is nearly a linear function of the temperature feedback, the slope of which is explained by the longwave radiative efficiency of WV (ϵ). Although ϵ remains constant across models in the global mean, it exhibits substantial spatial variations and is particularly weak in Antarctica, where near-surface inversions decouple the surface from the free troposphere. We introduce a surface–free troposphere temperature difference (SFTD) metric, showing that positive SFTD (e.g., high lifting condensation level) amplifies ϵ, while negative SFTD (e.g., strong surface inversion) suppresses it. These findings provide a clear explanation of how local climate conditions modulate the radiative compensation between temperature and WV feedbacks.

climate feedback↗
DOE OSTI2025

Physics or Knob‐Tuning? Tropical Anvil Peak Is Captured by GCMs

Tropical anvil clouds peak near 200 hPa and significantly impact Earth's climate, yet its physical realism in coarse‐resolution General Circulation Models (GCMs) remains debated. We examine anvil cloud formation by performing simulations with a GCM with a hierarchy of cloud fraction schemes ranging from a complex prognostic Tiedtke scheme to a simple binary scheme. All schemes consistently reproduce the anvil peak. The robust anvil peak arises because extremely cold temperatures at the upper troposphere facilitate frequent saturation events, producing clouds that disproportionately influence mean cloud fraction. Sensitivity experiments with enhanced evaporation of cloud condensate unexpectedly show increased anvil coverage, highlighting how slight evaporative moistening reinforces local saturation in cold upper‐tropospheric conditions. These results demonstrate that the tropical anvil cloud peak emerges from fundamental thermodynamic constraints, rather than specific cloud fraction parameterization choices.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Limited surface impacts of the January 2021 sudden stratospheric warming

Subseasonal weather prediction can reduce economic disruption and loss of life, especially during “windows of opportunity” when noteworthy events in the Earth system are followed by characteristic weather patterns. Sudden stratospheric warmings (SSWs), breakdowns of the winter stratospheric polar vortex, are one such event. They often precede warm temperatures in Northern Canada and cold, stormy weather throughout Europe and the United States - including the most recent SSW on January 5th, 2021. Here we assess the drivers of surface weather in the weeks following the SSW through initial condition “scrambling” experiments using the real-time CESM2(WACCM6) Earth system prediction framework. We find that the SSW itself had a limited impact, and that stratospheric polar vortex stretching and wave reflection had no discernible contribution to the record cold in North America in February. Instead, the tropospheric circulation and bidirectional coupling between the troposphere and stratosphere were dominant contributors to variability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Aircraft observations in a tropical supercluster over the equatorial Indian Ocean during MISO-BOB field campaign

The Monsoon Intra-Seasonal Oscillations in the Bay of Bengal (MISO-BOB) field campaign was conducted in the Indian Ocean during the 2018 and 2019 summer monsoon seasons. WC-130J aircraft of the 53rd Weather Reconnaissance Squadron of the US Air Force participated in the campaign in June 2018. The dropsonde observations across a tropical supercluster showed zonal wind variations in association with the structure of the convectively coupled Kelvin wave (CCKW). Within the supercluster, easterlies (westerlies) were observed in the upper (lower) troposphere; this transformation occurred just below the 0° C level. The cold pool had an easterly component throughout, and it was coldest (by 2.5° C) at the center of the supercluster, deepest (~ 1000 m) at its rear/western end, and shallowest (~ 300 m) at the front/eastern end. The level of free convection (LFC) at the front end was at 897 m altitude. At the eastern flank of the supercluster, zonal convergence in the lower troposphere occurred between 500-1500 m levels above the surface between the westerlies within the supercluster and opposing ambient easterlies. Thus, the uplifting of conditionally unstable air parcels above LFC to the east of the supercluster was likely to occur due to this convergence rather than the cold pool influence. Conversely, the western flank of the supercluster had low-level zonal divergence. These observations support the notion of ‘self-similarity’ among the mesoscale convective systems and large-scale waves.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Primary photodissociation mechanisms of pyruvic acid on S 1 : observation of methylhydroxycarbene and its chemical reaction in the gas phase

Pyruvic acid, a representative alpha-keto carboxylic acid, is one of the few organic molecules destroyed in the troposphere by solar radiation rather than by reactions with free radicals. To date, only its stable final products were identified, often with contribution from secondary chemistry, making it difficult to elucidate photodissociation mechanisms following excitation to the lowest singlet excited-state (S 1 ) and the role of the internal hydrogen bond in the most-stable Tc conformer. Using multiplexed photoionization mass spectrometry we report the first direct experimental evidence, via the observation of singlet methylhydroxycarbene (MHC) following 351 nm excitation, supporting the decarboxylation mechanism previously proposed. Decarboxylation to MHC + CO 2 represents 97–100% of product branching at 351 nm. We observe vinyl alcohol and acetaldehyde, which we attribute to isomerization of MHC. We also observe a 3 ± 2% yield of the Norrish Type I photoproducts CH 3 CO + DOCO, but only from d 1 -pyruvic acid. At 4 Torr pressure, we measure a photodissociation quantum yield of $1.0^{+0}_{–0.4}$, consistent with IUPAC recommendations. However, our measured product branching fractions disagree with IUPAC. In light of previous calculations, these results support a mechanism in which hydrogen transfer on the S 1 excited state occurs at least partially by tunneling, in competition with intersystem crossing to the T 1 state. Here, we present the first evidence of a bimolecular reaction of MHC in the gas phase, where MHC reacts with pyruvic acid to produce a C 4 H 8 O 2 product. This observation implies that some MHC produced from pyruvic acid in Earth's troposphere will be stabilized and participate in chemical reactions with O 2 and H 2 O, and should be considered in atmospheric modeling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Chemical imaging of individual stratospheric particles sampled over North America

The increasing size, severity, and frequency of wildfires have led to dramatic increases in particulate matter concentrations in the troposphere. Severe wildfires can generate intense convective systems capable of transporting large quantities of biomass burning organic aerosols (BBOA) to the upper troposphere and lower stratosphere (UTLS). Chemically complex organic matter and light-absorbing carbonaceous material is introduced into stratospheric regions that were historically isolated from direct surface emissions. In this study, stratospheric particles were sampled over North America during the Dynamics and Chemistry of the Summer Stratosphere (DCOTSS) campaign, an aircraft-based research project designed to characterize convective perturbation in the UTLS. Particle samples collected from six research flights during summer 2022 were analyzed using Computer-Controlled Scanning Electron Microscopy and Scanning Transmission X-ray Microscopy to investigate particle size distributions, morphology, chemical composition, and mixing state of stratospheric particles along transects across the continental United States and adjacent Pacific Ocean airspace. Analysis revealed that all sampled particles contained detectable levels of carbon, with most exhibiting organic volume fractions of 0.37 ± 0.20. Notably, about 5% of the particles also contained soot inclusions, which indicates the presence of refractory black carbon transported to stratospheric altitudes and provides direct evidence of wildfire-derived black carbon reaching the UTLS. Typical particle morphology exhibits organic shells over soot and inorganic cores and suggests secondary processing and aging of BBOA during transport to and within the UTLS. These findings provide compelling evidence that wildfire emissions play a critical role in affecting the long-term composition and radiative properties of stratospheric particles.

Sharpe, Steven [Purdue Univ., West Lafayette, IN (↗
DOE OSTI2020

Direct kinetic measurements and theoretical predictions of an isoprene-derived Criegee intermediate

Isoprene has the highest emission into Earth’s atmosphere of any nonmethane hydrocarbon. Atmospheric processing of alkenes, including isoprene, via ozonolysis leads to the formation of zwitterionic reactive intermediates, known as Criegee intermediates (CIs). Direct studies have revealed that reactions involving simple CIs can significantly impact the tropospheric oxidizing capacity, enhance particulate formation, and degrade local air quality. Methyl vinyl ketone oxide (MVK-oxide) is a four-carbon, asymmetric, resonance-stabilized CI, produced with 21 to 23% yield from isoprene ozonolysis, yet its reactivity has not been directly studied. We present direct kinetic measurements of MVK-oxide reactions with key atmospheric species using absorption spectroscopy. Direct UV-Vis absorption spectra from two independent flow cell experiments overlap with the molecular beam UV-Vis-depletion spectra reported recently [M. F. Vansco, B. Marchetti, M. I. Lester, J. Chem. Phys. 149, 44309 (2018)] but suggest different conformer distributions under jet-cooled and thermal conditions. Comparison of the experimental lifetime herein with theory indicates only the syn -conformers are observed; anti -conformers are calculated to be removed much more rapidly via unimolecular decay. We observe experimentally and predict theoretically fast reaction of syn -MVK-oxide with SO 2 and formic acid, similar to smaller alkyl-substituted CIs, and by contrast, slow removal in the presence of water. We determine products through complementary multiplexed photoionization mass spectrometry, observing SO 3 and identifying organic hydroperoxide formation from reaction with SO 2 and formic acid, respectively. The tropospheric implications of these reactions are evaluated using a global chemistry and transport model.

09 BIOMASS FUELS↗
DOE OSTI2021

Transport and chemistry of isoprene and its oxidation products in deep convective clouds

Deep convective clouds can transport trace gases from the planetary boundary layer into the upper troposphere where subsequent chemistry may impact aerosol particle formation and growth. In this modelling study, we investigate processes that affect isoprene and its oxidation products injected into the upper troposphere by an isolated deep convective cloud in the Amazon. We run a photochemical box model with coupled cloud microphysics along hundreds of individual air parcel trajectories sampled from a cloud-resolving model simulation of a convective event. The box model simulates gas-phase chemical reactions, gas scavenging by liquid and ice hydrometeors, and turbulent dilution inside a deep convective cloud. The results illustrate the potential importance of gas uptake to anvil ice in regulating the intensity of the isoprene oxidation and associated low volatility organic vapour concentrations in the outflow. Isoprene transport and fate also depends on the abundance of lightning-generated nitrogen oxide radicals (NO x = NO + NO 2 ). If gas uptake on ice is efficient and lightning activity is low, around 30% of the boundary layer isoprene will survive to the cloud outflow after approximately one hour of transport, while all the low volatile oxidation products will be scavenged by the cloud hydrometeors. If lightning NO x is abundant and gas uptake by ice is inefficient, then all isoprene will be oxidised during transport or in the immediate outflow region, while several low volatility isoprene oxidation products will have elevated concentrations in the cloud outflow. Reducing uncertainties associated with the uptake of vapours on ice hydrometeors, especially HO 2 and oxygenated organics, is essential to improve predictions of isoprene and its oxidation products in deep convective outflows and their potential contribution to new particle formation and growth.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Cyclone Fani: the tug-of-war between regional warming and anthropogenic aerosol effects

Before Cyclone Amphan took place in 2020, Cyclone Fani (May 2019) is the strongest pre-monsoon cyclone in the Bay of Bengal (BOB) since 1991, killing 90 people in eastern India and Bangladesh while causing US$1.81 billion of damages. Fani developed during a period of high concentration of anthropogenic aerosols in the BOB with abnormally high sea surface temperature (SST), thereby presenting an opportunity to understand the compound effects of atmospheric aerosols and regional climate warming on a tropical cyclone. A quantitative attribution analysis was conducted using the Weather Research and Forecasting model with chemistry (WRF-Chem) run at the convection-permitting (4 km) grid spacing, accompanied by an ensemble of coarser-resolution simulations to quantify the uncertainty. The removal of post-1990 trends in the tropospheric variables and SST from WRF-Chem's initial conditions (IC) and boundary conditions (BC, including the lateral and lower boundary conditions) resulted in a reduction of cyclone precipitation by about 51% during the 5 d of April 28-May 2. The removal of tropospheric warming shows approximately twice as strong an effect on Fani (39% reduction in precipitation) as that of SST warming (22% reduction). When aerosol's direct and indirect effects were removed from the simulations, i.e., no aerosol influence on radiation and cloud microphysics, Fani initially strengthened but later weakened, as measured by geopotential height and precipitation amounts. These results suggest that aerosol and its interaction with the atmosphere acted to mitigate the strengthening effect of anthropogenic warming on Fani, but was not strong enough to entirely counteract it. Although the ensemble of coarser simulations appears to overestimate Cyclone Fani in terms of precipitation, the direction of the effects is in agreement with that obtained from the 4 km simulations. Given the increasing anthropogenic aerosols in the BOB, future attribution studies using more sophisticated dynamical aerosol models on BOB tropical cyclones are urged.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

CO 2 fertilization of crops offsets yield losses due to future surface ozone damage and climate change

Tropospheric ozone (O 3 ) is harmful to plant productivity and negatively impacts crop yields. O3 concentrations are projected to decrease globally in the optimistic Representative Concentration Pathway of 2.6 W m –2 (RCP2.6) but increase globally following the high-emission scenario under the RCP8.5, with substantial implications for global food security. The damaging effect of O 3 on future crop yield is affected by CO 2 fertilization and climate change, and their interactions for RCP scenarios have yet to be quantified. In this study, we used the Joint UK Land Environment Simulator modified to include crops (JULES-crop) to quantify the impacts, and relative importance of present-day and future O3, CO 2 concentration and meteorology on crop production at the regional scale until 2100 following RCP2.6 and RCP8.5 scenarios. We focus on eight major crop-producing regions that cover the production of wheat, soybean, maize, and rice. Our results show that CO 2 alone has the largest effect on regional yields, followed by climate and O 3 . However, the CO 2 fertilization effect is offset by the negative impact of tropospheric O 3 in regions with high O 3 concentrations, such as South Asia and China. Simulated crop yields in 2050 were compared with Food and Agriculture Organisation (FAO) statistics to investigate the differences between a socioeconomic and a biophysical process-based approach. Results showed that FAO estimates are closer to our JULES-crop RCP8.5 scenario. This study demonstrates that air pollution could be the biggest threat to future food production and highlights an urgent policy need to mitigate the threat of climate change and O 3 pollution on food security.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Regional infrasonic observations from surface explosions – influence of atmospheric variations and realistic terrain

SUMMARY A pair of 1-ton, conventional surface explosions were conducted at the Nevada National Security Site in the fall of 2020 producing seismoacoustic signatures observable hundreds of kilometres from the source location. Regional infrasonic observations include tropospheric ducting at large distances to the south, a wide stratospheric waveguide with signals observed more than 700 km to the east, and anomalous arrivals in the stratospheric shadow zone. Notable differences in propagation between the events are identified despite the explosions being conducted just two days apart due to a sharp temporal shift in the tropospheric winds as well as structural changes in the stratospheric winds. Propagation simulations of the two events have been completed using a combination of ray tracing and parabolic equation (PE) methods. Simulations have been conducted to quantify the impact of the temporal variations in the atmosphere as well as the influence of terrain on propagation. Temporal variations in reflection locations are found to produce notable changes in downrange propagation due to spatially varying terrain features. Finite frequency effects modelled by the PE are found to predict ensonification not included in corresponding 2D ray tracing simulations. Notable variations in predicted signal amplitude are found due to focusing by along-path and cross-path terrain gradients; though, the later of these is only modelled using fully 3D ray tracing analysis.

58 GEOSCIENCES↗
DOE OSTI2022

MSU Disentanglement Analysis Software

This software is used to disentangle the forced-versus-unforced components of tropospheric temperature change over the satellite era (after 1979) using maps of surface temperature change as a predictor. In general, the software assembles training datasets (from pre-computed surface temperature trend maps and domain averaged tropospheric warming rates), trains statistical/machine learning (ML) algorithms, applies the trained statistical/ML model to climate model data and observations, and then saves the results. A leave-one-out approach is used in which the statistical/ML models are iteratively trained on (N- 1) climate model and then applied to the remaining climate model (and observations). Each model includes a large ensemble (i.e., >10) of model simulations. The software relies on scikit-learn ridge regression, PLS regression, and neural network algorithms.

Po-Chedley, StephenD↗
DOE OSTI2021

Impacts of Varying Concentrations of Cloud Condensation Nuclei on Deep Convective Cloud Updrafts—A Multimodel Assessment

Abstract This study presents results from a model intercomparison project, focusing on the range of responses in deep convective cloud updrafts to varying cloud condensation nuclei (CCN) concentrations among seven state-of-the-art cloud-resolving models. Simulations of scattered convective clouds near Houston, Texas, are conducted, after being initialized with both relatively low and high CCN concentrations. Deep convective updrafts are identified, and trends in the updraft intensity and frequency are assessed. The factors contributing to the vertical velocity tendencies are examined to identify the physical processes associated with the CCN-induced updraft changes. The models show several consistent trends. In general, the changes between the High-CCN and Low-CCN simulations in updraft magnitudes throughout the depth of the troposphere are within 15% for all of the models. All models produce stronger (~+5%–15%) mean updrafts from ~4–7 km above ground level (AGL) in the High-CCN simulations, followed by a waning response up to ~8 km AGL in most of the models. Thermal buoyancy was more sensitive than condensate loading to varying CCN concentrations in most of the models and more impactful in the mean updraft responses. However, there are also differences between the models. The change in the amount of deep convective updrafts varies significantly. Furthermore, approximately half the models demonstrate neutral-to-weaker (~−5% to 0%) updrafts above ~8 km AGL, while the other models show stronger (~+10%) updrafts in the High-CCN simulations. The combination of the CCN-induced impacts on the buoyancy and vertical perturbation pressure gradient terms better explains these middle- and upper-tropospheric updraft trends than the buoyancy terms alone.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A Theory for the Response of Tropical Moist Convection to Mechanical Orographic Forcing

Abstract Spatial patterns of tropical rainfall are strongly influenced by mountains. Although theories for precipitation induced by convectively stable upslope ascent exist for the midlatitudes, these do not represent the interaction of moist convection with orographic forcing. Here, we present a theory for convective precipitation produced by the mechanical interaction of a tropical ridge with a basic-state horizontal wind. Deviations from this basic state are represented as the sum of a “dry” perturbation, due to the stationary orographic gravity wave, and a “moist” perturbation that carries the convective response. The moist component dynamics are subject to the weak temperature gradient approximation; they are forced by the dry mode’s influence on lower-tropospheric moisture and temperature. Analytical solutions provide estimates of the precipitation distribution, including peak precipitation, upstream extent, and rain shadow extent. The theory can be used with several degrees of complexity depending on the technique used to compute the dry mode, which can be drawn from linear mountain wave theory or full numerical simulations. To evaluate the theory, we use a set of convection-permitting simulations with a flow-perpendicular ridge in a long channel. The theory makes a good prediction for the cross-slope precipitation profile, indicating that the organization of convective rain by orography can be quantitatively understood by considering the effect of stationary orographic gravity waves on a lower-tropospheric convective quasi-equilibrium state.

54 ENVIRONMENTAL SCIENCES↗