Search NASA⌕ Search

SEARCH · Search NASA

Results for “climate simulation”

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

Evaluation of native Earth system model output with ESMValTool v2.6.0

Earth system models (ESMs) are state-of-the-art climate models that allow numerical simulations of the past, present-day, and future climate. To extend our understanding of the Earth system and improve climate change projections, the complexity of ESMs heavily increased over the last decades. As a consequence, the amount and volume of data provided by ESMs has increased considerably. Innovative tools for a comprehensive model evaluation and analysis are required to assess the performance of these increasingly complex ESMs against observations or reanalyses. One of these tools is the Earth System Model Evaluation Tool (ESMValTool), a community diagnostic and performance metrics tool for the evaluation of ESMs. Input data for ESMValTool needs to be formatted according to the CMOR (Climate Model Output Rewriter) standard, a process that is usually referred to as “CMORization”. While this is a quasi-standard for large model intercomparison projects like the Coupled Model Intercomparison Project (CMIP), this complicates the application of ESMValTool to non-CMOR-compliant climate model output. In this paper, we describe an extension of ESMValTool introduced in v2.6.0 that allows seamless reading and processing of “native” climate model output, i.e., operational output produced by running the climate model through the standard workflow of the corresponding modeling institute. This is achieved by an extension of ESMValTool's preprocessing pipeline that performs a CMOR-like reformatting of the native model output during runtime. Thus, the rich collection of diagnostics provided by ESMValTool is now fully available for these models. For models that use unstructured grids, a further preprocessing step required to apply many common diagnostics is regridding to a regular latitude–longitude grid. Extensions to ESMValTool's regridding functions described here allow for more flexible interpolation schemes that can be used on unstructured grids. Currently, ESMValTool supports nearest-neighbor, bilinear, and first-order conservative regridding from unstructured grids to regular grids. Example applications of this new native model support are the evaluation of new model setups against predecessor versions, assessing of the performance of different simulations against observations, CMORization of native model data for contributions to model intercomparison projects, and monitoring of running climate model simulations. For the latter, new general-purpose diagnostics have been added to ESMValTool that are able to plot a wide range of variable types. Currently, five climate models are supported: CESM2 (experimental; at the moment, only surface variables are available), EC-Earth3, EMAC, ICON, and IPSL-CM6. As the framework for the CMOR-like reformatting of native model output described here is implemented in a general way, support for other climate models can be easily added.

58 GEOSCIENCES↗
DOE OSTI2023

Future Projections of the El Niño—Southern Oscillation and Tropical Pacific Mean State in CMIP6

Abstract The El Niño—Southern Oscillation (ENSO) is an important mode of tropical Pacific atmosphere‐ocean variability that drives teleconnections with weather and climate globally. However, prior studies using state‐of‐the‐art climate models lack consensus regarding future ENSO projections and are often impacted by tropical Pacific sea‐surface temperature (SST) biases. We used 173 simulations from 29 climate models participating in the Coupled Model Intercomparison Project, version 6 (CMIP6) to analyze model biases and future ENSO projections. We analyzed two ENSO indices, namely the ENSO Longitude Index (ELI), which measures zonal shifts in tropical Pacific deep convection and accounts for changes in background SST, and the Niño 3.4 index, which measures SST anomalies in the central‐eastern equatorial Pacific. We found that the warm eastern tropical‐subtropical Pacific SST bias typical of previous generations of climate models persists into many of the CMIP6 models. Future projections of ENSO shift toward more El Niño‐like conditions based on ELI in 48% of simulations and 55% of models, in association with a future weakening of the zonal equatorial Pacific SST gradient. On the other hand, none of the models project a significant shift toward La Niña‐like conditions. The standard deviation of the Niño 3.4 index indicates a lack of consensus on whether an increase or decrease in ENSO variability is expected in the future. Finally, we found a possible relationship between historical SST and low‐level cloud cover biases in the ENSO region and future changes in ELI; however, this result may be impacted by limitations in data availability.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Simulation of African Easterly Waves in a Global Climate Model

