Search NASA⌕ Search

SEARCH · Search NASA

Results for “Surface Radiation Budget”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Cloud radiative effect dominates variabilities of surface energy budget in the dark Arctic

Climate models simulate a wide range of temperatures in the Arctic. Here we investigate one of the main drivers of changes in surface temperature: the net surface heat flux in the models. We show that in the winter months of the dark Arctic, there is a more than two-fold difference in the net surface heat fluxes among the models, and this difference is dominated by the downward infrared radiation from clouds. Owing to the small amount of water vapor in the winter Arctic, infrared radiation from clouds transmits more easily to the surface in the Arctic than at other latitudes, resulting in large cloud radiative effect at the surface. The dominant role of the cloud effect is also found in the transient variability of the net surface heat flux. Results demonstrate that accurate simulation of clouds is crucial for determining the net surface heat flux, which in turn affects surface temperature and sea ice properties in the Arctic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Low-Level Liquid-Bearing Clouds Contribute to Seasonal Lower Atmosphere Stability and Surface Energy Forcing over a High-Mountain Watershed Environment

Measurements of atmospheric structure and surface energy budgets distributed along a high-altitude moun tain watershed environment near Crested Butte, Colorado, from two separate, but coordinated, field campaigns, Surface Atmosphere Integrated field Laboratory (SAIL) and Study of Precipitation, the Lower Atmosphere, and Surface for Hydrometeorology (SPLASH), are analyzed. This study identifies similarities and differences in how clouds influence the radiative budget over one snow-free summer season (2022) and two snow-covered seasons (2021/22; 2022/23) for this alpine location. A relationship between lower-tropospheric stability stratification and longwave radiative flux from the presence or absence of clouds is identified. When low clouds persisted, often with signatures of supercooled liquid in win ter, the lower troposphere experienced weaker stability, while radiatively clear skies that are less likely to be influenced by liquid droplets were associated with appreciably stronger lower-tropospheric stratification. Corresponding surface turbu lent heat fluxes partitioned differently based upon the cloud–stability stratification regime derived from early morning radiosounding profiles. Combined with the differences in the radiative budget largely resulting from dramatic seasonal dif ferences in surface albedo, the lower atmosphere stratification, surface energy budget, and near-surface thermodynamics are shown to be modified by the effective longwave radiative forcing of clouds. The diurnal evolution of thermodynamics and surface energy components varied depending on the early morning stratification state. Thus, the importance of quies cent versus synoptically active large-scale meteorology is hypothesized as a critical forcing for cloud properties and associ ated surface energy budget variations. The physical relationships between clouds, radiation, and stratification can provide a useful suite of metrics for process understanding and to evaluate numerical models in such an undersampled, highly com plex terrain environment.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Integrating Reanalysis and Satellite Cloud Information to Estimate Surface Downward Long-Wave Radiation

The estimation of downward long-wave radiation (DLR) at the surface is very important for the understanding of the Earth’s radiative budget with implications in surface–atmosphere exchanges, climate variability, and global warming. Theoretical radiative transfer and observationally based studies identify the crucial role of clouds in modulating the temporal and spatial variability of DLR. In this study, a new machine learning algorithm that uses multivariate adaptive regression splines (MARS) and the combination of near-surface meteorological data with satellite cloud information is proposed. The new algorithm is compared with the current operational formulation used by the European Organization for the Exploitation of Meteorological Satellites (EUMETSAT) Satellite Application Facility on Land Surface Analysis (LSA-SAF). Both algorithms use near-surface temperature and dewpoint temperature along with total column water vapor from the latest European Centre for Medium-range Weather Forecasts (ECMWF) reanalysis ERA5 and satellite cloud information from the Meteosat Second Generation. The algorithms are trained and validated using both ECMWF-ERA5 and DLR acquired from 23 ground stations as part of the Baseline Surface Radiation Network (BSRN) and the Atmospheric Radiation Measurement (ARM) user facility. Results show that the MARS algorithm generally improves DLR estimation in comparison with other model estimates, particularly when trained with observations. When considering all the validation data, root mean square errors (RMSEs) of 18.76, 23.55, and 22.08 W·m –2 are obtained for MARS, operational LSA-SAF, and ERA5, respectively. The added value of using the satellite cloud information is accessed by comparing with estimates driven by ERA5 total cloud cover, showing an increase of 17% of the RMSE. The consistency of MARS estimate is also tested against an independent dataset of 52 ground stations (from FLUXNET2015), further supporting the good performance of the proposed model.

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↗
DOE OSTI2023

Sensitivity of Arctic Surface Temperature to Including a Comprehensive Ocean Interior Reflectance to the Ocean Surface Albedo Within the Fully Coupled CESM2

