SEARCH · Engineering Papers

Results for “Cloud parameterizations”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

All sources NASA NTRS DOE OSTI

NASA NTRS1988

Shortwave radiative parameterization of large atmospheric aerosols - Dust and water clouds

A parameterization of the single scattering properties of an aerosol for a dust layer in the solar spectral region is presented. The parameterization scheme accounts for the microphysical changes of an aerosol layer and is applied to dust and water clouds. A scaling parameter is presented which describes, simply and effectively, the optical properties of an aerosol layer with a large effective size parameter.

Ackerman, Steven A.↗

DOE OSTI2022

Investigation on cloud vertical structures based on Ka-band cloud radar observations at Wuhan in Central China

Based on vertically pointing observations of Ka-band frequency-modulated continuous millimeter wave cloud radar (MMCR) in 2020, we investigate cloud vertical structures at Wuhan in central China. Case comparison indicates that the cloud top locations in the MMCR observation are largely consistent with those in satellite data. The precipitation has the occurrence frequency of 9.5% over the year, with the maximum value of 16.3% in July and the submaximal value of 14.7% in January, showing the significant influences of monsoons on precipitation at Wuhan. As the humid, cloudy and rainy climate, the probability of precipitation and non-precipitation clouds increases gradually from 53.9% in January to a large value of 71.5% in July. The non-precipitation clouds have the probability of 44.5% in the year, with two prominent peaks at the reflectivities of -35 to -30 dBZ around 12 km in summer and - 37– -32 dBZ around 10 km in spring. The cloud top can reach the heights of 14.49, 16.62, 13.38 and 11.73 km from spring to winter, respectively. The monthly mean cloud layer thickness is the largest (2.18 km) in May and September, and the smallest (1.36 km) in November. The seasonally averaged base and top heights and thickness of clouds have small diurnal variations. The high cloud is predominant in the warm months from May to September, and in the other cool months, the middle cloud has the most fraction. The quantitative investigation of clouds is helpful in improving the cloud parameterization in climate models.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2000

Comparison of Cirrus Cloud Models: A Project of the GEWEX Cloud System Study (GCSS) Working Group on Cirrus Cloud Systems

The GEWEX Cloud System Study (GCSS, GEWEX is the Global Energy and Water Cycle Experiment) is a community activity aiming to promote development of improved cloud parameterizations for application in the large-scale general circulation models (GCMs) used for climate research and for numerical weather prediction. The GCSS strategy is founded upon the use of cloud-system models (CSMs). These are "process" models with sufficient spatial and temporal resolution to represent individual cloud elements, but spanning a wide range of space and time scales to enable statistical analysis of simulated cloud systems. GCSS also employs single-column versions of the parametric cloud models (SCMs) used in GCMs. GCSS has working groups on boundary-layer clouds, cirrus clouds, extratropical layer cloud systems, precipitating deep convective cloud systems, and polar clouds.

Starr, David O'C.↗

NASA NTRS2000

Top-of-Atmosphere Albedo Estimation from Angular Distribution Models using Scene Identification from Satellite Cloud Property Retrievals