Abstract African easterly waves (AEWs) exert significant influence on local and downstream high-impact weather including tropical cyclone (TC) genesis over the Atlantic. Accurate representation of AEWs in climate and weather prediction models therefore is necessary for skillful predictions. In this study, we examine simulated AEWs, including their evolution, vertical structure, and linkage to tropical cyclone genesis, in the NASA Goddard Earth Observing System Model, version 5 (GEOS-5) atmospheric global climate model. Identified by the leading empirical orthogonal function mode of time-filtered precipitation, the observed westward propagating AEWs along the southern track over the Atlantic are largely captured in GEOS-5, but with a slower phase speed and significantly weaker amplitude downstream off the West Africa coast. The weak downstream development of AEWs in GEOS-5 is accompanied by much reduced TC genesis over the main development region. Further analyses suggest that the slow westward propagation and weaker AEW amplitude downstream can be ascribed to a weak African easterly jet, while overestimated negative (positive) meridional potential vorticity (PV) gradients appear to the north (south) of 10°N in GEOS-5. The greatly overestimated positive meridional PV gradient to the south of 10°N is expected to generate strong horizontal stretching in the AEW wave pattern in the model, which hinders organization of convection and its feedback to sustain the AEW development. Persistent and vigorous AEW precipitation in the Guinea Highlands of the West Africa coast could also be responsible for reduced westward propagation of AEWs in the model.

Meteorology & Atmospheric Sciences↗
DOE OSTI2021

A Study of Cost-Saving Potential of Load Flexibility Measures in Grid-Interactive Multifamily Buildings

With recent advances in smart technologies, more and more smart devices are penetrating the residential and commercial buildings market. The introduction of these smart devices is also helping IoT companies emerge with load aggregator roles in the sector. With more utility companies on the track of supporting OpenADR protocols, the aggregators could play a significant role in providing load flexibilities by automatically responding to demand response (DR) events and coordinating load flexibility measures between customers. This would benefit utility companies by reducing stress on the grid during critical peak demand hours as well as customers by allowing them to utilize utility rate structures advantageous to those able to reduce electric usage during high-demand hours. This study evaluates cost and energy savings from adopting multiple load flexibility measures in multifamily buildings. Combinations of different load flexibility measures, including space temperature floating, light dimming, automatic window shading, and water heater temperature floating, are considered. The simulations are performed using OpenStudio®, an open-source U.S. Department of Energy (DOE) simulation platform. For the case study, we used a midrise apartment building with weather conditions from Denver, Colorado. To compare climate zone differences, simulations were also performed for Los Angeles, California, and Chicago, Illinois. Initial results indicate that the application of automated load flexibility measures without careful consideration of dispatching strategies and DR program enrollments could significantly affect the savings. To get meaningful cost savings, aggregators need to encourage tenant awareness to curtail energy usage through occupant behavior in addition to dispatching automatic load flexibility measures. The outcomes from this study are believed to help load aggregators understand the risks and benefits of load flexibility opportunities.

building energy modeling↗
DOE OSTI2025

Insights Into Nature-Based Climate Solutions: Managing Forests for Climate Resilience and Carbon Stability

Successful implementation of forest management as a nature-based climate solution is dependent on the durability of management-induced changes in forest carbon storage and sequestration. As forests face unprecedented stability risks in the face of ongoing climate change, much remains unknown regarding how management will impact forest stability, or how interactions with climate might shift the response of forests to management across spatiotemporal scales. Here, we used a process-based model to simulate multidecadal projections of forest dynamics in response to changes in management and climate. Simulations were conducted across gradients in forest type, edaphic factors, and management intensity under two alternate radiative forcing scenarios (RCP4.5 and RCP8.5). This allowed for the quantification of forest stability shifts in response to climate change, and the role of management in modulating that response, where ecosystem stability is characterized as the resilience and temporal stability of net primary production, aboveground biomass, and soil carbon. Our results indicate that forest structure is primarily shaped by management, but the same management strategy often produced divergent structures over time, due to interactions with regional climate change. We found that management can be used to increase stability and minimize the release of stored carbon by reducing mortality, but also highlight the regional dependency of management-induced changes in resilience to climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Natural variability contributes to model–satellite differences in tropical tropospheric warming