Abstract Almost all current climate models simplify the ocean surface albedo (OSA) by assuming the reflected solar energy without the ocean interior contribution. In this study, an improved ocean surface albedo scheme is incorporated into the Community Earth System Model version 2 (CESM2) to assess the sensitivity of Arctic surface temperature to including ocean interior reflectance to the OSA. Fully coupled CESM2 simulations with and without ocean interior reflectance are subsequently performed, we focus on the analysis of Arctic surface temperature responses. Incorporating ocean interior reflectance increases absorbed solar radiation and warms the ocean, enhancing seasonal heat storage and release across the Arctic Ocean, and increasing sea ice reduction and positive climate feedbacks that elevates Arctic surface temperature. Seasonal variations in air‐surface temperature differences induce changes in turbulent heat flux patterns, concurrently modifying dynamic advection and moisture processes that affect boundary layer humidity and low clouds, especially in winter. Based on partitioning physical processes in the thermodynamic energy equation, surface air warming is induced primarily through positive heating anomalies of vertical advection, latent heat release, and longwave radiative forcing. Through an examination of the surface energy budget, skin temperature warming is driven predominantly by increased downward longwave radiation, positive surface albedo feedback in summer, and increased conductive heat transport from the ocean particularly in winter. Significant effects of ocean interior reflectance on the Arctic Ocean, including sea surface warming and sea ice reduction, justify the importance of ocean interior reflectance in climate models for better understanding of ongoing Arctic climate changes.

Meteorology & Atmospheric Sciences↗
DOE OSTI2022

SNICAR-ADv4: a physically based radiative transfer model to represent the spectral albedo of glacier ice

Abstract. Accurate modeling of cryospheric surface albedo is essential for our understanding of climate change as snow and ice surfaces regulate the global radiative budget and sea-level through their albedo and mass balance. Although significant progress has been made using physical principles to represent the dynamic albedo of snow, models of glacier ice albedo tend to be heavily parameterized and not explicitly connected with physical properties that govern albedo, such as the number and size of air bubbles, specific surface area (SSA), presence of abiotic and biotic light absorbing constituents (LACs), and characteristics of any overlying snow. Here, we introduce SNICAR-ADv4, an extension of the multi-layer two-stream delta-Eddington radiative transfer model with the adding–doubling solver that has been previously applied to represent snow and sea-ice spectral albedo. SNICAR-ADv4 treats spectrally resolved Fresnel reflectance and transmittance between overlying snow and higher-density glacier ice, scattering by air bubbles of varying sizes, and numerous types of LACs. SNICAR-ADv4 simulates a wide range of clean snow and ice broadband albedo (BBA), ranging from 0.88 for (30 µm) fine-grain snow to 0.03 for bare and bubble-free ice under direct light. Our results indicate that representing ice with a density of 650 kg m−3 as snow with no refractive Fresnel layer, as done previously, generally overestimates the BBA by an average of 0.058. However, because most naturally occurring ice surfaces are roughened “white ice”, we recommend modeling a thin snow layer over bare ice simulations. We find optimal agreement with measurements by representing cryospheric media with densities less than 650 kg m−3 as snow and larger-density media as bubbly ice with a Fresnel layer. SNICAR-ADv4 also simulates the non-linear albedo impacts from LACs with changing ice SSA, with peak impact per unit mass of LACs near SSAs of 0.1–0.01 m2 kg−1. For bare, bubble-free ice, LACs actually increase the albedo. SNICAR-ADv4 represents smooth transitions between snow, firn, and ice surfaces and accurately reproduces measured spectral albedos of a variety of glacier surfaces. This work paves the way for adapting SNICAR-ADv4 to be used in land ice model components of Earth system models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Radiative feedbacks on land surface change and associated tropical precipitation shifts

Changes in land surface albedo and land surface evaporation modulate the atmospheric energy budget by changing temperatures, water vapor, clouds, snow and ice cover, and the partitioning of surface energy fluxes. Here idealized perturbations to land surface properties are imposed in a global model to understand how such forcings drive shifts in zonal mean atmospheric energy transport and zonal mean tropical precipitation. For a uniform decrease in global land albedo, the albedo forcing and a positive water vapor feedback contribute roughly equally to increased energy absorption at the top of the atmosphere (TOA), while radiative changes due to the temperature and cloud cover response provide a negative feedback and energy loss at TOA. Decreasing land albedo causes a northwards shift in the zonal mean intertropical convergence zone (ITCZ). The combined effects on ITCZ location of all atmospheric feedbacks roughly cancel for the albedo forcing; the total ITCZ shift is comparable to that predicted for the albedo forcing alone. For an imposed increase in evaporative resistance that reduces land evaporation, low cloud cover decreases in the northern mid-latitudes and more energy is absorbed at TOA there; longwave loss due to warming provides a negative feedback on the TOA energy balance and ITCZ shift. Imposed changes in land albedo and evaporative resistance modulate fundamentally different aspects of the surface energy budget. However, the pattern of TOA radiation changes due to the water vapor and air temperature responses are highly correlated for these two forcings because both forcings lead to near-surface warming.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Polar primary aerosols across the ocean-sea ice-snow-atmosphere interface: From sources to impacts