The next generation of Earth radiation budget satellite instruments will routinely merge estimates of global top-of-atmosphere radiative fluxes with cloud properties. This information will offer many new opportunities for validating radiative transfer models and cloud parameterizations in climate models. In this study, five months of POLarization and Directionality of the Earth's Reflectances (POLDER) 670 nm radiance measurements are considered in order to examine how satellite cloud property retrievals can be used to define empirical Angular Distribution Models (ADMs) for estimating top-of-atmosphere (TOA) albedo. ADMs are defined for 19 scene types defined by satellite retrievals of cloud fraction and cloud optical depth. Two approaches are used to define the ADM scene types: The first assumes there are no biases in the retrieved cloud properties and defines ADMs for fixed discrete intervals of cloud fraction and cloud optical depth (fixed-tau approach). The second approach involves the same cloud fraction intervals, but uses percentile intervals of cloud optical depth instead (percentile-tau approach). Albedos generated using these methods are compared with albedos inferred directly from the mean observed reflectance field. Albedos based on ADMs that assume cloud properties are unbiased (fixed-tau approach) show a strong systematic dependence on viewing geometry. This dependence becomes more pronounced with increasing solar zenith angle, reaching approximately equals 12% (relative) between near-nadir and oblique viewing zenith angles for solar zenith angles between 60 deg and 70 deg. The cause for this bias is shown to be due to biases in the cloud optical depth retrievals. In contrast, albedos based on ADMs built using percentile intervals of cloud optical depth (percentile-tau approach) show very little viewing zenith angle dependence and are in good agreement with albedos obtained by direct integration of the mean observed reflectance field (less than 1% relative error). When the ADMs are applied separately to populations consisting of only liquid water and ice clouds, significant biases in albedo with viewing geometry are observed (particularly at low sun elevations), highlighting the need to account for cloud phase both in cloud optical depth retrievals and in defining ADM scene types. ADM-derived monthly mean albedos determined for all 5 deg x 5 deg latitude/longitude regions over ocean are in good agreement (regional RMS relative errors less than 2%) with those obtained by direct integration when ADM albedos inferred from specific angular bins are averaged together. Albedos inferred from near-nadir and oblique viewing zenith angles are the least accurate, with regional RMS errors reaching approximately 5-10% (relative). Compared to an earlier study involving ERBE ADMs, regional mean albedos based on the 19 scene types considered here show a factor of 4 reduction in bias error and a factor of 3 reduction in RMS error.

Loeb, N. G.↗

NASA NTRS2012

Constraints on the Profiles of Total Water PDF in AGCMs from AIRS and a High-Resolution Model

Atmospheric general circulation model (AGCM) cloud parameterizations generally include an assumption about the subgrid-scale probability distribution function (PDF) of total water and its vertical profile. In the present study, the Atmospheric Infrared Sounder (AIRS) monthly-mean cloud amount and relative humidity fields are used to compute a proxy for the second moment of an AGCM total water PDF called the RH01 diagnostic, which is the AIRS mean relative humidity for cloud fractions of 0.1 or less. The dependence of the second moment on horizontal grid resolution is analyzed using results from a high-resolution global model simulation.The AIRS-derived RH01 diagnostic is generally larger near the surface than aloft, indicating a narrower PDF near the surface, and varies with the type of underlying surface. High-resolution model results show that the vertical structure of profiles of the AGCM PDF second moment is unchanged as the grid resolution changes from 200 to 100 to 50 km, and that the second-moment profiles shift toward higher values with decreasing grid spacing.Several Goddard Earth Observing System, version 5 (GEOS-5), AGCM simulations were performed with several choices for the profile of the PDF second moment. The resulting cloud and relative humidity fields were shown to be quite sensitive to the prescribed profile, and the use of a profile based on the AIRS-derived proxy results in improvements relative to observational estimates. The AIRS-guided total water PDF profiles, including their dependence on underlying surface type and on horizontal resolution, have been implemented in the version of the GEOS-5 AGCM used for publicly released simulations.

Molod, Andrea↗

NASA NTRS1999

EAULIQ: The Next Generation

This report summarizes the design of a new version of the stratiform cloud parameterization called Eauliq; the new version is called Eauliq NG. The key features of Eauliq NG are: (1) a prognostic fractional area covered by stratiform cloudiness, following the approach developed by M. Tiedtke for use in the ECMWF model; (2) separate prognostic thermodynamic variables for the clear and cloudy portions of each grid cell; (3) separate vertical velocities for the clear and cloudy portions of each grid cell, allowing the model to represent some aspects of observed mesoscale circulations; (4) cumulus entrainment from both the clear and cloudy portions of a grid cell, and cumulus detrainment into the cloudy portion only; and (5) the effects of the cumulus-induced subsidence in the cloudy portion of a grid cell on the cloud water and ice there. In this paper we present the mathematical framework of Eauliq NG; a discussion of cumulus effects; a new parameterization of lateral mass exchanges between clear and cloudy regions; and a theory to determine the mesoscale mass circulation, based on the hypothesis that the stratiform clouds remain neutrally buoyant through time and that the mesoscale circulations are the mechanism which makes this possible. An appendix also discusses some time-differencing methods.

Randall, David A.↗

DOE OSTI2025

Stable Machine‐Learning Parameterization of Subgrid Processes in a Comprehensive Atmospheric Model Learned From Embedded Convection‐Permitting Simulations

