Search NASA⌕ Search

SEARCH · Search NASA

Results for “climate index”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Integrated assessment model diagnostics: key indicators and model evolution

Integrated assessment models (IAMs) form a prime tool in informing climate mitigation strategies. Diagnostic indicators that allow to compare these models can help to describe and explain differences in model projections. This also increases transparency and comparability. Earlier, the IAM community has developed an approach to diagnose models (Kriegler et al., 2015). Here we build on this, by proposing a selected set of well-defined indicators as a community standard, similar to metrics used for other modeling communities such as climate models. These indicators are the relative abatement index (RAI), emission reduction type index (ERT), inertia timescale (IT), fossil fuel reduction (FFR), transformation index (TI) and cost per abatement value (CAV). We apply the approach to 17 IAMs, including both older version as well as their latest versions, as applied in the IPCC 6th Assessment Report (AR6). The study shows that the approach can be easily applied and allows for comparison of model versions in time. The indicators and their trends can often be explained in terms of model characteristics and changes. We show that together, the set of six indicators can provide an useful indication of the main traits of the model and can roughly indicate the general model behavior. The results also show that there is often a considerable spread across the models. Interestingly, the diagnostic values often change for different model versions, but there does not seem to be a distinct trend across the different models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Attribution of the Australian bushfire risk to anthropogenic climate change

Abstract. Disastrous bushfires during the last months of 2019 and January 2020 affected Australia, raising the question to what extent the risk of these fires was exacerbated by anthropogenic climate change. To answer the question for southeastern Australia, where fires were particularly severe, affecting people and ecosystems, we use a physically based index of fire weather, the Fire Weather Index; long-term observations of heat and drought; and 11 large ensembles of state-of-the-art climate models. We find large trends in the Fire Weather Index in the fifth-generation European Centre for Medium-Range Weather Forecasts (ECMWF) Atmospheric Reanalysis (ERA5) since 1979 and a smaller but significant increase by at least 30 % in the models. Therefore, we find that climate change has induced a higher weather-induced risk of such an extreme fire season. This trend is mainly driven by the increase of temperature extremes. In agreement with previous analyses we find that heat extremes have become more likely by at least a factor of 2 due to the long-term warming trend. However, current climate models overestimate variability and tend to underestimate the long-term trend in these extremes, so the true change in the likelihood of extreme heat could be larger, suggesting that the attribution of the increased fire weather risk is a conservative estimate. We do not find an attributable trend in either extreme annual drought or the driest month of the fire season, September–February. The observations, however, show a weak drying trend in the annual mean. For the 2019/20 season more than half of the July–December drought was driven by record excursions of the Indian Ocean Dipole and Southern Annular Mode, factors which are included in the analysis here. The study reveals the complexity of the 2019/20 bushfire event, with some but not all drivers showing an imprint of anthropogenic climate change. Finally, the study concludes with a qualitative review of various vulnerability and exposure factors that each play a role, along with the hazard in increasing or decreasing the overall impact of the bushfires.

van Oldenborgh, Geert Jan (ORCID:0000000268989535)↗
DOE OSTI2023

The effect of the vertical gradients of photosynthetic parameters on the CO 2 assimilation and transpiration of a Panamanian tropical forest

Terrestrial biosphere models (TBMs) include the representation of vertical gradients in leaf traits associated with modeling photosynthesis, respiration, and stomatal conductance. However, model assumptions associated with these gradients have not been tested in complex tropical forest canopies. Here we compared TBM representation of the vertical gradients of key leaf traits with measurements made in a tropical forest in Panama, and then quantified the impact of the observed gradients on simulated canopy scale CO 2 and water fluxes. Comparison between observed and TBM trait gradients showed divergence that impacted canopy scale simulations of water vapor and CO 2 exchange. Notably, the ratio between the dark respiration rate and the maximum carboxylation rate was lower near the ground than at the top-of-canopy, leaf-level water-use efficiency was markedly higher at the top-of-canopy, and the decrease in maximum carboxylation rate from the top-of-canopy to the ground was less than TBM assumptions. The representation of the gradients of leaf traits in TBMs is typically derived from measurements made within-individual plants, or, for some traits, assumed constant due to a lack of experimental data. Our work shows that these assumptions are not representative of the trait gradients observed in species-rich, complex tropical forests.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Soil organic matter diagenetic state informs boreal forest ecosystem feedbacks to climate change

