Search NASA⌕ Search

SEARCH · Search NASA

Results for “Climate feedbacks”

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 127 records · Page 7

Cloud Radiation Forcings and Feedbacks: General Circulation Model Tests and Observational Validation

Using an atmospheric general circulation model (the National Center for Atmospheric Research Community Climate Model: CCM2), the effects on climate sensitivity of several different cloud radiation parameterizations have been investigated. In addition to the original cloud radiation scheme of CCM2, four parameterizations incorporating prognostic cloud water were tested: one version with prescribed cloud radiative properties and three other versions with interactive cloud radiative properties. The authors' numerical experiments employ perpetual July integrations driven by globally constant sea surface temperature forcings of two degrees, both positive and negative. A diagnostic radiation calculation has been applied to investigate the partial contributions of high, middle, and low cloud to the total cloud radiative forcing, as well as the contributions of water vapor, temperature, and cloud to the net climate feedback. The high cloud net radiative forcing is positive, and the middle and low cloud net radiative forcings are negative. The total net cloud forcing is negative in all of the model versions. The effect of interactive cloud radiative properties on global climate sensitivity is significant. The net cloud radiative feedbacks consist of quite different shortwave and longwave components between the schemes with interactive cloud radiative properties and the schemes with specified properties. The increase in cloud water content in the warmer climate leads to optically thicker middle- and low-level clouds and in turn to negative shortwave feedbacks for the interactive radiative schemes, while the decrease in cloud amount simply produces a positive shortwave feedback for the schemes with a specified cloud water path. For the longwave feedbacks, the decrease in high effective cloudiness for the schemes without interactive radiative properties leads to a negative feedback, while for the other cases, the longwave feedback is positive. These cloud radiation parameterizations are empirically validated by using a single-column diagnostic model. together with measurements from the Atmospheric Radiation Measurement program and from the Tropical Ocean Global Atmosphere Combined Ocean-Atmosphere Response Experiment. The inclusion of prognostic cloud water produces a notable improvement in the realism of the parameterizations, as judged by these observations. Furthermore, the observational evidence suggests that deriving cloud radiative properties from cloud water content and microphysical characteristics is a promising route to further improvement.

Lee,Wan-Ho↗

Long-Term Monitoring of Global Climate Forcings and Feedbacks

A workshop on Long-Term Monitoring of Global Climate Forcings and Feedbacks was held February 3-4, 1992, at NASA's Goddard Institute for Space Studies to discuss the measurements required to interpret long-term global temperature changes, to critique the proposed contributions of a series of small satellites (Climsat), and to identify needed complementary monitoring. The workshop concluded that long-term (several decades) of continuous monitoring of the major climate forcings and feedbacks is essential for understanding long-term climate change.

Hansen, J.↗

Arctic Cloud‐Base Ice Precipitation Properties Retrieved Using Bayesian Inference

Abstract Cloud‐climate feedbacks are still the greatest source of uncertainty in current climate projections. Arctic clouds, which are predominantly stratiform and supercooled, often long‐lived, and nearly continuously precipitate ice particles, contribute roughly 10% of the uncertainty attributed to the global cloud feedback. This Arctic cloud uncertainty is driven by incomplete observational and theoretical knowledge required to estimate and explain the state and active processes occurring in those clouds. A focus on ice precipitation properties at Arctic cloud base rather than the surface deconfounds the product of cloud condensate sink processes from the influence of the atmospheric thermodynamic state below cloud base, rendering cloud‐base properties a more appealing target for inference and evaluation of model simulations. Here I describe an inverse model for the estimation of cloud base ice precipitation properties over Utqiagvik, North Slope of Alaska, using the synthesis of ground‐based radar and lidar measurements. By leveraging a Markov Chain Monte Carlo algorithm as the core of the inverse model, a wide range of particle size distributions are sampled, and different combinations of ice habit models are examined, both of which are typically fixed in other retrieval methods. Results show intriguing links between different cloud base thermodynamic and ice precipitation properties. Apparent ice number concentration enhancements at temperatures of −5 and −15°C suggest possible secondary ice production (SIP). The analysis alludes to an overestimation of SIP occurrence and intensity, especially in studies relying only on radar or lidar measurements. Finally, reflectivity‐dependent ice precipitation rate and ice water content parameterizations are presented.

