Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Urban”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Differences in urban plant community compositions across an urban-rural gradient in Knoxville, TN

Urban forests, or vegetation in areas under heavy human influence, provide many ecosystem services to urban residents such as localized cooling via evapotranspiration, shade, filtering of air pollution, and the associated health benefits of natural spaces. In order to quantify the magnitude of localized cooling by trees growing in varying levels of urbanization (based on % impervious surfaces, e.g., buildings, pavement), urban forest species composition, tree size, and tree density must be characterized. As a part of Oak Ridge National Laboratory’s (ORNL) urban forest temperature study, we conducted tree censuses in five Knoxville city parks where ORNL meteorological stations are deployed. Moreover, we measured every woody plant ≥ 5 cm diameter at breast height (DBH) within a 50 m radius of each site’s meteorological station for its DBH and species identification. When possible, individuals were identified down to species. Certain genera (Quercus spp., Carya spp., Pinus spp.) were identified down to genera in interest of time. Individual and total site basal area were calculated from measured DBH data. Results show notable differences in urban plant community compositions and total woody plant basal area across sites, with more urban sites closer to downtown (West View and SEEED) having lower tree basal area than the more suburban sites (West Hills, Cumberland Estates, and Victor Ashe). We identified 54 species across all sites, with West Hills and Victor Ashe having the highest species diversity. Our results show differences in forest compositions and sizes across Knoxville, which are currently informing ORNL’s evapotranspiration estimates for each site. Data Summary: Census data for West Hills (WH), Cumberland Estates (CE), Victor Ashe (VA), West View (WV), and Socially Equal Energy Efficient Development or SEEED (SD) urban forests in Knoxville, TN, USA, including tree size based on diameter at breast height (DBH; 1.3 m), species identification (Latin and common names), and basal area per stem (BA=π×[.5*DBH]^2). Field data are summarized in this file: “Community_Composition_Data.CSV”. Site-specific data detailing each site’s coordinates, number of stems measured at DBH, average tree DBH, α-diversity (number of species present), and total site basal area (sum of individual basal areas per site) are in this file: “Site_Comparisons.CSV”.

Warren, Jeffrey [ORNL] (ORCID:0000000206804697)↗
DOE OSTI2026

Urban Parameters Arizona Urban Corridor 100m

132 Urban parameters based on building physical dimensions and location were generated for the cities in six Arizona Counties at 100m resolution using the NATURF model. To use the binary file with WRF, the binary file and the index file must be placed in their own directory in WRF_GEOG and accessed in the same way NUDAPT44 would be accessed.

Dumas, Melissa [ORNL] (ORCID:0000000233190846)↗
DOE OSTI2024

Urban morphology and urban water demand: a case study in the land constrained Los Angeles region using urban growth modeling

The interactions between population growth, urban morphology, and water demand have important implications for water resources and supply in urban regions. Water use for irrigation comprises a significant fraction of urban water demand, and is potentially influenced by long-term changes in urban morphology. To investigate this, we used spatially explicit projections of urban land development intensity (fraction impervious area) generated from a 30 m resolution urban growth model for the Los Angeles (LA) region. Recent historical data on water use and high resolution landcover were used to establish relationships between green area, urban development intensity, and outdoor water demand. These relationships were then used to project outdoor and total water demand in 2100 using the urban growth model outputs. We considered two different population scenarios informed by the shared socioeconomic pathway (SSP) projections for the region (SSP3 and SSP5), and three scenarios of urban development intensification. Our analysis is resolved for over 80 water providers in the region, from the urban core to suburban fringe, and highlights diverse demand responses influenced by initial urban form and water demand attributes. Assumptions about outdoor water use factors based on recent water supply data were found to be nearly as influential on future outdoor demand as the urban growth scenario settings. Compared to previous studies, our work is unique in coherently linking high resolution SSP population scenarios, urban land cover evolution, and urban water demand projections, demonstrating the approach for the LA region—the largest population center in the western United States.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Projected Urban Morphology of the Los Angeles Area by the Year 2100