The fate of soil organic carbon (SOC) in boreal forests is dependent on the integrative ecosystem response to climate change. For example, boreal forest productivity is often nitrogen (N) limited, and climate warming can enhance N cycling and primary productivity. However, the net effect of this feedback on the SOC reservoir and its longevity with climate change remain unclear due to difficulty in detecting small differences between large and variable carbon (C) fluxes needed to determine net changes in soil reservoirs. The diagenetic state of SOC – resulting from the physicochemical and biological transformations that alter the original biomolecular composition of detrital inputs to soil over time – is useful for tracing the net response of SOC at the timescales relevant to climate change not usually discernible from fluxes and stocks alone. Here, we test two hypotheses using a mesic boreal forest climate transect: (1) the SOC diagenetic state is maintained across this climosequence, and (2) the maintenance of the SOC diagenetic state is a consequence of coupled soil C and N cycling, signifying the role of enhanced N cycling supporting SOC inputs that maintain SOC stocks within the warmer-climate forests. Shifts in nonvascular to vascular plant inputs with climate observed in these and other boreal forests highlighted the need to carefully separate biogeochemical indicators of SOC source from those signifying diagenetic alteration. We thus evaluated and applied lignin biomarkers to assess the diagenetic alteration of SOC in these boreal forest organic soils and directly compared the lignin diagenetic state with that of soil organic nitrogen (SON) assessed through amino acid composition. The lignin diagenetic state remained constant across the climate transect, indicating a balance between the input and removal of lignin in these mesic boreal forests. When combined with previous knowledge of these forest ecosystems, including the diagenetic state of SON and direct measures of C fluxes and stocks, the results indicate a coupled increase in C and N cycling with climate warming that supports forest productivity and maintains SOC stocks. This balance could markedly shift as other factors begin to limit forest productivity (e.g., trace nutrients, water) with further climate change or affect forest nutrient allocation (e.g., forest age or compositional change). Further application of the approach presented here could be used to detect the limits of this and other ecosystem–climate feedbacks, by providing a tractable and parameterizable index of the lignin state across large spatial scales, necessary for ecosystem-scale parameterizations.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1987

Evaluating North American net primary productivity with satellite observations

An ecological model is developed to estimate annual net primary productivity (NPP) in 12 North American biomes. The model combines existing models which address canopy photosynthesis in response to light, temperature, and moisture availability, and account for respiration. Climate data, solar radiation data, and spectral vegetation index data are utilized. Estimates of NPP from the model compare well with data in the literature, but a systematic error is suspected. Difficulties encountered in specifying certain model parameters are discussed as possible sources of this error. The results of this study suggest the promise of remotely sensed measurements for macroscale evaluation and modeling of NPP.

Goward, Samuel N.↗
NASA NTRS2002

Is There Really an Intermittent Biennial Oscillation in the Great Plains Low-Level Jet Over Texas?

In the 15-year GEOS-1 reanalysis data set, a maximum of interannual variance of low- level meridional flow for the warm season (May through August) occurs over southeast Texas. This variance maximum seems to be dominated by a marked biennial oscillation that occurs only during the first 6 (or possibly 8) years of the reanalysis period (1980-85 or possibly 1980-1987) and then completely disappears by the 9th year. This biennial oscillation seems to be associated with interannual fluctuations in ground wetness, surface temperature and surface pressure gradients over Texas. The periods of drier soil lead to warmer surface temperatures, lower surface pressures, stronger pressure gradients between Texas and the Gulf of Mexico and stronger southerly winds. This intermittent biennial oscillation is also evident in corresponding fields for the the NCEP/NCAR reanalysis data set for the years 1978-1985 (and possibly from 1978- 1987) and 1995-2000, but not during other periods. There are also obvious biennial oscillations evident during these periods in U.S. Climate Division records for the Palmer Drought Severity Index (PDSI) for Texas. Month-by-month correlations of this index with certain el Nino related indices are as high as .45 for the first period and as high as .55 or .6 for the second period for some regions in Texas. The seasonal cycle of the biennial signal in the PDSI and precipitation for the first period suggest that the drought in Texas and Mexico is ended (caused) by a reversal in the sign of anomalies in precipitation rate for the fall/winter season. Analysis of tropical Pacific SST patterns shows a .5 to .75 K biennial oscillation of SSTs along the precipitation-free track to the southwest of the Mexican coast during the fall and winter months of the 1978 to 1985 period that might explain the reversal in precipitation anomalies and hence the entire intermittent biennial oscillation in ground hydrology and low-level flow.