54 ENVIRONMENTAL SCIENCES↗

Arctic Mixed-Phase Cloud Base Ice Precipitation Properties Over the NSA Site

Cloud-climate feedbacks are still the greatest source of uncertainty in current climate projections. Arctic clouds, which are predominantly stratiform and supercooled, often long-lived, and nearly continuously precipitate ice particles, contribute roughly 10% of the uncertainty attributed to the global cloud feedback. This arctic cloud uncertainty is driven by incomplete observational and theoretical knowledge required to estimate and explain the state and active processes occurring in those clouds. A focus on ice precipitation properties at arctic cloud base rather than the surface deconfounds the product of cloud condensate sink processes from the influence of the atmospheric thermodynamic state below cloud base, rendering cloud-base properties a more appealing target for inference and evaluation of model simulations. This data set provides more than 1800 samples of cloud-base ice precipitation properties over Utqiagvik, North Slope of Alaska, all of which were retrieved using the synthesis of ARM radar and lidar measurements. The retrieved ice precipitation variables in this data set include, among others, the ice number concentration, water content, PSD parameters, precipitation rate, mass-weighted fall velocity, vertical air motion, and effective radius, all of which are highly valuable for model evaluation and a general understanding of polar cloud sink processes. Each variable sample includes its mean value and associated uncertainty. Additional variables based on ARM measurements (liquid layer statistics, etc.) are included in this data set. The retrieval algorithm and analysis of this data set are described in Silber (JGR, 2023, https://doi.org/10.1029/2022JD038202).

54 ENVIRONMENTAL SCIENCES↗

Arctic Mixed-Phase Cloud Base Ice Precipitation Properties During the M-PACE Field Campaign

Cloud-climate feedbacks are still the greatest source of uncertainty in current climate projections. Arctic clouds, which are predominantly stratiform and supercooled, often long-lived, and nearly-continuously precipitate ice particles, contribute roughly 10% of the uncertainty attributed to the global cloud feedback. This Arctic cloud uncertainty is driven by incomplete observational and theoretical knowledge required to estimate and explain the state and active processes occurring in those clouds. A focus on ice precipitation properties at Arctic cloud base rather than the surface deconfounds the product of cloud condensate sink processes from the influence of the atmospheric thermodynamic state below cloud base, rendering cloud-base properties a more appealing target for inference and evaluation of model simulations. This dataset provides a set of 25 samples from the M-PACE field campaign, all of which were retrieved using the synthesis of ARM radar and lidar measurements. The retrieved ice precipitation variables in this dataset include, among others, the ice number concentration, water content, PSD parameters, precipitation rate, mass-weighted fall velocity, vertical air motion, and effective radius, all of which are highly valuable for model evaluation and a general understanding of polar cloud sink processes. Each variable sample includes its mean value and associated uncertainty. Additional variables based on ARM measurements (liquid layer statistics, etc.) are included in this dataset as well. The retrieval algorithm and analysis of this dataset are described in Silber (JGR, 2023, https://doi.org/10.1029/2022JD038202).

54 ENVIRONMENTAL SCIENCES↗

Dispersal and fire limit Arctic shrub expansion

Abstract Arctic shrub expansion alters carbon budgets, albedo, and warming rates in high latitudes but remains challenging to predict due to unclear underlying controls. Observational studies and models typically use relationships between observed shrub presence and current environmental suitability (bioclimate and topography) to predict shrub expansion, while omitting shrub demographic processes and non-stationary response to changing climate. Here, we use high-resolution satellite imagery across Alaska and western Canada to show that observed shrub expansion has not been controlled by environmental suitability during 1984–2014, but can only be explained by considering seed dispersal and fire. These findings provide the impetus for better observations of recruitment and for incorporating currently underrepresented processes of seed dispersal and fire in land models to project shrub expansion and climate feedbacks. Integrating these dynamic processes with projected fire extent and climate, we estimate shrubs will expand into 25% of the non-shrub tundra by 2100, in contrast to 39% predicted based on increasing environmental suitability alone. Thus, using environmental suitability alone likely overestimates and misrepresents shrub expansion pattern and its associated carbon sink.

54 ENVIRONMENTAL SCIENCES↗

Clouds and Convective Self-Aggregation in a Multimodel Ensemble of Radiative-Convective Equilibrium Simulations

The Radiative-Convective Equilibrium Model Intercomparison Project (RCEMIP) is an intercomparison of multiple types of numerical models configured in radiative-convective equilibrium (RCE). RCE is an idealization of the tropical atmosphere that has long been used to study basic questions in climate science. Here, we employ RCE to investigate the role that clouds and convective activity play in determining cloud feedbacks, climate sensitivity, the state of convective aggregation, and the equilibrium climate. RCEMIP is unique among intercomparisons in its inclusion of a wide range of model types, including atmospheric general circulation models (GCMs), single column models (SCMs), cloud-resolving models (CRMs), large eddy simulations (LES), and global cloud-resolving models (GCRMs). The first results are presented from the RCEMIP ensemble of more than 30 models. While there are large differences across the RCEMIP ensemble in the representation of mean profiles of temperature, humidity, and cloudiness, in a majority of models anvil clouds rise, warm, and decrease in area coverage in response to an increase in sea surface temperature (SST). Nearly all models exhibit self-aggregation in large domains and agree that self-aggregation acts to dry and warm the troposphere, reduce high cloudiness, and increase cooling to space. The degree of self-aggregation exhibits no clear tendency with warming. There is a wide range of climate sensitivities, but models with parameterized convection tend to have lower climate sensitivities than models with explicit convection. In models with parameterized convection, aggregated simulations have lower climate sensitivities than unaggregated simulations.

54 ENVIRONMENTAL SCIENCES↗

A data-driven global soil heterotrophic respiration dataset and the drivers of its inter-annual variability

Soil heterotrophic respiration (SHR), one of the primary carbon fluxes from terrestrial ecosystems to the atmosphere, is important for carbon-climate feedbacks because of its sensitivity to available litter and soil carbon, climatic conditions, and nutrient availability. However, until recently limited SHR data were available, and most published global SHR estimates have either a short time span, coarse spatial resolution, or reply on overly-simple model formulations. To better understand and quantify the global distribution of SHR and its sensitivity to climate variability, we produced a new global SHR dataset using Random Forest algorithms, up-scaling 455 point data from the Global Soil Respiration Database (SRDB 4.0) with gridded fields of climatic, edaphic and productivity as explanatory variables. We estimated a global total SHR of 46.8 Pg C yr-1 over 1985-2013 (95% confidence interval: 38.6-56.3 Pg C yr-1), with a significant increasing trend of 0.03 Pg C yr-2 during this period. We found that the choice of soil moisture datasets contributes more to the difference among these data-driven SHR members rather than that of productivity, temperature and precipitation data sources. We also analyzed the influence of climatic variables on the inter-annual variability (IAV) of our SHR product. Water availability was the dominant driver of IAV at global scales, although the inferred sensitivity depends on the choice of the soil moisture gridded dataset. At the ecosystem scale, temperature strongly controls the IAV of SHR in tropical forests, while water availability dominates in extra-tropical forest and semi-arid regions. Our machine-learning gridded SHR dataset and outputs from process-based land surface models (TRENDYv6) show agreement for a strong association between water variability and SHR IAV at the global scale, but the two approaches lead to different temporal trend globally and different controlling variables for IAV at the ecosystem scale. Our study provides evidence for the pervasive and important role of water availability in driving SHR, indicating both a direct effect limiting decomposition rates and an indirect effect through the amount of fresh organic matter made available to SHR from productivity. In consideration of potential limitations and uncertainties remaining in our data-driven SHR datasets, we call for a more scientifically designed observation network for SHR, more observation data compilation, and increased use of deep learning methods making maximum use of observation data in hand. This will benefit process-based models, and improve our understanding of SHR response to future anomalous environmental conditions.

Yao, Yitong↗

Radiative, Hydrologic, and Circulation Responses to Warming in Cess‐Potter Simulations Using the Global 3.25‐km SCREAM

Using the global 3.25-km Simple Cloud Resolving E3SM Atmosphere Model (SCREAM 3 km), a pair of 13-month Cess-Potter simulations are performed to quantify the radiative feedbacks and the hydrologic and circulation responses to warming. Large-scale aspects of SCREAM 3 km's top-of-atmosphere radiative fluxes, precipitation rates, and circulations are in good agreement with observations and reanalysis, with notable differences, including a drier lower free-troposphere in the Tropics, reduced precipitation and humidity over the Tropical West Pacific, and poleward shifted Southern Hemisphere midlatitude jet. In response to warming, SCREAM 3 km predicts a total radiative feedback within the top 15% of the CMIP5 and CMIP6 models, which puts it substantially higher than the feedback reported by other kilometer-scale models. SCREAM 3 km's high radiative feedback stems from a strongly positive shortwave cloud feedback, most prominent over the mid- and high-latitudes. SCREAM 3 km's high precipitation response also puts it among the highest of CMIP models, whereas its circulation response is within the spread of CMIP models. An ensemble of five perturbed initial condition Cess-Potter simulations with a 12 km version of SCREAM (SCREAM 12 km) is performed to characterize uncertainty and resolution sensitivity. It suggests that the uncertainty from analyzing a pair of 1-year simulations is small compared to the inter-model spread in feedbacks and precipitation response. SCREAM 12 km also produces a strong precipitation response to warming but a much lower cloud feedback and total radiative feedback. The results from these experiments suggest that the spread in climate feedbacks will likely persist in the next generation of kilometer-scale models.

54 ENVIRONMENTAL SCIENCES↗

Regional environmental controllers influence continental scale soil carbon stocks and future carbon dynamics

Abstract Understanding the influence of environmental factors on soil organic carbon (SOC) is critical for quantifying and reducing the uncertainty in carbon climate feedback projections under changing environmental conditions. We explored the effect of climatic variables, land cover types, topographic attributes, soil types and bedrock geology on SOC stocks of top 1 m depth across conterminous United States (US) ecoregions. Using 4559 soil profile observations and high-resolution data of environmental factors, we identified dominant environmental controllers of SOC stocks in 21 US ecoregions using geographically weighted regression. We used projected climatic data of SSP126 and SSP585 scenarios from GFDL-ESM 4 Earth System Model of Coupled Model Intercomparison Project phase 6 to predict SOC stock changes across continental US between 2030 and 2100. Both baseline and predicted changes in SOC stocks were compared with SOC stocks represented in GFDL-ESM4 projections. Among 56 environmental predictors, we found 12 as dominant controllers across all ecoregions. The adjusted geospatial model with the 12 environmental controllers showed an R 2 of 0.48 in testing dataset. Higher precipitation and lower temperatures were associated with higher levels of SOC stocks in majority of ecoregions. Changes in land cover types (vegetation properties) was important in drier ecosystem as North American deserts, whereas soil types and topography were more important in American prairies. Wetlands of the Everglades was highly sensitive to projected temperature changes. The SOC stocks did not change under SSP126 until 2100, however SOC stocks decreased up to 21% under SSP585. Our results, based on environmental controllers of SOC stocks, help to predict impacts of changing environmental conditions on SOC stocks more reliably and may reduce uncertainties found in both, geospatial and Earth System Models. In addition, the description of different environmental controllers for US ecoregions can help to describe the scope and importance of global and local models.

58 GEOSCIENCES↗

Process Drivers, Inter-model Spread, and a Path Forward: What Do We Know About Arctic Amplification?

Over the recent decades, rapid changes have been observed in the Arctic resulting from enhanced surface warming in the Arctic relative to the rest of the globe, a phenomenon called Arctic Amplification. While surface-based amplification of Arctic warming in response to increased CO2 has been established. The driving mechanisms, relevant processes, and the relative roles of local and remote mechanisms to Arctic Amplification remain unclear. Research over the last decade has made it clear that coupling between the atmosphere, sea ice, and ocean play a key role in Arctic Amplification. This presentation summarizes recent research on the fundamental role that atmosphere-sea ice-ocean interactions and energy exchanges play in driving Arctic Amplification, how these interactions manifest as radiative climate feedbacks, and how these processes influence the differences in the projected Arctic climate change across CMIP5 and 6. We hypothesize that models that simulate larger warming do so by generating a larger local atmospheric circulation response that more effectively redistributes energy drawn from the ocean in sea ice retreat regions Arctic-wide and contributing to a larger increase in downwelling long wave radiation. This process is also fundamentally tied to the Arctic lapse rate feedback.

Patrick Charles Taylor↗

Process drivers, inter-model spread, and a path forward: What do we know about Arctic Amplification?

Over the recent decades, rapid changes have been observed in the Arctic resulting from enhanced surface warming in the Arctic relative to the rest of the globe, a phenomenon called Arctic Amplification. While surface-based amplification of Arctic warming in response to increased CO2 has been established. The driving mechanisms, relevant processes, and the relative roles of local and remote mechanisms to Arctic Amplification remain unclear. Research over the last decade has made it clear that coupling between the atmosphere, sea ice, and ocean play a key role in Arctic Amplification. This presentation summarizes recent research on the fundamental role that atmosphere-sea ice-ocean interactions and energy exchanges play in driving Arctic Amplification, how these interactions manifest as radiative climate feedbacks, and how these processes influence the differences in the projected Arctic climate change across CMIP5 and 6. We hypothesize that models that simulate larger warming do so by generating a larger local atmospheric circulation response that more effectively redistributes energy drawn from the ocean in sea ice retreat regions Arctic-wide and contributing to a larger increase in downwelling long wave radiation. This process is also fundamentally tied to the Arctic lapse rate feedback.

Patrick C. Taylor↗

Influence of Tundra Polygon Type and Climate Variability on CO 2 and CH 4 Fluxes Near Utqiagvik, Alaska

Arctic tundra has the potential to generate significant climate feedbacks, but spatial complexity makes it difficult to quantify the impacts of climate on ecosystem-atmosphere fluxes, particularly in polygonal tundra comprising wetter and drier polygon types on the scale of tens of meters. We measured CO 2 , CH 4 , and energy fluxes using eddy covariance for 7 yr (April to November, 2013–2019) in polygonal tundra near Utqiagvik, Alaska. This period saw the earliest snowmelt, latest snow accumulation, and hottest summer on record. To estimate fluxes by polygon type, we combined a polygon classification with a flux-footprint model. Methane fluxes were highest in the summer months but were also large during freeze-up and increased with the warming trend in August–November temperatures. While CO 2 respiration had a consistent, exponential relationship with temperature, net ecosystem exchange was more variable among years. CO 2 and CH 4 exchange (June–September) ranged between -0.83 (Standard error [SE] = 0.03) and -1.32 (SE = 0.04) μmol m -2 s-1 and 13.92 (SE = 0.26)—23.42 (SE = 0.45) nmol m -2 s -1 , respectively, and varied interannually (p ≤ 0.05). The maximum-influence method effectively attributed fluxes to polygon types. Areas dominated by low-centered polygons had higher CO 2 fluxes except in 2016–2017. Methane fluxes were highest in low-centered polygons 2013–2015 and in flat-centered polygons in subsequent years, possibly due to increasing temperature and precipitation. Sensible and latent heat fluxes also varied significantly among polygon types. Accurate characterization of Arctic fluxes and their climate dependencies requires spatial disaggregation and long term observations.

54 ENVIRONMENTAL SCIENCES↗

Cloud feedback - A stabilizing effect for the early earth

The effect of variations in cloud cover, optical properties, and fractional distribution with altitude on the mean surface temperature of a model of the early earth has been investigated. In all cases examined, cloud-climate feedbacks result in temperatures greater than those in models with no cloud feedbacks. If the model of hydrospheric feedback effects is correct, then cloud feedbacks are as important to the climate as changes in solar luminosity and atmospheric composition during the earth's atmospheric evolution. In particular, the early earth need not become completely ice-covered if strong negative cloud feedbacks occur. However, until a proper understanding of cloud feedbacks is available, conclusions regarding conditions in the early atmosphere must remain in doubt.

Rossow, W. B.↗

Using Satellites to Understand Climate and Climate Change

This viewgraph presentation reviews the measurement of climate with the use of satellites. The basic greenhouse effect, Ice-albedo feedback, climate models and observations, aerosol-cloud interactions, and the Antarctic are discussed, along with the human effect on climate change.

remote sensing↗

The Interaction Between Climate Forcing and Feedbacks

A Perturbed Parameter Ensemble (PPE) with the Community Atmosphere Model version 6 (CAM6) is used to better understand the sensitivity of aerosol forcing and cloud feedbacks to changes in model processes. Aerosol forcing through aerosol-cloud interactions is mostly negative (a cooling) due to shortwave radiation, while feedbacks are positive or negative in different regions due to contrasting longwave and shortwave effects. Both forcing and feedbacks are related to the mean climate state. Higher magnitude cloud radiative effects generally mean larger magnitude net negative forcing and larger magnitude net positive feedback. Aerosol forcing is broadly related to the susceptibility of clouds to drop number. Feedbacks also related to susceptibility, but to a lesser extent and in different regions to aerosol forcing. Aerosol forcing and cloud feedbacks are anti-correlated in the CAM6 PPE such that stronger negative forcing is associated with stronger positive feedbacks. Even the processes governing forcing and feedback sensitivity in the PPE are similar. These include the warm rain formation process, ice loss processes and deep convective intensity.

54 ENVIRONMENTAL SCIENCES↗

Combining Emergent Constraints for Climate Sensitivity

A method is proposed for combining information from several emergent constraints into a probabilistic estimate for a climate sensitivity proxy Y such as equilibrium climate sensitivity (ECS). The method is based on fitting a multivariate Gaussian PDF for Y and the emergent constraints using an ensemble of global climate models (GCMs); it can be viewed as a form of multiple linear regression of Y on the constraints. The method accounts for uncertainties in sampling this multidimensional PDF with a small number of models, for observational uncertainties in the constraints, and for overconfidence about the correlation of the constraints with the climate sensitivity. Its general form (Method C) accounts for correlations between the constraints. Method C becomes less robust when some constraints are too strongly related to each other; this can be mitigated using regularization approaches such as ridge regression. An illuminating special case, Method U, neglects any correlations between constraints except through their mutual relationship to the climate proxy; it is more robust to small GCM sample size and is appealingly interpretable. These methods are applied to ECS and the climate feedback parameter using a previously published set of 11 possible emergent constraints derived from climate models in the Coupled Model Intercomparison Project (CMIP). The ±2σ posterior range of ECS for Method C with no overconfidence adjustment is 4.3 ± 0.7 K. For Method U with a large overconfidence adjustment, it is 4.0 ± 1.3 K. This study adds confidence to past findings that most constraints predict higher climate sensitivity than the CMIP mean.

54 ENVIRONMENTAL SCIENCES↗

Accurate assessment of land–atmosphere coupling in climate models requires high-frequency data output

Land–atmosphere (L–A) interactions are important for understanding convective processes, climate feedbacks, the development and perpetuation of droughts, heatwaves, pluvials, and other land-centered climate anomalies. Local L–A coupling (LoCo) metrics capture relevant L–A processes, highlighting the impact of soil and vegetation states on surface flux partitioning and the impact of surface fluxes on boundary layer (BL) growth and development and the entrainment of air above the BL. A primary goal of the Climate Process Team in the Coupling Land and Atmospheric Subgrid Parameterizations (CLASP) project is parameterizing and characterizing the impact of subgrid heterogeneity in global and regional Earth system models (ESMs) to improve the connection between land and atmospheric states and processes. A critical step in achieving that aim is the incorporation of L–A metrics, especially LoCo metrics, into climate model diagnostic process streams. However, because land–atmosphere interactions span timescales of minutes (e.g., turbulent fluxes), hours (e.g., BL growth and decay), days (e.g., soil moisture memory), and seasons (e.g., variability in behavioral regimes between soil moisture and latent heat flux), with multiple processes of interest happening in different geographic regions at different times of year, there is not a single metric that captures all the modes, means, and methods of interaction between the land and the atmosphere. And while monthly means of most of the LoCo-relevant variables are routinely saved from ESM simulations, data storage constraints typically preclude routine archival of the hourly data that would enable the calculation of all LoCo metrics. Here, we outline a reasonable data request that would allow for adequate characterization of sub-daily coupling processes between the land and the atmosphere, preserving enough sub-daily output to describe, analyze, and better understand L–A coupling in modern climate models. A secondary request involves embedding calculations within the models to determine mean properties in and above the BL to further improve characterization of model behavior. Higher-frequency model output will (i) allow for more direct comparison with observational field campaigns on process-relevant timescales, (ii) enable demonstration of inter-model spread in L–A coupling processes, and (iii) aid in targeted identification of sources of deficiencies and opportunities for improvement of the models.

54 ENVIRONMENTAL SCIENCES↗