Primary aerosols play a critical role in polar climate systems, influencing cloud formation, precipitation, radiative balance, and surface energy budgets. This paper provides a comprehensive synthesis of primary aerosol sources, transformation and removal processes, and broader atmospheric impacts in polar regions, emphasizing their links to ocean and sea ice biogeochemistry. These aerosols (including sea salt, primary organic aerosol, and primary biological aerosol particles) originate from marine and cryospheric environments and are emitted through physical processes, such as wave breaking, bubble bursting, and blowing snow. Emission sources include seawater, sea ice, snow, and freshwater from river discharge and glacial runoff. Once airborne, these particles can serve as a chemical reservoir, influencing atmospheric composition and reactivity, and as seeds for cloud droplet and ice crystal formation, influencing cloud microphysics and polar climate. Despite their importance, many of the processes governing primary aerosol emissions and transformations remain poorly constrained. The most pressing knowledge gaps pertain to emission processes, limited spatiotemporal observational coverage, instrumentation constraints, parameterization development, and the integration of interdisciplinary expertise. To improve our understanding of primary aerosol drivers and their response to climate, future research efforts should prioritize strategically coordinated and cross-disciplinary process studies, advancements in measurement technologies and coverage, and close collaboration between modelers and observational scientists to inform and refine model parameterizations. As polar regions continue to undergo profound changes marked by increased precipitation, reduced sea and land ice, freshening oceans, and shifting ecosystem dynamics, characterizing present-day primary aerosol populations is vital. Improved understanding will be essential for anticipating future changes in aerosol-radiation and aerosol-cloud interactions and their implications for polar and global climate systems.

Aerosol-cloud↗
DOE OSTI2024

The Effect of Physically Based Ice Radiative Processes on Greenland Ice Sheet Albedo and Surface Mass Balance in E3SM

Abstract A significant portion of surface melt on the Greenland Ice Sheet (GrIS) is due to dark ice regions in the ablation zone, where solar absorption is influenced by the physical properties of the ice, light absorbing constituents (LACs), and the overlying crustal surface or melt ponds. Earth system models (ESMs) typically prescribe the albedo of ice surfaces as a constant value in the visible and near‐infrared spectral regions. This work advances ESM ice radiative transfer modeling by (a) incorporating a physically based radiative transfer model (SNow, ICe and Aerosol Radiation model Adding‐Doubling Version 4; SNICAR‐ADv4) into the Energy Exascale Earth System Model (E3SM), (b) determining spatially and temporally varying bare ice physical properties over the GrIS ablation zone from satellite observations to inform SNICAR‐ADv4, and (c) assessing the impacts on simulated GrIS albedo and surface mass balance associated with modeling of more realistic bare ice albedo. GrIS‐wide bare ice albedo in E3SMv2 is overestimated by ∼4% in the visible and ∼7% in the near‐infrared wavelengths compared to the Moderate Resolution Imaging Spectroradiometer. Our bare ice physical property retrieval method found that LACs, ice crustal surfaces, and melt ponds reduce visible albedo by 30% in the bare ice region of the GrIS ablation zone. The realistic bare ice albedo reduces surface mass balance by ∼145 Gt, or 0.4 mm of sea‐level equivalent between 2000 and 2021 compared to the default E3SM. This work highlights the importance of simulating bare ice albedo accurately and realistically to improve our ability to quantify changes in the GrIS surface mass and radiative energy budgets.

58 GEOSCIENCES↗
DOE OSTI2026

The two radiative states of the Arctic atmosphere and their impacts on the surface energy budget of sea ice