Helfand, H. Mark↗
NASA NTRS2020

Montana Water Resources: Developing a Composite Moisture Index Utilizing NASA Earth Observations for Drought Monitoring in the Missouri River Basin

The Missouri River Basin provides irrigation water for a substantial part of the domestic agricultural sector within the United States. Drought events pose a significant threat to the economic livelihoods of dependent individuals, industries, and ecosystems (e.g. farmers, local tribes, hydroelectric power, wildlife). In 2017 alone, the Missouri River Basin experienced a severe drought that resulted in a $2.6 billion loss to the U.S. Northern Plains. In response to such events, organizations throughout the basin, such as the Montana Climate Office, have dedicated efforts for drought monitoring and communicating relevant information to local stakeholders. In an effort to aid regional decision-making capabilities, the project team partnered with the Montana Climate Office, NOAA National Weather Service (NWS) Missouri Basin River Forecast Center, and NOAA Regional Climate Services, Central Region to create a monthly Composite Moisture Index (CMI) that relies on NASA Earth observations from the Terra Moderate Resolution Imaging Spectroradiometer (MODIS) and the Soil Moisture Active Passive (SMAP) mission. From these satellites, as well as the NOAA NWS National Operational Hydrologic Remote Sensing Center’s Snow Data Assimilation System (SNODAS), our team aggregated climate datasets including soil moisture, snow cover, snow depth, and snow water equivalent to compute a CMI that indicates regional moisture conditions during the winter months. The March CMI values produced over the Missouri River headwater subbasin strongly correlate (r = 0.75) with spring and early summer stream discharge, demonstrating the use of this metric to indicate moisture conditions for the snowmelt and growing seasons.

DEVELOP Project Summary↗
NASA NTRS2021

Montana Water Resources: Developing a Composite Moisture Index Utilizing NASA Earth Observations for Drought Monitoring in Montana and the Missouri River Basin

The Missouri River Basin provides irrigation water for a substantial part of the domestic agricultural sector within the United States. Drought events pose a significant threat to economic livelihoods of dependent individuals, industries, and ecosystems (e.g., farmers, local tribes, hydroelectric power, wildlife). In 2017 alone, the Missouri River Basin experienced a severe drought that resulted in a $2.6 billion loss to the U.S. Northern Plains. In response to such events, organizations throughout the basin, such as the Montana Climate Office, have dedicated efforts for drought monitoring and communicating relevant information to local stakeholders. In an effort to aid regional decision-making capabilities, the project team partnered with the Montana Climate Office, NOAA National Weather Service (NWS) Missouri Basin River Forecast Center, and NOAA Regional Climate Services, Central Region to create a monthly Composite Moisture Index (CMI) that relies on NASA Earth observations from the Terra Moderate Resolution Imaging Spectroradiometer (MODIS) and the Soil Moisture Active Passive (SMAP) mission. From these satellites, as well as the NOAA NWS National Operational Hydrologic Remote Sensing Center Snow Data Assimilation System (SNODAS), our team aggregated climate datasets including soil moisture, snow cover, snow depth, and snow water equivalent to compute a CMI that indicates regional moisture conditions during the winter months. The March CMI values produced over the Missouri River headwater subbasin strongly correlate (r = 0.75) with spring and early summer stream discharge, demonstrating the use of this metric to indicate moisture conditions for the snowmelt and growing seasons.

DEVELOP Tech Paper↗
DOE OSTI2023

The Dichotomy of Wet and Dry Trends Over India by Aerosol Indirect Effects in CMIP5 Models