Significance Climate models have, on average, simulated substantially more tropical tropospheric warming than satellite data, with few simulations matching observations. It has been suggested that this discrepancy arises because climate models are overly sensitive to greenhouse gas increases. Tropical tropospheric temperature trends from a large ensemble of simulations performed with a single climate model span a wide range that is solely due to natural climate variability. A subset of these simulations have warming rates in accord with satellite observations. Simulations with diminished tropical tropospheric warming due to climate variability exhibit characteristic patterns of surface warming that are similar to the observed record. Our results indicate that multidecadal variability can explain current model–observational differences in the rate of tropical tropospheric warming.

58 GEOSCIENCES↗
DOE OSTI2022

A Multi-Model Framework for Assessing Long- and Short-Term Climate Influences on the Electric Grid

Climate change influences many aspects of the electric grid, but prior work and industry practices often ignore the potential effects of changing climate, or they only consider a single effect or individual effects in isolation. Challenges vary with each grid and include adapting to long-term trends such as changing temperature and precipitation or shorter-term events such as drought or storms that could increase in frequency or intensity. Here we present a multi-model framework designed to analyze the effects of long and short-term climate impacts in combination. This framework couples capacity expansion and production cost models with hydrologic models and future climate scenario data to analyze alternative climate and energy futures at high spatial, temporal, and process resolutions. Furthermore, we constructed and evaluated the results of a suite of simulated scenarios exploring climate impacts on capacity investment and stress-tested the resulting future infrastructures using hourly dispatch modeling under alternative drought and load conditions. We demonstrate the approach through a case study of the U.S. Western Interconnection, where climate impacts depend on interactions between temperature-induced load, water availability for hydropower, technology competitiveness, and demand flexibility. Changes in 2038 generating capacity range from -8.5-16.6 GW, and changes in 2038 transmission capacity range from -1-2 GW. Capacity increases are driven by higher load from higher temperatures, while capacity reductions can be achieved in scenarios with higher future hydropower availability and increased demand flexibility. Scenarios requiring additional capacity cost an additional $\$5$-$\$17$ billion (discounted) from 2018 to 2038; however, scenarios with capacity reductions cost $\$1$-$\$18$ billion less. Stress tests on four 2038 infrastructures demonstrated that the identified systems were able to serve at least 99.999% of load and 99.96% of reserves. However, drought and unexpected high-load conditions can result in reduced capacity to respond to contingency events we did not model. Although these results are system and scenario specific, they highlight the importance of considering multiple climate change impacts simultaneously in long-term planning efforts and demonstrate a multi-model, multiscale approach that can be flexibly applied to any system and set of climate change concerns.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Modes of Variability in E3SM and CESM Large Ensembles

An adequate characterization of internal modes of climate variability (MoV) is a prerequisite for both accurate seasonal predictions and climate change detection and attribution. Assessing the fidelity of climate models in simulating MoV is therefore essential; however, doing so is complicated by the large intrinsic variations in MoV and the limited span of the observational record. Large ensembles (LEs) provide a unique opportunity to assess model fidelity in simulating MoV and quantify intermodel contrasts. Here, in this work, these goals are pursued in four recently produced LEs: the Energy Exascale Earth System Model (E3SM) versions 1 and 2 LEs, and the Community Earth System Model (CESM) versions 1 and 2 LEs. In general, the representation of global coupled modes is found to improve across successive E3SM and CESM versions in conjunction with the fidelity of the base state climate while the patterns of extratropical modes are well simulated across the ensembles. Various persistent shortcomings for all MoV are however identified and discussed. The results both demonstrate the successes of these recent model versions and suggest the potential for continued improvement in the representation of MoV with advances in model physics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Simulating Observations of Southern Ocean Clouds and Implications for Climate