This dataset provides projections of urban building morphologies for the Los Angeles urban area at 30-meter spatial resolution. It contains 192 raster files that detail two primary building attributes: building footprint fractions (ranging from 0 to 1) and average building heights (ranging from 0 to 75 meters). The projections account for a wide range of future pathways, covering two Shared Socioeconomic Pathway (SSP) scenarios (SSP3 and SSP5), two population scenarios, two developed land intensification scenarios, and four distinct levels of intensification. The dataset was created using dual Generative Adversarial Networks (GANs) trained on 2015 land cover and building properties from the National Land Cover Database (NLCD) and Model America datasets. Supporting information on the dataset has been described in the LAUrbanAreaMorphologyProjections2100_README.txt file.

Pandey, Bhartendu↗
DOE OSTI2026

Will cities keep getting hotter? The interplay of urban expansion and greening reshapes future urban heat trajectories

Urban heat islands (UHIs) pose growing risks to public health, infrastructure, and resilience. While often assumed to intensify with urban growth, dynamic changes in urban expansion and vegetation greenness complicate UHI trajectories, which remain poorly understood. This study investigated the interplay of urban expansion and greenness change on UHI spatial profile across 36 Chinese megacities during 2003–2018 using multiple satellite products. We introduce a framework that classifies urban areas into four dynamic development pathways based on impervious surface area (ISA) and enhanced vegetation index (EVI) trends: urbanized-greening, urbanized-browning, urbanizing-greening (UingG), and urbanizing-browning (UingB). While most urbanized centers exhibited greening driven by targeted initiatives and urbanizing suburbs showed browning due to vegetation loss, about 30% urban areas showed the reverse pattern, revealing overlooked complexity in urban development. Urban expansion and browning strengthened UHI in suburban areas, whereas greening initiatives mitigated UHI in urban center and mitigated UHI enhancement in suburban areas. Slowed warming in urban centers together with accelerated warming in suburban areas flattened the temperature gradient between urban centers and suburbs. This dynamic expanded the spatial extent of elevated temperatures and reshaped the classic urban-to-rural UHI profile into a flatter form. In UingB areas, UHI intensification was jointly driven by increase in ISA, vegetation loss, and their interaction, while in UingG areas, EVI increases and a negative interaction together offset over half of the warming driven by urban expansion. These findings reveal that UHI evolution is not unidirectional but depends on localized urbanization and greening dynamics, offering pathways for strategic heat mitigation.

greening initiative↗
DOE OSTI2026

Importance of Spatially Continuous Urban Surface Properties in Urban‐Resolving Earth System Modeling

Accurate representation of urban properties and processes at higher resolutions in global modeling systems is essential for advancing our ability to capture the complexities of urban systems and informing effective resilience strategies. However, the prescription of coarse global-scale urban properties in most state-of-the-art Earth system models (ESMs) is limiting their potential for capturing urban signals as they advance toward kilometer-scale simulation capabilities. To bridge this gap in inadequate urban property representation and to advance urban-resolving Earth system modeling, this work integrates the newly-developed global 1 km-resolution facet-level urban surface property data set, U-Surf, into the land component of Community Earth System Model (CESM)—Community Terrestrial System Model (CTSM). The land-only CTSM simulations are validated against satellite measurements, ground-based urban weather stations, flux tower observations, and reanalysis data. Results demonstrate that the enhanced urban properties allow improved simulations of urban meteorology and surface energy fluxes compared to the default coarse-resolution categorical urban canopy parameters. Spatial scaling analysis reveals regime-dependent information loss during resolution aggregation, as well as substantial scale-dependent variations in urban surface energy flux representation. Furthermore, these findings have critical implications for coupled Earth system modeling when including the effect of land-atmosphere interaction. This work establishes a foundation for future urban-resolving kilometer-scale ESM development, which will enable systematic intra- and inter-city comparisons that inform urban adaptation strategies across diverse global urban environments.