Modern climate projections often suffer from inadequate spatial and temporal resolution due to computational limitations, resulting in inaccurate representations of sub-grid processes. A promising technique to address this is the multiscale modeling framework (MMF), which embeds a kilometer-resolution cloud-resolving model (CRM) within each atmospheric column of a host climate model to replace traditional convection and cloud parameterizations. Machine learning offers a unique opportunity to make MMF more accessible by emulating the embedded CRM and reducing its substantial computational cost. Although many studies have demonstrated proof-of-concept success of achieving stable hybrid simulations, it remains a challenge to achieve near operational-level success with real geography and comprehensive variable emulation that includes, for example, explicit cloud condensate coupling. In this study, we present a stable hybrid model capable of integrating for at least 5 years with near operational-level complexity, including coarse-grid geography, seasonality, explicit cloud condensate and wind predictions, and land coupling. Our model demonstrates skillful online performance, achieving a 5-year zonal mean tropospheric temperature bias within 2 K, water vapor bias within 1 g/kg, and a precipitation root mean square error of 0.96 mm/day. Key factors contributing to our online performance include an expressive U-Net architecture and physical thermodynamic constraints for microphysics. With microphysical constraints mitigating unrealistic cloud formation, our work is the first to demonstrate realistic multi-year cloud condensate climatology under the MMF framework. Despite these advances, online diagnostics reveal persistent biases in certain regions, highlighting the need for innovative strategies to further optimize online performance.