Southern Ocean (S. Ocean) clouds are important for climate prediction. Yet previous global climate models failed to accurately represent cloud phase distributions in this observation-sparse region. In this study, data from the Southern Ocean Clouds, Radiation, Aerosol, Transport Experimental Study (SOCRATES) experiment is compared to constrained simulations from a global climate model (the Community Atmosphere Model, CAM). Nudged versions of CAM are found to reproduce many of the features of detailed in situ observations, such as cloud location, cloud phase, and boundary layer structure. The simulation in CAM6 has improved its representation of S. Ocean clouds with adjustments to the ice nucleation and cloud microphysics schemes that permit more supercooled liquid. Comparisons between modeled and observed hydrometeor size distributions suggest that the modeled hydrometeor size distributions represent the dual peaked shape and form of observed distributions, which is remarkable given the scale difference between model and observations. Comparison to satellite observations of cloud physics is difficult due to model assumptions that do not match retrieval assumptions. Some biases in the model's representation of S. Ocean clouds and aerosols remain, but the detailed cloud physical parameterization provides a basis for process level improvement and direct comparisons to observations. This is crucial because cloud feedbacks and climate sensitivity are sensitive to the representation of S. Ocean clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Climate impacts in scenarios: time to close the loop?

Reaching a full understanding of the consequences of climate change for society and ecosystems, and the ensuing needs for adaptation, requires a consideration of the interactions between human and Earth systems. Currently, however, climate research largely separates the influence of society on climate from the influence of climate on society; that is, it doesn’t “close the loop.” A primary example of this approach is the generation and use of earth system model (ESM) simulations in the climate change research community. Large-scale socio-economic models, known as integrated assessment models (IAMs), are used to project emissions and land use change which serve as input to ESMs. ESM projections then serve as input to models of impacts on society and ecosystems. But, according to this modeling chain, those impacts do not affect the emissions and land use that drove the ESMs in the first place. Previous work has not drawn firm conclusions on whether this feedback would be large enough to warrant explicitly accounting for it. Two prominent possibilities, however, are that emissions and land use scenarios representing the high and low ends of the plausible range of future climate change are both too extreme. The high-end scenario may miss damaging impacts that would reduce economic activity, and therefore emissions, while the low-end scenario may ignore climate feedbacks that would make large-scale land-based carbon removal ineffective and therefore would hamper mitigation at the level assumed by the scenario. In this piece, we identify the opportunities and challenges that implementing such feedback loops would face. We argue that recent developments in climate impact research, human system modeling and ESM emulation make the time ripe to use IAMs in a structured model intercomparison exercise. Model intercomparison projects have benefitted the climate modeling community for decade snow, and more recently have also benefitted the impact modeling community. An IAM intercomparison focused on integrating impacts could make large strides in testing the implications of these feedbacks and assessing whether closing the loop would fundamentally change our outlook on future climate changes and their consequences.

Tebaldi, Claudia↗
DOE OSTI2023

Vertical Structure and Ice Production Processes of Shallow Convective Postfrontal Clouds over the Southern Ocean in MARCUS. Part II: Modeling Study

Abstract Part I of this series presented a detailed overview of postfrontal mixed-phase clouds observed during the Measurements of Aerosols, Radiation, and Clouds over the Southern Ocean (MARCUS) field campaign. In Part II, we focus on a multiday (23–26 February 2018) case with the aim of understanding ice production as well as model sensitivity to ice process parameterizations using the Weather Research and Forecasting (WRF) Model. The control simulation with the Predicted Particle Properties (P3) microphysics scheme underestimates the ice content and overestimates the supercooled liquid water content, contrary to the bias common in global climate models. The simulations targeted at ice production processes show negligible sensitivity to cloud droplet number concentrations. Further, neither increasing ice nuclei particle (INP) concentrations to an unrealistic level nor adjusting it to MARCUS field estimations alone guarantees more ice production in the model. However, the simulated clouds are found to be highly sensitive to the implementation of immersion freezing, the thresholding of condensation/deposition freezing initiation, and the rime splintering process. By increasing immersion freezing of cloud droplets, relaxing thresholds for condensation/deposition freezing, or removing rime splintering thresholds, the model significantly improves its performance in producing ice. The relaxation of the immersion freezing temperature threshold to the observed cloud-top temperature suggests an in-cloud seeder–feeder mechanism. The results of this work call for an increase in observations of INP, especially over the remote Southern Ocean and at relatively high temperatures, and measurements of ice particle size distributions to better constrain ice nucleating processes in models.

Meteorology & Atmospheric Sciences↗
DOE OSTI2025