The surface energy budget (SEB) is a central regulator of Arctic climate and sea ice evolution, yet its processes remain poorly constrained due to sparse observations and complex, coupled surface-atmosphere interactions. This study leverages year-long measurements from the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) to provide the most comprehensive assessment to date of the central Arctic SEB and its modulation by atmospheric variability. Ship- and ice-based observations from October 2019 to September 2020 were used to directly measure or tightly constrain each term of the SEB, leading to exceptional energetic closure with the seasonal snow and ice mass balance. The analysis reveals strong seasonal transitions in atmosphere-surface energy transfer that are modulated by the atmospheric state and constrained by the ability of the surface temperature to respond. Classification of the atmosphere into its two dominant radiative states—the semi-transparent (ST) and opaque (OP)—highlights the central role of synoptic-scale variability in clouds. The ST atmospheric state dominated the long winter ice growth season, with limited cloudiness supporting persistent surface radiative cooling and ice growth. The OP state, associated with liquid-containing or thick ice clouds, became dominant in spring, with the combination of increased solar heating and cloud surface longwave warming driving ice and snow melt. Eddy covariance versus bulk approaches for deriving surface turbulent heat fluxes provide vastly different perspectives on the role of turbulence in modulating the SEB. These results establish a high-quality benchmark dataset for Arctic SEB studies and demonstrate how the balance of atmospheric radiative states exerts a first-order control on the annual evolution of the sea ice. The findings have broad implications for advancing observing technologies, understanding Arctic amplification, improving climate models, and predicting future sea ice change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Trends in Downwelling Longwave Radiance Over the Southern Great Plains

Downwelling longwave radiation is an important part of the surface energy budget. Spectral trends in the downwelling longwave radiance (DLR) provide insight into the radiative drivers of climate change. In this research, we process and analyze a 23-year DLR record measured by the Atmospheric Emitted Radiance Interferometer (AERI) at the U.S. Department of Energy Atmospheric Radiation Program Southern Great Plains (SGP) site. Two AERIs were deployed at SGP with an overlapping observation period of about 10 years, which allows us to examine the consistency and accuracy of the measurements and to account for discrepancies between them due to errors associated with the instruments themselves. We then analyzed the all-sky radiance trends in DLR, which are associated with the surface warming trend at SGP during this same period and also the complex changes in meteorological conditions. For instance, the observed radiance in the CO2 absorption band follows closely the near-surface air temperature variations. The significant positive radiance trends in weak absorption channels, such as in the wings of the CO 2 band and in the weak absorption channels in the H 2 O vibration-rotational band, show earlier detectability of climate change. The magnitude of the radiance trend uncertainty in the DLR record mainly results from internal climate variability rather than from measurement error, which highlights the importance of continuing the DLR spectral measurements to unambiguously detect and attribute climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Distinct roles of land cover in regulating spatial variabilities of temperature responses to radiative effects of aerosols and clouds

Surface temperature responses to radiative perturbations due to aerosols and clouds are complicated by the land surface properties. To disentangle these complexities, this study, from a terrestrial surface energy budget perspective, isolates the underlying biophysical processes from the instantaneous radiative effects of aerosols and clouds on surface temperature using the National Center for Atmospheric Research Community Earth System Model version 1.2.1. It is found that in comparison with the global heterogeneous distributions of instantaneous radiative perturbations at the surface induced by aerosols and clouds, the spatial variations of the corresponding surface temperature responses to aerosol direct radiative effects (DRE) during the daytime and cloud radiative effects (CRE) during the nighttime are amplified. It is because of the consistent global distribution of the local surface climate sensitivity (a function of land cover properties such as surface roughness and Bowen ratio) with those of daytime DRE and nighttime CRE. By applying identical anthropogenic aerosol and precursor emissions over eight major past, present and projected future anthropogenic aerosol emitting regions (i.e. Brazil, China, East Africa, India, Indonesia, South Africa, United States and Western Europe), surface temperature responses to aerosol radiative cooling in the daytime and cloud radiative warming in the nighttime over these regions positively regulated by local surface climate sensitivities are prominent.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Subtropical Eastern North Pacific SST bias in earth system models

Karthik Balaguru (PNNL), Luke Van Roekel (LLNL), L. Ruby Leung (PNNL) and Milena Veneziani (LLNL) Abstract: This study systematically evaluates the warm SST bias in the Subtropical Eastern North Pacific (SENP), a problem plaguing most CMIP6 models, using the Energy Exascale Earth System Model version 1 (E3SM). Compared to the model at its standard resolution (1o atmosphere, 30-60 km ocean), the SST bias improves in several regions at high resolution (0.25o atmosphere, 18-6 km ocean) except in subtropical eastern boundary regions, such as the SENP. A mixed layer heat budget analysis reveals that errors in the surface radiative fluxes and horizontal heat advection are primarily responsible for the bias. Analysis of HighResMIP models indicates that the lack of improvement in the SENP SST bias with resolution is not native to E3SM alone. A larger CMIP6 inter-model correlation of SST bias with biases in surface winds than radiation highlights the important role of large-scale circulation biases in the SENP SST bias in coupled models.