Abstract Aerosol‐cloud interactions, also known as aerosol indirect effect (AIE), substantially impact rainfall frequency and intensity. Here, we analyze NEX‐GDDP, a multimodel ensemble of high‐resolution (0.25°) historical simulations and future projections statistically downscaled from 21 CMIP5 models, to quantify the importance of AIE on extreme climate indices, specifically consecutive dry days (CDD), consecutive wet days (CWD), and simple daily intensity index (SDII). The 21 NEX‐GDDP CMIP5 models are classified into models with reliable (REM) and unreliable (UREM) monsoon climate simulated over India based on their simulations of the climate indices. The REM group is further decomposed based on whether the models represent only the direct (REM ADE ) or the direct and indirect (REM ALL ) aerosol effects. Compared to REM ADE , including all aerosol effects significantly improves the model skills in simulating the observed historical trends of all three climate indices over India. Specifically, AIE enhances dry days and reduces wet days in India in the historical period, consistent with the observed changes. However, by the middle and end of the 21st century, there is a relative decrease in dry days and an increase in wet days and precipitation intensity. Moreover, the REM ALL simulated future CWD and CDD changes are mostly opposite to those in REM ADE , indicating the substantial role of AIE in the future projection of dry and wet climates. These findings underscore the crucial role of AIE in future projections of the Indian hydroclimate and motivate efforts to accurately represent AIE in climate models.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2000

A Footprint of Solar Variability Variability on Climate

An 11-year variation is shown here to exist in an index of the atmospheric variability defined as the difference between the North Atlantic Oscillation and Arctic Oscillation indices.

UV North Atlantic Oscillation El Nino↗
DOE OSTI2021

Heat Vulnerability Index Development and Mapping

Extreme heat is one of the leading causes of weather-related deaths in the U.S. Exposure to extreme heat will be exacerbated due to global climate change. It is thus crucial to design a key performance indicator, heat vulnerability index (HVI), to represent overall heat risk which can help identify susceptible regions and sub-populations in cities in the face of heatwaves. Most existing HVI tools only consider outdoor heat exposure. This paper developed an HVI web map tool incorporating both the outdoor and the indoor heat exposure, as well as population sensitivity and adaptation capability across census tracts in the city of Fresno, California. The tool can assist the planning of infrastructure and resources to reduce residents’ vulnerability to extreme heat events.

Xu, Yujie↗
NASA NTRS2000

Retrieving the Vertical Structure of the Effective Aerosol Complex Index of Refraction from a Combination of Aerosol in Situ and Remote Sensing Measurements During TARFOX

The largest uncertainty in estimates of the effects of atmospheric aerosols on climate stems from uncertainties in the determination of their microphysical properties, including the aerosol complex index of refraction, which in turn determines their optical properties. A novel technique is used to estimate the aerosol complex index of refraction in distinct vertical layers from a combination of aerosol in situ size distribution and remote sensing measurements during the Tropospheric Aerosol Radiative Forcing Observational Experiment (TARFOX). In particular, aerosol backscatter measurements using the NASA Langley LASE (Lidar Atmospheric Sensing Experiment) instrument and in situ aerosol size distribution data are utilized to derive vertical profiles of the "effective" aerosol complex index of refraction at 815 nm (i.e., the refractive index that would provide the same backscatter signal in a forward calculation on the basis of the measured in situ particle size distributions for homogeneous, spherical aerosols). A sensitivity study shows that this method yields small errors in the retrieved aerosol refractive indices, provided the errors in the lidar-derived aerosol backscatter are less than 30% and random in nature. Absolute errors in the estimated aerosol refractive indices are generally less than 0.04 for the real part and can be as much as 0.042 for the imaginary part in the case of a 30% error in the lidar-derived aerosol backscatter. The measurements of aerosol optical depth from the NASA Ames Airborne Tracking Sunphotometer (AATS-6) are successfully incorporated into the new technique and help constrain the retrieved aerosol refractive indices. An application of the technique to two TARFOX case studies yields the occurrence of vertical layers of distinct aerosol refractive indices. Values of the estimated complex aerosol refractive index range from 1.33 to 1.45 for the real part and 0.001 to 0.008 for the imaginary part. The methodology devised in this study provides, for the first time, a complete set of vertically resolved aerosol size distribution and refractive index data. yielding the vertical distribution of aerosol optical properties required for the determination of aerosol-induced radiative flux changes.