An analysis of Potential Benefits of Adding Diurnal Thermal Storage in Geothermal Heat Pump Systems

Geothermal heat pumps (GHP) are highly efficient for space heating and cooling, utilizing the subsurface of the ground as seasonal energy storage through borehole heat exchangers (i.e., heat stored in summer is used to provide heating in winter, and cold stored in winter is used to provide cooling in summer). On the other hand, diurnal thermal energy storage can actively be charged and discharged for load shifting. Conventional thermal energy storage uses tanks filled with thermal storage media (such as phase change materials) and heat exchangers. It occupies building floor space or land areas outside the building, which may not be available especially in residential buildings. This paper will introduce a novel configuration that can utilize existing GHP system for both diurnal and seasonal thermal storage. This system can actively charge the ground on daily basis by injecting heat or cold in borehole heat exchangers. A simulation model was created for the novel configuration. Computer simulations were conducted using the model for the new GHP system serving a typical single-family house in three U.S. cities representing hot, warm, and cold climates. Based on simulation results, the potential benefits of load shifting and downsize of the expensive ground heat exchangers used by the GHP system were assessed.

Anees, Fady [ORNL]↗
DOE OSTI2022

Hygrothermal simulation of exterior retrofits in a cold climate

This work presents results from moisture modeling as part of a project undertaken by Pacific Northwest National Laboratory and its research partners, the University of Minnesota, and the Oak Ridge National Laboratory. The research goal is to identify exterior wall retrofit systems for cold climates that are low cost, energy efficient, and do not result in moisture durability problems. A base case wall was identified representing a typical wood frame construction, circa 1950. Retrofit options were selected based on input from industry, academia, and published work and seven options were constructed and installed at the University of Minnesota’s cold weather exposure facility in Cloquet, MN. Measurements were carried out during winter months and the data was used for model validation. Hygrothermal simulations were then carried out using WUFI (Version 6.4) in accordance with standard ANSI/ASRHAE 160-2016, Criteria for Moisture-Control Design Analysis in Buildings. Simulations were run out to three years and results show that the exterior retrofits improve thermal performance and do not negatively impact moisture durability of the existing wall.

Aldykiewicz Jr, Antonio↗
DOE OSTI2022

Forest regeneration within Earth system models: current process representations and ways forward

Earth system models must predict forest responses to global change in order to simulate future global climate, hydrology, and ecosystem dynamics. These models are increasingly adopting vegetation demographic approaches that explicitly represent tree growth, mortality, and recruitment, enabling advances in the projection of forest vulnerability and resilience, as well as evaluation with field data. To date, simulation of regeneration processes has received far less attention than simulation of processes that affect growth and mortality, in spite of their critical role maintaining forest structure, facilitating turnover in forest composition over space and time, enabling recovery from disturbance, and regulating climate-driven range shifts. Here, our critical review of regeneration process representations within current Earth system vegetation demographic models reveals the need to improve parameter values and algorithms for reproductive allocation, dispersal, seed survival and germination, environmental filtering in the seedling layer, and tree regeneration strategies adapted to wind, fire, and anthropogenic disturbance regimes. These improvements require synthesis of existing data, specific field data-collection protocols, and novel model algorithms compatible with global-scale simulations. Vegetation demographic models offer the opportunity to more fully integrate ecological understanding into Earth system prediction; regeneration processes need to be a critical part of the effort.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Seasonal Dependent Impact of Ice Cloud Longwave Scattering on the Polar Climate

Most climate models neglect cloud longwave (LW) scattering because scattering is considered negligible compared to strong LW absorption by clouds and greenhouse gases. While this rationale is valid for simulating extrapolar regions, it is questionable for the polar regions, where the atmosphere is dry and hence has weak absorption, and ice clouds that have strong scattering capability frequently occur. Using the slab-ocean Community Earth System Model, we show that ice cloud LW scattering can warm winter surface air temperature by 0.8–1.8 K in the Arctic and 1.3–1.9 K in the Antarctic, while this warming becomes much weaker in polar summer. Such scattering effect cannot be correctly assessed when sea surface temperature and sea ice are prescribed as this effect is manifested through a surface-atmosphere coupling. Furthermore, cloud LW scattering is a necessity for the correct simulation of polar climate and surface radiation budget, especially in the winter.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Observational constraints on low cloud feedback reduce uncertainty of climate sensitivity