Hu, Zeyuan [NVIDIA Corporation, Santa Clara, CA (U↗

NASA NTRS2000

Comparison of Cirrus Cloud Models: A Project of the GEWEX Cloud System Study (GCSS) Working Group on Cirrus Cloud Systems

The GEWEX Cloud System Study (GCSS, GEWEX is the Global Energy and Water Cycle Experiment) is a community activity aiming to promote development of improved cloud parameterizations for application in the large-scale general circulation models (GCMs) used for climate research and for numerical weather prediction (Browning et al, 1994). The GCSS strategy is founded upon the use of cloud-system models (CSMs). These are "process" models with sufficient spatial and temporal resolution to represent individual cloud elements, but spanning a wide range of space and time scales to enable statistical analysis of simulated cloud systems. GCSS also employs single-column versions of the parametric cloud models (SCMs) used in GCMs. GCSS has working groups on boundary-layer clouds, cirrus clouds, extratropical layer cloud systems, precipitating deep convective cloud systems, and polar clouds.

Starr, David OC.↗

DOE OSTI2023

An Ensemble of Neural Networks for Moist Physics Processes, Its Generalizability and Stable Integration

Abstract With the recent advances in data science, machine learning has been increasingly applied to convection and cloud parameterizations in global climate models (GCMs). This study extends the work of Han et al. (2020, https://doi.org/10.1029/2020MS002076 ) and uses an ensemble of 32‐layer deep convolutional residual neural networks, referred to as ResCu‐en, to emulate convection and cloud processes simulated by a superparameterized GCM, SPCAM. ResCu‐en predicts GCM grid‐scale temperature and moisture tendencies, and cloud liquid and ice water contents from moist physics processes. The surface rainfall is derived from the column‐integrated moisture tendency. The prediction uncertainty inherent in deep learning algorithms in emulating the moist physics is reduced by ensemble averaging. Results in 1‐year independent offline validation show that ResCu‐en has high prediction accuracy for all output variables, both in the current climate and in a warmer climate with +4K sea surface temperature. The analysis of different neural net configurations shows that the success to generalize in a warmer climate is attributed to convective memory and the 1‐dimensional convolution layers incorporated into ResCu‐en. We further implement a member of ResCu‐en into CAM5 with real world geography and run the neural‐network‐enabled CAM5 (NCAM) for 5 years without encountering any numerical integration instability. The simulation generally captures the global distribution of the mean precipitation, with a better simulation of precipitation intensity and diurnal cycle. However, there are large biases in temperature and moisture in high latitudes. These results highlight the importance of convective memory and demonstrate the potential for machine learning to enhance climate modeling.

Meteorology & Atmospheric Sciences↗

NASA NTRS1993

Nonhydrostatic effects in numerical modeling of mesoscale convective systems and baroclinic waves

The present investigation is concerned with the role of convection upon mesoscale modeling results, particularly when the grid resolution becomes small enough that there is not a clear scale separation between the explicitly resolved circulations and the parameterized clouds. In those situations, the vertical accelerations in explicitly resolved circulations become strong enough that the hydrostatic assumption may no longer be valid. These concerns arise from interests in improving mesoscale modeling per se and in improving the subgrid-scale parameterizations in global models. The hydrostatic and the nonhydrostatic options of the Colorado State University Regional Atmospheric Modeling System were used to simulate dry gravity currents in two dimensions, using several different horizontal grid sizes. With horizontal grid intervals of 10 km or less, nonhydrostatic simulations produce wider and colder heads and weaker but wider forced updrafts than do the hydrostatic simulations. Comparing the hydrostatic and nonhydrostatic models show that the difference between the vertical mass fluxes is much less than the difference between the vertical velocities. When the grid is fine enough to resolve the head of the gravity current, horizontal convergence at the gust front extends upwards almost to the head of the cold air. Vertical mass flux in the forced updraft at the front varies with horizontal grid size mainly as a function of the height of the simulated head. For coarser grids, which do not resolve the head, vertical mass flux at all heights decreases with increasing horizontal grid size. A comparison on nonhydrostatic simulations with horizontal grid intervals of 1 km and 2 km illustrates how decreasing the grid size does not necessarily increase the intensity of the resolved circulation. The smaller grid enables the simulated gravity current to entrain a bubble of warm air behind the head, which results in a weaker circulation with a shorter head and weaker updraft.

Cohen, Charles↗

NASA NTRS2015

How Difficult is it to Reduce Low-Level Cloud Biases With the Higher-Order Turbulence Closure Approach in Climate Models?

Low-level clouds cover nearly half of the Earth and play a critical role in regulating the energy and hydrological cycle. Despite the fact that a great effort has been put to advance the modeling and observational capability in recent years, low-level clouds remains one of the largest uncertainties in the projection of future climate change. Low-level cloud feedbacks dominate the uncertainty in the total cloud feedback in climate sensitivity and projection studies. These clouds are notoriously difficult to simulate in climate models due to its complicated interactions with aerosols, cloud microphysics, boundary-layer turbulence and cloud dynamics. The biases in both low cloud coverage/water content and cloud radiative effects (CREs) remain large. A simultaneous reduction in both cloud and CRE biases remains elusive. This presentation first reviews the effort of implementing the higher-order turbulence closure (HOC) approach to representing subgrid-scale turbulence and low-level cloud processes in climate models. There are two HOCs that have been implemented in climate models. They differ in how many three-order moments are used. The CLUBB are implemented in both CAM5 and GDFL models, which are compared with IPHOC that is implemented in CAM5 by our group. IPHOC uses three third-order moments while CLUBB only uses one third-order moment while both use a joint double-Gaussian distribution to represent the subgrid-scale variability. Despite that HOC is more physically consistent and produces more realistic low-cloud geographic distributions and transitions between cumulus and stratocumulus regimes, GCMs with traditional cloud parameterizations outperform in CREs because tuning of this type of models is more extensively performed than those with HOCs. We perform several tuning experiments with CAM5 implemented with IPHOC in an attempt to produce the nearly balanced global radiative budgets without deteriorating the low-cloud simulation. One of the issues in CAM5-IPHOC is that cloud water content is much higher than in CAM5, which is combined with higher low-cloud coverage to produce larger shortwave CREs in some low-cloud prevailing regions. Thus, the cloud-radiative feedbacks are exaggerated there. The turning exercise is focused on microphysical parameters, which are also commonly used for tuning in climate models. The results will be discussed in this presentation.

Xu, Kuan-Man↗

DOE OSTI2024

Locally narrow droplet size distributions are ubiquitous in stratocumulus clouds

Marine stratocumulus clouds are the “global reflectors,” sharply contrasting with the underlying dark ocean surface and exerting a net cooling on Earth’s climate. The magnitude of this cooling remains uncertain in part owing to the averaged representation of microphysical processes, such as the droplet-to-drizzle transition in global climate models (GCMs). Current GCMs parameterize cloud droplet size distributions as broad, cloud-averaged gammas. Here, using digital holographic measurements of discrete stratocumulus cloud volumes, we found cloud droplet size distributions to be narrower at the centimeter scale, never resembling the cloud average. These local distributions tended to form pockets of similar-looking cloud regions, each characterized by a size distribution shape that is diluted to varying degrees. These observations open the way for new modeling representations of microphysical processes.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2007

Spectral Retrieval of Latent Heating Profiles from TRMM PR data: Moistening Estimates over Tropical Ocean Regions - Part 3

The global hydrological cycle is central to the Earth's climate system, with rainfall and the physics of precipitation formation acting as the key links in the cycle. Two-thirds of global rainfall occurs in the tropics with the associated latent heating (LH) accounting for threefourths of the total heat energy available to the Earth's atmosphere. In the last decade, it has been established that standard products of LH from satellite measurements, particularly TRMM measurements, would be a valuable resource for scientific research and applications. Such products would enable new insights and investigations concerning the complexities of convection system life cycles, the diabatic heating controls and feedbacks related to rne-sosynoptic circulations and their forecasting, the relationship of tropical patterns of LH to the global circulation and climate, and strategies for improving cloud parameterizations In environmental prediction models. However, the LH and water vapor profile or budget (called the apparent moisture sink, or Q2) is closely related. This paper presented the development of an algorithm for retrieving Q2 using 'TRMM precipitation radar. Since there is no direct measurement of LH and Q2, the validation of algorithm usually applies a method called consistency check. Consistency checking involving Cloud Resolving Model (CRM)-generated LH and 42 profiles and algorithm-reconstructed is a useful step in evaluating the performance of a given algorithm. In this process, the CRM simulation of a time-dependent precipitation process (multiple-day time series) is used to obtain the required input parameters for a given algorithm. The algorithm is then used to "~econsti-LKt~h"e heating and moisture profiles that the CRM simulation originally produced, and finally both sets of conformal estimates (model and algorithm) are compared each other. The results indicate that discrepancies between the reconstructed and CM-simulated profiles for Q2, especially at low levels, are larger than those for latent heat. Larger discrepancies in Q2 at low levels are due to moistening for non-precipitating region that algorithm cannot reconstruct. Nevertheless, the algorithm-reconstructed total Q2 profiles are in good agreement with the CRM-simulated ones.

Shige, S.↗

NASA NTRS2001

Parameterizations of Depositional Growth of Cloud Ice in a Bulk Microphysical Scheme

The conversion of cloud ice to snow by depositional growth, designated P(sub SFI), in the Goddard Cumulus Ensemble Model cloud physics parameterization is examined. The original formulation of P(sub SFI) is shown to produce excessive conversion of cloud ice to snow because of an implicit assumption that the relative humidity is 100% with respect to water even though the air may actually be quite less humid. Two possible corrections to this problem are proposed, the first involving application of a relative humidity dependent correction factor to the original formulation of P(sub SFI), and the second involving a new formulation of P(sub SFI) based on the equation for depositional growth of cloud ice.

Braun, Scott A.↗

NASA NTRS2015

Exploring the Inner Edge of the Habitable Zone with Fully Coupled Oceans

The role of rotation in planetary atmospheres plays an important role in regulating atmospheric and oceanic heat flow, cloud formation and precipitation. Using the Goddard Institute for Space Studies (GISS) three dimension General Circulation Model (3D-GCM) we demonstrate how varying rotation rate and increasing the incident solar flux on a planet are related to each other and may allow the inner edge of the habitable zone to be much closer than many previous habitable zone studies have indicated. This is shown in particular for fully coupled ocean runs -- some of the first that have been utilized in this context. Results with a 100m mixed layer depth and our fully coupled ocean runs are compared with those of Yang et al. 2014, which demonstrates consistency across models. However, there are clear differences for rotations rates of 1-16x present earth day lengths between the mixed layer and fully couple ocean models, which points to the necessity of using fully coupled oceans whenever possible. The latter was recently demonstrated quite clearly by Hu & Yang 2014 in their aquaworld study with a fully coupled ocean when compared with similar mixed layer ocean studies and by Cullum et al. 2014. Atmospheric constituent amounts were also varied alongside adjustments to cloud parameterizations (results not shown here). While the latter have an effect on what a planet's global mean temperature is once the oceans reach equilibrium they do not qualitatively change the overall relationship between the globally averaged surface temperature and incident solar flux for rotation rates ranging from 1 to 256 times the present Earth day length. At the same time this study demonstrates that given the lack of knowledge about the atmospheric constituents and clouds on exoplanets there is still a large uncertainty as to where a planet will sit in a given star's habitable zone.

solar flux↗

NASA NTRS1990

A scheme for parameterizing ice-cloud water content in general circulation models

The optical properties of ice clouds are a primary issue for climate and climate change. Evaluating these optical properties in three-dimensional models for studying climate will require a method to calculate the ice water content of such clouds. A procedure is developed to parameterize ice water content as a function of large-scale meteorological characteristics for use in circulation models in which the ice water content is not calculated by means of a three-dimensional prognostic equation for condensed water. The technique identifies large-scale flows in which ice clouds exist and calculates their ice water content by reconstructing the trajectory associated with cloud formation. As the cloud forms, its ice content changes both by deposition of ice from water vapor and by ice removal by sedimentation. The sedimentation process is found to modify significantly the ice water content expected from deposition alone. Ice water contents predicted by the parameterization are compared with aircraft observations collected in the middle latitudes and the tropics, and show reasonable agreement over four orders-of-magnitude of ice water content. A parameterization for the sublimation of ice crystals settling into ice-subsaturated environments is also presented.

Heymsfield, Andrew J.↗

DOE OSTI2022

Annual cycle observations of aerosols capable of ice formation in central Arctic clouds

The Arctic is warming faster than anywhere else on Earth, prompting glacial melt, permafrost thaw, and sea ice decline. These severe consequences induce feedbacks that contribute to amplified warming, affecting weather and climate globally. Aerosols and clouds play a critical role in regulating radiation reaching the Arctic surface. However, the magnitude of their effects is not adequately quantified, especially in the central Arctic where they impact the energy balance over the sea ice. Specifically, aerosols called ice nucleating particles (INPs) remain understudied yet are necessary for cloud ice production and subsequent changes in cloud lifetime, radiative effects, and precipitation. Here, we report observations of INPs in the central Arctic over a full year, spanning the entire sea ice growth and decline cycle. Further, these observations are size-resolved, affording valuable information on INP sources. Our results reveal a strong seasonality of INPs, with lower concentrations in the winter and spring controlled by transport from lower latitudes, to enhanced concentrations of INPs during the summer melt, likely from marine biological production in local open waters. This comprehensive characterization of INPs will ultimately help inform cloud parameterizations in models of all scales.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Marine Shallow Cloud Adjustments to the Presence of Shortwave-Absorbing Aerosols: Advancing Understanding Through a Combined Analysis of Lasic Datasets and Process Modeling (Final Report)

This is a final report for this award, whose purpose was to support the analysis of data emanating from the LASIC (Layered Atlantic Smoke Interaction with Clouds) campaign. (Layered Atlantic Smoke Interactions with Clouds) is a strategy to improve our understanding of aged carbonaceous aerosol, its seasonal evolution, and the mechanisms by which clouds adjust to the presence of the aerosol. The observational strategy centers on deploying the AMF1 cloud, aerosol, and atmospheric profiling instrumentation to Ascension Island, located within the trade-wind shallow cumulus regime (14.50W, 80S) 3000 km offshore of continental Africa. The location is within the latitude zone of the maximum outflow of aerosol, with the deepening boundary layer known to entrain free-tropospheric smoke. The primary activities for LASIC are: 1) to improve current knowledge on aged biomass burning aerosol and its radiative properties as a function of the seasonal cycle; 2) to use surface-based remote sensing to sensitively interrogate the atmosphere for the relative vertical location of aerosol and clouds; 3) to improve our understanding of the cloud adjustments to the presence of shortwave-absorbing aerosol within the vertical column, both through aerosol-radiation and through aerosol-cloud interactions; 4) to aid low cloud parameterization efforts for climate models. This award serves to further these activities and includes a sub-award to Dr. Pablo Saide at UCLA for activities 3) and 4).

54 ENVIRONMENTAL SCIENCES↗