Cheng, Yifan [University of Illinois Urbana-Champa↗
DOE OSTI2025

U-Surf: a global 1 km spatially continuous urban surface property dataset for kilometer-scale urban-resolving Earth system modeling

High-resolution urban climate modeling has faced substantial challenges due to the absence of a globally consistent, spatially continuous, and accurate dataset to represent the spatial heterogeneity of urban surfaces and their biophysical properties. This deficiency has long obstructed the development of urban-resolving Earth system models (ESMs) and ultra-high-resolution urban climate modeling, over large domains. Here, we present U-Surf, a first-of-its-kind 1 km resolution present-day (circa 2020) global continuous urban surface parameter dataset. Using the urban canopy model (UCM) in the Community Earth System Model as a base model for satisfying dataset requirements, U-Surf leverages the latest advances in remote sensing, machine learning, and cloud computing to provide the most relevant urban surface biophysical parameters, including radiative, morphological, and thermal properties, for UCMs at the facet and canopy level. Generated using a systematically unified workflow, U-Surf ensures internal consistency among key parameters, making it the first globally coherent urban canopy surface dataset. U-Surf significantly improves the representation of the urban land heterogeneity both within and across cities globally; provides essential, high-fidelity surface biophysical constraints to urban-resolving ESMs; enables detailed city-to-city comparisons across the globe; and supports next-generation kilometer-resolution Earth system modeling across scales. U-Surf parameters can be easily converted or adapted to various types of UCMs, such as those embedded in weather and regional climate models, as well as air quality models. The fundamental urban surface constraints provided by U-Surf can also be used as features for machine learning models and can have other broad-scale applications for socioeconomic, public health, and urban planning contexts. We expect U-Surf to advance the research frontier of urban system science, climate-sensitive urban design, and coupled human–Earth systems in the future. The dataset is publicly available at https://doi.org/10.5281/zenodo.11247598 (Cheng et al., 2024).

Cheng, Yifan [Univ. of Illinois at Urbana-Champaig↗
DOE OSTI2024

Urban morphology and urban water demand evolution in the Los Angeles region

Detailed description of the dataset sources used in this study, the experimental workflow, and plotting for the paper figures provided at the associated GitHub Meta Repo: https://github.com/IMMM-SFA/Ferencz_et_al_2024_ERL The future water demand projections from this study are hypothetical future water demands that reflect the population and urban land cover changes represented by the scenarios considered. The intent and emphasis of this work is investigating the interactions between population change, evolution of urban morphology, and water demand. These projections are not meant to be likely future demands for specific water providers or the LA region and should not be interpreted as such. The folders contain input and output data for each step of the "Recreate my Experiment" workflow described in the associated GitHub meta-repository as well as data used for plotting Figures for the paper that this dataset supports. Description of each folder's contents and use: Step_1a: Inputs to the associated python script provided on the GitHub repo. Step_1b: Inputs (downscaled population rasters) used by the associated python script provided on the GitHub repo. Original 1-km squared rasters that were downscaled also provided. Step_1c: Urban growth projection rasters corresponding to SSP3 and SSP5 population scenarios are provided in separate subfolders as well as the water provider boundaries used for analysis. Outputs of data processing also provided. Associated python script provided on GitHub. Step_1d: Description of Inputs used by the QGIS Model Builder GUI that automates geospatial processing and clipping the of the high-resolution 60 cm land cover data for each urban land class footprint within a defined polygon boundary. The Model Builder is provided on the GitHub repo and can be used by QGIS. The outputs of this step are in "Clipped Provider Hi Res Landcover". If the user wants to use The Model Builder for different regions of LA or to test our outputs, they will need to download the hi resolution landcover raster listed in the Readme and in Ref [2] of the GitHub Page. Step_1e: All necessary inputs to generate average monthly demand over the 2017-2021 period and the minimum and maximum demands over the 2014-2021 for each water provider. Associated python scripts are on GitHub. Step 2: Output data about land cover metrics (areas and fractions) for each urban land class for each water provider. Associated python script on GitHub. Uses outputs from Step 1d "Clipped Provider Hi Res Landcover" Step 3: Both the Inputs for and Outputs from the urban projection raster analysis Python script on GitHub. The inputs are urban land class rasters for specific SSP and zoning scenarios (low, medium, high) from Step 1c. The outputs are rasters of urban pixels that were converted to a higher land class and the number of land class units that changed (Values of 1, 2, or 3). For example, a value of 2 could be LC 21 -> 23 or LC 22 -> 24. These maps are label "intensification." The other outputs are "urban growth" rasters showing the conversion of non urban to urban land, which are indicated by pixel values of 1. These are used for the urban growth change maps in Figure 3. Step 4: Output projections of indoor and outdoor annual and monthly demands for each water provider for the average, minimum, and maximum monthly demand scenarios for each of the four urban growth scenarios (SSP3 med, SSP5 low, SSP5 med, and SSP5 high). The outputs also include metrics on each water provider used for the demand sensitivity analysis presented in Figure 8. Outputs from Step 4 are used for Figures 4 - 8 of the paper. Figures: This folder has data used for plotting Figures 1 through 5, and 8. Data for Figures 6 and 7 are sourced directly from folders associated with the Processing and Analysis Steps 1 - 4. The GitHub meta repository provides descriptions of how each figure was made and the associated plotting scripts used.

Los Angeles↗
DOE OSTI2026

Terrestrial laser scanning data (Levels 0 and 1) from Urban Biogeochemistry Pilot Project sites, Knoxville, Tennessee, Jul 2024 - Jul 2025

This data package contains data from terrestrial laser scanning (TLS) at five urban park sites in Knoxville, Tennessee, USA. All parks include open-grown and/or closed-canopy trees and mixed nearby land use. These study sites were established as part of the Urban Biogeochemistry Pilot Project, which has an overall goal of better understanding how hydrobiogeochemical cycling is altered within the human environment. These five sites represent a gradient of urbanization, and were instrumented to understand hydrological and biogeochemical cycling (e.g., soil moisture, soil physical properties and biogeochemistry, tree transpiration, species type). The TLS data archived here were collected to provide detailed, three-dimensional information about forest structure. Specifically, data were collected to allow tree- and stand-level characterization of woody structure and leaf area. TLS scans were placed to capture the area around trees with sap flow sensors, and as much of a 50 m radius area around the meteorological station as possible given site property limits. Derived products will allow upscaling of water content and transpiration data. This data package contains the following data: - High-level files document further details of the campaign and data package: 1_CampaignSummary.csv provides details about the campaign and study site, 2_ScanAreasDetail.csv provides details about each separate scan area (groups of scans post-processed into a single point cloud), 3_TerrestrialLidarSensor.csv provides further technical details about the Riegl VZ-400i TLS sensor, TLS_CSV_dd.csv is a CSV Data Dictionary providing information about the fields in CSV files following the ESS-DIVE CSV File Formatting Guidelines Reporting Format, TLS_flmd.csv is a File Level Metadata file providing information about each file in the data package following the ESS-DIVE File Level Metadata Reporting Format, and README.txt is a text file describing the overall project and file structure. - Level 0 data are the raw data (.PROJ folders) as recorded by the Riegl VZ-400i TLS instrument before scan co-registration and post-processing with the Riegl's proprietary RiSCAN PRO software, which requires a license. - Level 1 data contain post-processed, co-registered data from each scan area. The "PointClouds" folder for each scan area contains a .las file with 1 cm resolution point cloud data exported from RiSCAN PRO. These are the main files likely to be of interest to most users and can be further processed with any software capable of manipulating .las files (e.g. Python, R CloudCompare). The "Project Information" folder contains log files from post-processing in RiSCAN PRO that may be of interest to users who want to see detailed records of post-processing, including all PDF reports generated by RiSCAN PRO. The "ScanPositions" folder contains information about the final position of all TLS scans, after post-processing, in multiple formats. The file ScanPositions_*.csv provides final geo-referenced scan positions, and the file SOP_backup_*.csv can be used in RiSCAN PRO to restore the co-registered scan positions if users wish to re-process raw data (Level 0 .PROJ folders) with RiSCAN PRO software (e.g., subsample to a different resolution, exclude a certain scan position, or apply different filters on reflectance or deviation values) without redoing time-consuming co-registration steps.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Data and Scripts Associated with "Modeling Ecohydrological Responses of Vegetation to Urban Microclimates Using the E3SM Land Model"

This dataset supports the study of vegetation ecohydrological responses to urban microclimates using the land component of the Energy Exascale Earth System Model (ELM) at four urban sites in Knoxville, Tennessee, USA. It includes the model inputs, simulation outputs, and associated scripts for running ELM simulations and analyzing the resulting data. The Model_Inputs folder includes static surface data, satellite-derived phenology (i.e., leaf area index), and atmospheric forcing data used to drive ELM simulations. Detailed descriptions of these datasets are provided in Section 2.3.2 of the associated manuscript. The Model_Outputs folder contains simulation results for the baseline, treatment, and ensemble experiments. Outputs from the baseline and treatment simulations are provided as raw ELM NetCDF files. Because the raw outputs from the 4,000-member ensemble are prohibitively large, the ensemble results are provided as summarized CSV files, which also serve as the source data for Figure 5 of the associated manuscript. The Scripts folder contains three components: E3SM, the core codebase of the Energy Exascale Earth System Model (E3SM); elm-olmt, the Offline Land Model Testbed (OLMT) used to perform the simulations; and knoxville_elm, which contains the analysis scripts used to process model outputs and generate the figures and results presented in the associated manuscript. Additional information is provided in Scripts_readme.txt within the Scripts directory.

Lu, Xiaoman [ORNL] (ORCID:0000000306698780)↗
DOE OSTI2026

Soil Characterizations of Five Urban Sites in Knoxville, Tennessee. 2024-2025

This dataset includes soil characterizations from five urban parks in Knoxville, Tennessee, USA: Cumberland Estates Park (CEP), Socially Equal Energy Efficient Development (SEED), West View Park (WVP), Victor Ashe Park (VAP), and West Hills Park (WHP). The dataset consists of seven CSV files reporting the data from the measurements of gravimetric moisture content, pH, total carbon and nitrogen, soil texture, dissolved organic carbon and nitrogen, and microbial biomass carbon and nitrogen derived from these soil cores. During five separate sampling events conducted in 2024 and 2025, five soil cores were collected at each site within 3 meters of the remote soil monitoring equipment. In 2025, an additional three soil cores were collected adjacent to the monitoring equipment at each site to assess soil bulk density. This dataset is part of a larger study investigating the effects of soil moisture and plant evapotranspiration on ambient temperature and relative humidity across multiple urban parks in Knoxville, Tennessee.

Mayes, Melanie A [ORNL] (ORCID:0000000163689210)↗
DOE OSTI2026

Urban Parameters Los Angeles County 100m version 2

132 Urban parameters based on building physical dimensions and location were generated for the city of Los Angeles at 100m resolution using the NATURF model. To use the binary file with WRF, the binary file and the index file must be placed in their own directory in WRF_GEOG and accessed in the same way NUDAPT44 would be accessed. The kmz file can be visualized on Google Earth.

Sweet-Breu, Levi [Baylor University]↗
DOE OSTI2025

Meteorological Conditions in Urban Sites in Knoxville, USA. 2023

This 2023 dataset, which contains four csv files, reports hourly air temperature, wind speed, solar radiation, incident short wave, reflected short wave, incident long wave, emitted long wave, net radiation, and photosynthetically active radiation (PAR) data measured in urban parks in Knoxville, Tennessee, USA. The sites include Cumberland Estates Park (CEP), Socially Equal Energy Efficient Development (SEEED), West Hills (WHP), and Victor Ashe Park (VAP). Air temperature, wind speed and direction, and solar radiation data were obtained from a METER ATMOS 41 All-in-One Weather Station (Pullman, Washington, USA). Incident and emitted radiation (shortwave and longwave) measurements were made using an Apogee (Logan, Utah, USA) net radiometer (Model SN-500-SS). The SQ-521 Full-Spectrum Quantum sensor (Apogee Instruments, Inc) recorded the PAR. The measurement data for CEP, SEEED, and VAP started on July 1, 2023, while WHP started on August 26, 2023. All measurement ended on December 31, 2023. This work is a part of a larger study which investigates the impact of soil moisture and plant evapotranspiration on ambient temperature and relative humidity in several city parks in Knoxville, Tennessee.

Salvador, Christian [ORNL] (ORCID:0000000283287777↗
DOE OSTI2026

Meteorological Conditions in Urban Sites in Knoxville, USA. 2025

This 2025 dataset, which contains five csv files, reports hourly air temperature, wind speed, solar radiation, incident short wave, reflected short wave, incident long wave, emitted long wave, net radiation, and photosynthetically active radiation (PAR) data measured in urban parks in Knoxville, Tennessee, USA. The sites include Cumberland Estates Park (CEP), Socially Equal Energy Efficient Development (SEEED), West Hills Park (WHP), West View Park (WVP), and Victor Ashe Park (VAP). Air temperature, wind speed, and solar radiation data were obtained from a METER ATMOS 41 All-in-One Weather Station (Pullman, Washington, USA). Incident and emitted radiation (shortwave and longwave) measurements were made using an Apogee (Logan, Utah, USA) net radiometer (Model SN-500-SS). The SQ-521 Full-Spectrum Quantum sensor (Apogee Instruments, Inc) recorded the PAR. The measurement data for all the sites started on January 1, 2025. and all measurement ended on December 31, 2025. This work is a part of a larger study which investigates the impact of soil moisture and plant evapotranspiration on ambient temperature and relative humidity in several city parks in Knoxville, Tennessee.

Salvador, Christian [ORNL] (ORCID:0000000283287777↗
DOE OSTI2026

Urban Parameters Washington, DC 100m

132 Urban parameters based on building physical dimensions and location were generated for Washington, District of Columbia at 100m resolution using the NATURF model. To use the binary file with WRF, the binary file and the index file must be placed in their own directory in WRF_GEOG and accessed in the same way NUDAPT44 would be accessed.

Dumas, Melissa [ORNL] (ORCID:0000000233190846)↗
DOE OSTI2026

NATURF: Urban Building Parameters for Chicago, Illinois, USA at a 100m resolution

132 Urban parameters based on building physical dimensions and location were generated for the city of Chicago at 100m resolution using the NATURF model. To use the binary file with WRF, the binary file and the index file must be placed in their own directory in WRF_GEOG and accessed in the same way NUDAPT44 would be accessed.

Vernon, Chris R [Pacific Northwest National Labora↗
DOE OSTI2024

Sensitivity of mesoscale modeling to urban morphological feature inputs and implications for characterizing urban sustainability

We examine the differences in meteorological output from the Weather Research and Forecasting (WRF) model run at 270 m horizontal resolution using 10 m, 100 m and 1 km resolution 3D neighborhood morphological inputs and with no morphological inputs. We find that the spatial variability in temperature, humidity, and other meteorological variables across the city can vary with the resolution and the coverage of the 3D urban morphological input, and that larger differences occur between simulations run without 3D morphological input and those run with some type of 3D morphology. We also find that the inclusion of input-building-defined roughness length calculations would improve simulation results further. We show that these inputs produce different patterns of heat wave spatial heterogeneity across the city of Washington, DC. These findings suggest that understanding neighborhood level urban sustainability under extreme heat waves, especially for vulnerable neighborhoods, requires attention to the representation of surface terrain in numerical weather models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Demographic Microsimulator for Integrated Urban Systems: Adapting Panel Survey of Income Dynamics to Capture the Continuum of Life

Agent-based models (ABMs) in transportation modeling simulate activity and travel decisions at the disaggregate level of households and individuals. To do this, ABMs require detailed and realistic information on agents’ socioeconomic and demographic characteristics. Various synthetic population generators have been proposed to address this need. However, most of those currently in practice are cross-sectional in nature and do not account for the dynamics within households and individuals as they progress through life events over time. This is a major shortcoming, as literature has shown that transportation decisions are affected by the transition between and co-occurrence of life cycle events. While some demographic evolution simulators have been proposed to address this issue, they are developed using cross-sectional data and capture only a small set of life cycle events and their interdependence. Addressing these drawbacks, we propose a demographic microsimulator (DEMOS) that captures the “continuum of life” by considering a range of household- and individual-level life cycle events. DEMOS is developed using the Panel Survey of Income Dynamics, one of the world’s longest-running longitudinal surveys. The DEMOS submodels consider key life cycle events that are influenced by agents’ demographic variables. DEMOS is applied to evolve the population of the San Francisco Bay Area over a 9-year horizon. Results demonstrate how DEMOS generates life trajectories and how DEMOS outputs match the observed demographic trends. DEMOS is expected to enable longitudinal analysis in the context of ABMs and expand ABMs analyses relating to dynamic processes such as household-level vehicle transactions.

Demographic evolution↗