Marine low clouds strongly cool the planet. How this cooling effect will respond to climate change is a leading source of uncertainty in climate sensitivity, the planetary warming resulting from CO 2 doubling. In this study, we observationally constrain this low cloud feedback at a near-global scale. Satellite observations are used to estimate the sensitivity of low clouds to interannual meteorological perturbations. Combined with model predictions of meteorological changes under greenhouse warming, this permits quantification of spatially resolved cloud feedbacks. We predict positive feedbacks from midlatitude low clouds and eastern ocean stratocumulus, nearly unchanged trade cumulus and a near-global marine low cloud feedback of 0.19 ± 0.12 W m -2 K -1 (90% confidence). These constraints imply a moderate climate sensitivity (~3 K). Despite improved midlatitude cloud feedback simulation by several current-generation climate models, their erroneously positive trade cumulus feedbacks produce unrealistically high climate sensitivities. Conversely, models simulating erroneously weak low cloud feedbacks produce unrealistically low climate sensitivities.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

SSP‐Based Land‐Use Change Scenarios: A Critical Uncertainty in Future Regional Climate Change Projections

Abstract To better understand the role projected land‐use changes (LUCs) may play in future regional climate projections, we assess the combined effects of greenhouse‐gas (GHG)‐forced climate change and LUCs in regional climate model (RCM) simulations. To do so, we produced RCM simulations that are complementary to the North‐American Coordinated Regional Downscaling Experiment (NA‐CORDEX) simulations, but with future LUCs that are consistent with particular Shared Socioeconomic Pathways (SSPs) and related to a specific Representative Concentration Pathway (RCP). We examine the state of the climate at the end of the 21st century with and without two urban and agricultural LUC scenarios that follow SSP3 and SSP5 using the Weather Research and Forecasting (WRF) model forced by one global climate model, the MPI‐ESM, under the RCP8.5 scenario. We find that LUCs following different societal trends under the SSPs can significantly affect climate projections in different ways. In regions of significant cropland expansion over previously forested area, projected annual mean temperature increases are diminished by around 0.5°C–1.0°C. Across all seasons, where urbanization is high, projected temperature increases are magnified. In particular, summer mean temperature projections are up to 4°C–5°C greater and minimum and maximum temperature projections are increased by 2.5°C–6°C, amounts that are on par with the warming due to GHG‐forced climate change. Warming is also enhanced in the urban surroundings. Future urbanization also has a large influence on precipitation projections during summer, increasing storm intensity, event length, and the overall amount over urbanized areas, and decreasing precipitation in surrounding areas.

Bukovsky, Melissa S.↗
DOE OSTI2023

Comparing the cooling effectiveness of operationalisable urban surface combination scenarios for summer heat mitigation

Extreme summer heat in cities exacerbates the vulnerability of urban communities to heatwaves. Vegetative and reflective urban surfaces can help reduce urban heat. This study investigated the impacts of urban trees, green roofs and cool roofs on heat mitigation during average and extreme summer conditions in temperate oceanic Melbourne, Australia. We simulated the city climate using ‘The Air Pollution Model’ (TAPM) at a 1 km spatial resolution over 10 years, which according to our review of the literature, was the most prolonged period for simulation in Melbourne. During a widespread heatwave event, some of the tested scenarios with combined surface parameters could reduce the extreme values of the energy budget components- sensible heat, latent heat, and storage heat fluxes up to seasonal averages compared to the existing situation for Melbourne (control). The scenario with the highest (reasonable maximum) ground-level vegetation, green roofs, and cool roofs could reduce air temperatures up to 2.4 °C. The simulations suggest that a combined strategy with vegetative and high-albedo surfaces will deliver higher effectiveness with maximum cooling benefits and cost-effectiveness than individual strategies in cities. These results suggest the importance of collaborative strategic planning of urban surfaces to make cities healthier, sustainable, and liveable.

54 ENVIRONMENTAL SCIENCES↗