Balaguru, Karthik↗
DOE OSTI2024

The Extraordinary March 2022 East Antarctica “Heat” Wave. Part II: Impacts on the Antarctic Ice Sheet

Abstract Between 15 and 19 March 2022, East Antarctica experienced an exceptional heat wave with widespread 30°–40°C temperature anomalies across the ice sheet. In Part I, we assessed the meteorological drivers that generated an intense atmospheric river (AR) that caused these record-shattering temperature anomalies. Here, we continue our large collaborative study by analyzing the widespread and diverse impacts driven by the AR landfall. These impacts included widespread rain and surface melt that was recorded along coastal areas, but this was outweighed by widespread high snowfall accumulations resulting in a largely positive surface mass balance contribution to the East Antarctic region. An analysis of the surface energy budget indicated that widespread downward longwave radiation anomalies caused by large cloud-liquid water contents along with some scattered solar radiation produced intense surface warming. Isotope measurements of the moisture were highly elevated, likely imprinting a strong signal for past climate reconstructions. The AR event attenuated cosmic ray measurements at Concordia, something previously never observed. Last, an extratropical cyclone west of the AR landfall likely triggered the final collapse of the critically unstable Conger Ice Shelf while further reducing an already record low sea ice extent. Significance Statement Using our diverse collective expertise, we explored the impacts from the March 2022 heat wave and atmospheric river across East Antarctica. One key takeaway is that the Antarctic cryosphere is highly sensitive to meteorological extremes originating from the midlatitudes and subtropics. Despite the large positive temperature anomalies driven from strong downward longwave radiation, this event led to huge amounts of snowfall across the Antarctic interior desert. The isotopes in this snow of warm airmass origin will likely be detectable in future ice cores and potentially distort past climate reconstructions. Even measurements of space activity were affected. Also, the swells generated from this storm helped to trigger the final collapse of an already critically unstable Conger Ice Shelf while further degrading sea ice coverage.

Wille, Jonathan D.↗
DOE OSTI2020

Ecosystem-Level Energy and Water Budgets Are Resilient to Canopy Mortality in Sparse Semiarid Biomes

Climate-driven woody vegetation mortality is a defining feature of semiarid biomes that drives fundamental changes in ecosystem structure. However, the observed impacts of woody mortality on ecosystem-scale energy and water budgets and the responses of surviving vegetation are highly variable among studies in water-limited environments. A previous girdling manipulation experiment in a piñon-juniper woodland suggested that although ecosystem-scale evapotranspiration was not altered by large-scale piñon mortality, soil water content decreased and the surviving juniper experienced greater water stress than juniper in an undisturbed woodland. Here we experimentally explored to what extent mortality-induced changes in energy balance components can explain these results. We compared energy fluxes measured above two adjacent piñon-juniper woodlands where piñon girdling was implemented at one site and the other subsequently experienced large-scale natural piñon mortality. We found that the mortality-induced decrease in canopy area was not sufficient to alter surface reflectance, roughness, and partitioning between energy budget components at both sites. A radiative transfer model estimated that because of the sparse premortality canopy, surface reflectance is more sensitive to a large increase in understory leaf area than further loss of crown area. Increased water stress in the remaining juniper following both mortality events can be explained by an increase in radiation on the ground that promoted higher soil temperature and evaporation. We found similar responses of ecosystem and tree-level functions to both girdling and natural mortality. This suggests that girdling is an appropriate approach to explore the impact of tree mortality on ecosystem structure, function, and energy balance.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Measurements of aerosol microphysical and chemical properties in the central Arctic atmosphere during MOSAiC

The Arctic environment is transforming rapidly due to climate change. Aerosols’ abundance and physicochemical characteristics play a crucial, yet uncertain, role in these changes due to their influence on the surface energy budget through direct interaction with solar radiation and indirectly via cloud formation. Importantly, Arctic aerosol properties are also changing in response to climate change. Despite their importance, year-round measurements of their characteristics are sparse in the Arctic and often confined to lower latitudes at Arctic land-based stations and/or short high-latitude summertime campaigns. Here, we present unique aerosol microphysics and chemical composition datasets collected during the year-long Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition, in the central Arctic. These datasets, which include aerosol particle number concentrations, size distributions, cloud condensation nuclei concentrations, fluorescent aerosol concentrations and properties, and aerosol bulk chemical composition (black carbon, sulfate, nitrate, ammonium, chloride, and organics) will serve to improve our understanding of high-Arctic aerosol processes, with relevance towards improved modelling of the future Arctic (and global) climate.

54 ENVIRONMENTAL SCIENCES↗