Redemann, J.↗
NASA NTRS2000

Retrieving the Vertical Structure of the Effective Aerosol Complex Index of Refraction from a Combination of Aerosol in Situ and Remote Sensing Measurements During TARFOX

The largest uncertainty in estimates of the effects of atmospheric aerosols on climate stems from uncertainties in the determination of their microphysical properties, including the aerosol complex index of refraction, which in turn determines their optical properties. A novel technique is used to estimate the aerosol complex index of refraction in distinct vertical layers from a combination of aerosol in situ size distribution and remote sensing measurements during the Tropospheric Aerosol Radiative Forcing Observational Experiment (TARFOX). In particular, aerosol backscatter measurements using the NASA Langley LASE (Lidar Atmospheric Sensing Experiment) instrument and in situ aerosol size distribution data are utilized to derive vertical profiles of the 'effective' aerosol complex index of refraction at 815 nm (i.e., the refractive index that would provide the same backscatter signal in a forward calculation on the basis of the measured in situ particle size distributions for homogeneous, spherical aerosols). A sensitivity study shows that this method yields small errors in the retrieved aerosol refractive indices, provided the errors in the lidar derived aerosol backscatter are less than 30% and random in nature. Absolute errors in the estimated aerosol refractive indices are generally less than 0.04 for the real part and can be as much as 0.042 for the imaginary part in the case of a 30% error in the lidar-derived aerosol backscatter. The measurements of aerosol optical depth from the NASA Ames Airborne Tracking Sunphotometer (AATS-6) are successfully incorporated into the new technique and help constrain the retrieved aerosol refractive indices. An application of the technique to two TARFOX case studies yields the occurrence of vertical layers of distinct aerosol refractive indices. Values of the estimated complex aerosol refractive index range from 1.33 to 1.45 for the real part and 0.001 to 0.008 for the imaginary part. The methodology devised in this study provides, for the first time a complete set of vertically resolved aerosol size distribution and refractive index data, yielding the vertical distribution of aerosol optical properties required for the determination of aersol-induced radiative flux changes

Redemann, J.↗
NASA NTRS2014

On the Relationship Between Global Land-Ocean Temperature and Various Descriptors of Solar-Geomagnetic Activity and Climate

Examined are sunspot cycle- (SC-) length averages of the annual January-December values of the Global Land-Ocean Temperature Index ( ) in relation to SC-length averages of annual values of various descriptors of solar-geomagnetic activity and climate, incorporating lags of 0-5 yr. For the overall interval SC12-SC23, the is inferred to correlate best against the parameter incorporating lag = 5 yr, where the parameter refers to the resultant aa value having removed that portion of the annual aa average value due to the yearly variation of sunspot number (SSN). The inferred correlation between the and is statistically important at confidence level cl > 99.9%, having a coefficient of linear correlation r = 0.865 and standard error of estimate se = 0.149 degC. Excluding the most recent cycles SC22 and SC23, the inferred correlation is stronger, having r = 0.969 and se = 0.048 degC. With respect to the overall trend in the , which has been upwards towards warmer temperatures since SC12 (1878-1888), solar-geomagnetic activity parameters are now trending downwards (since SC19). For SC20-SC23, in contrast, comparison of the against SC-length averages of the annual value of the Mauna Loa carbon dioxide ( ) index is found to be highly statistically important (cl >> 99.9%), having r = 0.9994 and se = 0.012 degC for lag = 2 yr. On the basis of the inferred preferential linear correlation between the and , the current ongoing SC24 is inferred to have warmer than was seen in SC23 (i.e., >0.526 degC), probably in excess of 0.68 degC (relative to the 1951-1980 base period).

Wilson, Robert M.↗
DOE OSTI2025

Causal relationships of vegetation productivity with root zone water availability and atmospheric dryness at the catchment scale

Abstract. This study explores the causal relationships between catchment water availability, vapor pressure deficit, and gross primary productivity (GPP) across 341 catchments in the contiguous US. Seasonal climatic, hydrological, and vegetation characteristics were represented using the Horton index, ecological aridity index, evaporative fraction index, and carbon uptake efficiency. Statistical methods, including circularity statistics, correlation analysis, and causality tests, were employed to determine the complex interactions between catchment wetness, atmospheric dryness, and vegetation carbon uptake. The results revealed a maximum lag of 2 months in the intra-annual variability of catchment water supply–productivity and atmospheric water demand–productivity relationships, with hysteresis patterns varying with the catchment's hydrological characteristics. In catchments not permanently under water-limited or energy-limited conditions, vegetation experiences hydrological stress during the peak growing period, coinciding with the highest gross primary productivity and carbon uptake efficiency being out of phase with the Horton index and in phase with the evaporative fraction index. Causality analysis highlights strong temporal continuity in GPP seasonal characteristics, with a cause–effect relationship between catchment water supply, atmospheric demand, and vegetation productivity spanning a maximum of 2 months. These findings underscore the need for a comprehensive functional framework that integrates catchment water supply, atmospheric demand, and vegetation productivity to enhance our understanding and predictive capabilities with regard to ecosystem responses to climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Explaining the trends and variability in the United States tornado records using climate teleconnections and shifts in observational practices

The annual frequency of tornadoes during 1950–2018 across the major tornado-impacted states were examined and modeled using anthropogenic and large-scale climate covariates in a hierarchical Bayesian inference framework. Anthropogenic factors include increases in population density and better detection systems since the mid-1990s. Large-scale climate variables include El Niño Southern Oscillation (ENSO), Southern Oscillation Index (SOI), North Atlantic Oscillation (NAO), Pacific Decadal Oscillation (PDO), Arctic Oscillation (AO), and Atlantic Multi-decadal Oscillation (AMO). The model provides a robust way of estimating the response coefficients by considering pooling of information across groups of states that belong to Tornado Alley, Dixie Alley, and Other States, thereby reducing their uncertainty. The influence of the anthropogenic factors and the large-scale climate variables are modeled in a nested framework to unravel secular trend from cyclical variability. Population density explains the long-term trend in Dixie Alley. The step-increase induced due to the installation of the Doppler Radar systems explains the long-term trend in Tornado Alley. NAO and the interplay between NAO and ENSO explained the interannual to multi-decadal variability in Tornado Alley. PDO and AMO are also contributing to this multi-time scale variability. SOI and AO explain the cyclical variability in Dixie Alley. This improved understanding of the variability and trends in tornadoes should be of immense value to public planners, businesses, and insurance-based risk management agencies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Tree Diversity, Site Index, and Carbon Storage Decrease With Aridity in Douglas-Fir Forests in Western Canada

Forests are important for biodiversity, timber production and carbon accumulation, but these ecosystem services may be impacted by climate change. Field data collected from individual forest types occurring across a climatic gradient can contribute to forecasting these consequences. We examined how changes in temperature, precipitation and aridity affect ecosystem services in 23 mature Douglas-fir (Pseudotsuga menziesii) forests in nine climatic regions across a 900 km gradient in British Columbia, Canada. Using Canadian National Forest Inventory methodology, we assessed richness and diversity of plant functional groups, site index, and above- and below-ground carbon stocks. As aridity increased, ecosystem-level tree species richness declined on average from four to one species, Douglas-fir site index declined from 30 to 15 m, and ecosystem carbon storage decreased from 565 to 222 Mg ha –1 . Tree species richness was positively and herb species richness negatively correlated with carbon storage. Carbon storage by ecosystem compartment was largest in aboveground live tree biomass, declining in the following order: mineral soils > coarse woody debris and dead standing trees > forest floor > small and fine woody debris > understory plants. Mineral soil carbon at depths of 0-15 cm, 15-35 cm, and 35-55 cm increased with increasing mean annual precipitation and decreasing aridity. Our results indicate that as aridity increases and precipitation decreases, tree species richness, site index and carbon storage in existing Douglas-fir forests declines. However, assisted or natural migration of Douglas-fir into more humid regions could be associated with more diverse, productive, carbon-rich forests. This study informs carbon stock vulnerability and provides empirical data essential for carbon stock forecasts.

54 ENVIRONMENTAL SCIENCES↗