
==== Front
Geohealth
Geohealth
10.1002/(ISSN)2471-1403
GH2
GeoHealth
2471-1403
John Wiley and Sons Inc. Hoboken

10.1029/2024GH001079
GH2559
2024GH001079
Geospatial Data Applications for Environmental Justice
Geospatial Data Applications for Environmental Justice
Atmospheric Composition and Structure
Pollution: Urban and Regional
Atmospheric Composition and Structure
Aerosols and Particles
Air/Sea Constituent Fluxes
Volcanic Effects
Biogeosciences
Urban Systems
Biogeosciences
Climate Dynamics
Modeling
Pollution: Urban, Regional and Global
Computational Geophysics
Modeling
Numerical Solutions
Cryosphere
Avalanches
Mass Balance
Modeling
Geodesy and Gravity
Mass Balance
Ocean Monitoring with Geodetic Techniques
Ocean/Earth/atmosphere/hydrosphere/cryosphere interactions
Global Change from Geodesy
Global Change
Abrupt/Rapid Climate Change
Climate Variability
Climate Dynamics
Earth System Modeling
Impacts of Global Change
Land/Atmosphere Interactions
Oceans
Regional Climate Change
Sea Level Change
Solid Earth
Water Cycles
Hydrology
Climate Impacts
Hydrological Cycles and Budgets
Land/Atmosphere Interactions
Modeling
Informatics
Modeling
Marine Geology and Geophysics
Gravity and Isostasy
Atmospheric Processes
Climate Change and Variability
Atmospheric Processes
Climatology
General Circulation
Land/Atmosphere Interactions
Ocean/Atmosphere Interactions
Regional Modeling
Theoretical Modeling
Oceanography: General
Climate and Interannual Variability
Marine Pollution
Numerical Modeling
Natural Hazards
Megacities and Urban Environment
Vulnerability
Natural Hazards
Atmospheric
Geological
Oceanic
Physical Modeling
Climate Impact
Risk
Disaster Risk Analysis and Assessment
Oceanography: Physical
Air/Sea Interactions
Decadal Ocean Variability
General Circulation
Ocean influence of Earth rotation
Sea Level: Variations and Mean
Surface Waves and Tides
Tsunamis and Storm Surges
Oceanography: Biological and Chemical
Aerosols
Paleoceanography
Abrupt/Rapid Climate Change
Aerosols
Policy Sciences
Benefit‐cost Analysis
Radio Science
Radio Oceanography
Seismology
Earthquake Ground Motions and Engineering Seismology
Volcano Seismology
Volcanology
Volcano/Climate Interactions
Atmospheric Effects
Volcano Monitoring
Effusive Volcanism
Mud Volcanism
Explosive Volcanism
Volcanic Hazards and Risks
Research Article
Research Article
Synergy of Urban Heat, Pollution, and Social Vulnerability in One of America's Most Rapidly Growing Cities: Houston, We Have a Problem
Blackford et al.
Blackford Andrew https://orcid.org/0000-0002-5407-5572
1
Cowan Trent 1
Nair Udaysankar https://orcid.org/0000-0002-3724-9630
1 nairu@uah.edu

Phillips Christopher https://orcid.org/0000-0001-7539-7362
2
Kaulfus Aaron https://orcid.org/0000-0002-8319-1126
3
Freitag Brian https://orcid.org/0000-0003-3642-7041
3
1 Department of Atmospheric and Earth Science The University of Alabama in Huntsville Huntsville AL USA
2 Earth System Science Center The University of Alabama in Huntsville Huntsville AL USA
3 National Aeronautical and Space Administration Marshall Space Flight Center Huntsville AL USA
* Correspondence to:
U. Nair,
nairu@uah.edu

04 9 2024
9 2024
8 9 10.1002/gh2.v8.9 e2024GH00107916 7 2024
22 4 2024
01 8 2024
© 2024 The Author(s). GeoHealth published by Wiley Periodicals LLC on behalf of American Geophysical Union.
https://creativecommons.org/licenses/by/4.0/ This is an open access article under the terms of the http://creativecommons.org/licenses/by/4.0/ License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Abstract

During the first two decades of the twenty‐first century, we analyze the expansion of urban land cover, urban heat island (UHI), and urban pollution island (UPI) in the Houston Metropolitan Area (HMA) using land cover classifications derived from Landsat and land/aerosol products from NASA’s Moderate Resolution Imaging Spectroradiometer. Our approach involves both direct utilization and fusion with in situ observations for a comprehensive characterization. We also examined how social vulnerability within the HMA changed during the study period and whether the synergy of UHI, UPI, and social vulnerability enhances environmental inequalities. We found that urban land cover within the HMA increased by 1,345.09 km2 and is accompanied by a 171.92 (73.93) % expansion of the daytime (nighttime) UHI. While the UPI experienced an overall reduction in particulate pollution, the magnitude of change is smaller compared to the surroundings. Further, the UPI showed localized enhancement in particulate pollution caused by increases in vehicular traffic. Our analysis found that the social vulnerability of the HMA urban regions increased during the study period. Overall, we found that the urban growth during the first two decades of the twenty‐first century resulted in a synergy of UHI, UPI, and social vulnerability, causing an increase in environmental inequalities within the HMA.

Key Points

The rapid urbanization of Houston expands the urban heat and pollution islands, raising exposure to high temperatures and poor air quality

Communities with highest social vulnerability are most impacted by the adverse effects of urban heat and pollution islands

The Clean Air Act’s efficacy is locally diminished by Houston’s urban pollution island

Houston
pollution
heat stress
urban heat island
urban pollution island
social vulnerability
NASA MSFC80MSFC22M0004 source-schema-version-number2.0
cover-dateSeptember 2024
details-of-publishers-convertorConverter:WILEY_ML3GV2_TO_JATSPMC version:6.4.8 mode:remove_FC converted:04.09.2024
Blackford, A. , Cowan, T. , Nair, U. , Phillips, C. , Kaulfus, A. , & Freitag, B. (2024). Synergy of urban heat, pollution, and social vulnerability in one of America's most rapidly growing cities: Houston, we have a problem. GeoHealth, 8 , e2024GH001079. 10.1029/2024GH001079
==== Body
pmc1 Introduction

The Houston Metropolitan Area (HMA), Texas, one of the top 10 most populous metropolitan areas in the United States (U.S.), is ethnically and racially diverse but also faces severe racial and economic segregation (Emerson et al., 2012). Approximately 24% and 37% of HMA's high‐income and low‐income households, respectively, are located within corresponding neighborhoods linked to historic redlining practices (Lane et al., 2022; Nicholson, 2017). Such segregation contributes to environmental inequities, including heightened heat stress, pollution, flood risk, and reduced green space access (Bullard, 1983; Fry & Taylor, 2012; Logan et al., 2021; Waren, 2013). Prior work (Ulpiani, 2021) has found that urban regions have increased air temperature and pollution compared to rural regions, termed Urban Heat Islands (UHI) and Urban Pollution islands (UPI) respectively. Additionally, social vulnerability refers to socioeconomic factors that increase the risk of adverse health effects from environmental conditions (Mah et al., 2023). Understanding synergies between the UHI, UPI, and social vulnerability are therefore important for equitable and healthy urban growth (H. Li et al., 2018; Logan et al., 2021; Streutker, 2003; Ulpiani, 2021). Additionally, urban growth leads to the expansion of the wildland‐urban interface, with emissions from prescribed burning around the HMA adding to urban pollution (Nowak & Walton, 2005).

Prior studies have examined UHI and UPI associated with the HMA (Moser et al., 2017; Streutker, 2002, 2003). However, the impact of the last two decades' population and urban land cover expansion on the UHI, UPI, social vulnerability, and their interactions remain unexamined. A major hurdle is the inadequacy of sparse meteorological and air quality networks for monitoring UHI and UPI growth. This concern is heightened by the HMA's complex spatial pollution variability, influenced by approximately 400 chemical manufacturing units, two of the U.S.'s biggest oil refineries, daily vehicular traffic exceeding 140 million miles a day, small urban sources such as gas stations, restaurants, dry cleaners, and construction, and episodic smoke transport from fires in Mexico, Central America, and surrounding states (O’Dell et al., 2020; Sexton & Linder, 2015; Sexton et al., 2007). Notably, pollution sensors on Google Street View cars in the HMA reveal fine‐scale maxima of air pollution undetected by traditional air quality monitoring networks (Miller et al., 2020).

Implications of Miller et al. (2020) are quite substantial, as the United Nations (2022) projects that by 2050, 87% of the American population and 68% of the world population will live in urbanized locations. This study thus addresses the above mentioned limitations by utilizing the National Aeronautical Space Administration (NASA) Terra and Aqua Earth‐observing platforms' aerosol and land temperature products which are combined with situ observations and socioeconomic data sets to quantify the UHI and UPI (specifically particulate pollution) growth over the HMA. This study aims to: (a) characterize urban growth that has occurred over the HMA during the last two decades, (b) characterize associated growth of the UHI and UPI, and (c) analyze how the combined effects of the UHI and UPI contribute to environmental inequalities in the region.

2 Study Area

The U.S. census‐designated area of Houston ‐The Woodlands‐Sugarland Metropolitan Area (HMA) is situated in southeastern Texas along the Gulf of Mexico (GoM). It includes Harris County, containing the city of Houston, and the surrounding counties—Galveston, Chambers, Liberty, Montgomery, Waller, Austin, Fort Bend, and Brazoria. The HMA approximately occupies 24,459 km2, surpassing the size of New Jersey. Houston, the fourth most populous U.S. city, had a city‐limits population of nearly 2.3 million and a metropolitan area population of 6.5 million as of 2021 (United Nations, 2022; U.S. Census Bureau, 2022).

3 Data and Methodology

The data sets and analysis methodology used in this study are described in Sections 3.1, 3.2, 3.3, 3.4, 3.5, 3.6.

3.1 Satellite Observations

To quantify multi‐decadal urban land cover growth in the HMA, we utilize the National Land Cover Database (NLCD) land cover classification for the contiguous U.S. (CONUS). Derived from Landsat satellite observations at 30 m spatial resolution, (C. G. Homer et al., 2007; C. H. Homer et al., 2012; C. Homer et al., 2015) NLCD is available at 2–3‐year intervals from 2001 to 2021. NLCD classifies locations on a 30 m CONUS grid into one of 34 land cover classes, with high‐density, medium‐density, low‐density, and open‐space urban classes relevant to this study. Our analysis uses the 2001 and 2019 NLCD data sets to examine and quantify multi‐decadal changes in spatial patterns and the extent of urban land cover.

We analyze alterations in Land Surface Temperature (LST), vegetation cover, and atmospheric aerosol loading using land and aerosol products derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor on NASA’s Terra and Aqua platforms. Both satellites follow sun‐synchronous orbits, passing over the HMA daily at around 10:30 a.m. and 1:30 p.m. local time. For consistency and data availability throughout the first two decades (Aqua operations commenced midway through 2002) and to minimize cloud contamination, we only utilized land products from Terra. The eight‐day composite day and night LST (MOD11A2) at 1 km spatial resolution with retrieval errors of 2.0 ± 0.5 K (Wan, 2014) is used to investigate UHI changes. Changes in urban greenness were examined using the 16‐day composite Normalized Difference Vegetation Index (NDVI) product at 250 m spatial resolution (MOD13Q1). NDVI is calculated as the ratio of the difference between near‐infrared and red reflectance channels normalized by the sum of reflectance values of the same channels. With values ranging from −1 to 1, NDVI signifies extremes of bare ground (zero green vegetation cover) and full coverage by green vegetation at a given location (Didan et al., 2015).

The Aerosol Optical Depth (AOD) product utilized is a daily 0.55 μm Terra and Aqua composite with 1 km spatial resolution, (MCD19A2). The AOD is calculated using the multi‐angle implementation of the atmospheric correction (MAIAC) algorithm (Lyapustin & Wang, 2018) and reflects the total aerosol mass within the atmospheric column over a specific location. To ensure data consistency and avoid potential anomalies from Aqua observations available only beginning in mid‐2002, the AOD analysis focuses solely on data from 2003 to 2019.

All of the satellite data were accessed and analyzed using Google Earth Engine (GEE). Note that GEE is a global cloud‐computing platform that provides access to satellite imagery and derived products acquired over the last 50 years (Gorelick et al., 2017).

3.2 Surface Meteorological and Air Quality Observations

We analyzed surface temperature, humidity, and particulate pollution trends in the HMA using in situ meteorological and air quality observations. Ten meteorological sites‐ four in the Houston city limits, three in the suburbs, and three on the outskirts‐ were also used for statistical regression modeling to develop high‐resolution spatial patterns of monthly maximum air temperature. Data were obtained from NOAA’s National Centers for Environmental Information (NCEI), with the station details in Table S5 of the Supporting Information S1, ensuring consistency in station locations.

For particulate air quality, PM2.5 concentrations from six EPA monitoring stations within the HMA were used. PM2.5 refers to inhalable particulates with a diameter of less than 2.5 μm, posing the greatest health risk.

3.3 Social Vulnerability Index

Social vulnerability, encompassing socioeconomic and demographic factors affecting community resilience to external stresses, (Flanagan et al., 2020; Mah et al., 2023) can be quantified through a social vulnerability index (SVI). The CDC’s SVI data set, derived from U.S. Census Bureau's census tract data, incorporates various factors such as socioeconomic status, household composition, disability, minority status, housing, and transportation attributes (Flanagan et al., 2020). The overall SVI ranking used here represents the nationwide percentile ranking of social vulnerability for a census tract. An SVI value of 1 indicates that all other census tracts have equal or lower social vulnerability, while a value of 0 means none have equal or lower social vulnerability compared to the considered census tract. This study employs SVI to assess the societal impacts of the UHI and UPI in terms of environmental inequalities.

3.4 Satellite Data Analysis of Urban Growth, LST, NDVI, and AOD

We analyzed NLCD LULC classifications for 2001 and 2019, quantifying transitions between high, medium, low‐density, and open‐space urban classes. Urban land cover masks for both years were used to calculate statistics for LST, NDVI, AOD, SVI, maximum air temperature, and PM2.5 separately in urban and non‐urban regions. To study the influence of urban growth, a mask identifying changes in urban extent between 2001 and 2019 was derived.

To assess changes in LST, NDVI, and AOD, we generated decadal average spatial maps for 2000–2009 and 2010–2019. This involved the computation of monthly averages of each variable (ρ‾m(λ,ϕ)) for each pixel across the HMA domain: (1) ρ‾m(λ,ϕ)=∑i=1Nρ(xy,i)N

where ρ(λ, ϕ, i) is the gridded satellite‐derived value of the variable at a given latitude (x) and longitude (y), i is a specific instance of observation/retrieval, and N the total number of instances. Subscript m refers to the monthly mean of the variable. The variable average (ρ‾d(λ,ϕ)) for the first and second decades, where d refers to the decadal mean, is then computed by averaging the monthly averages for the periods of 2000–2009 and 2010–2019, except for AOD. For the first decade, the AOD average is computed using 2003–2009 data to account for the unavailability of Aqua observations before mid‐2002. Note that we utilized decadal average spatial maps, as only persistent features will be retained in long‐term averages, thereby increasing the confidence in change analysis.

To quantify temporal changes in urban regions, we spatially averaged decadal mean maps of day and night LST, NDVI, and AOD over census tracts in the HMA region with more than 50% urban land cover (referred to as urban census tracts). Urban census tracts for the 2000–2009 period were defined using 2000 census tracts and 2001 NLCD classifications, while those for the 2010–2019 period used 2020 census tracts and 2019 NLCD classifications. Trend analysis for the 2000–2019 period is conducted using time series created from spatial averages of monthly mean LST, NDVI, and AOD maps for urban census tracts. The 2020 census tracts and 2019 NLCD classification were used for this computation to ensure a consistent average for meaningful trend analysis. The seasonal Mann‐Kendall test determined statistically significant monotonic trends, and the Thiel‐Sen slope quantified the trend (Kendall, 1948; Mann, 1945; Sen, 1968; Thiel, 1950). Note that for the processing of LST for computing decadal averages and time series, the methodology of L. Hu and Brunsell (2013) was applied to minimize cloud contamination in our results.

3.5 Analysis of Surface Meteorological and Air Pollution Observations

To link satellite analysis of LST with near‐surface air temperature and particulate pollution, we performed a time series analysis of monthly averaged values for daily maximum air temperature and PM2.5 concentrations. The time series was examined for monotonic trends using the Mann‐Kendall seasonal test, and the trends were quantified using the Sen slope.

We employed statistical regression modeling to fuse surface observations with satellite data, creating high‐resolution estimates of monthly average maximum air temperature. Multiple linear regression models were developed for monthly average maximum air temperature (T max), using the following predictors chosen based on prior research (de Souza et al., 2022; Yuvaraj, 2020): monthly average satellite‐derived LST and season (expressed categorically as DJF‐1, MAM‐2, JJA‐3, SON‐4). Regression modeling utilized monthly averaged maximum temperature from 10 meteorological stations, paired with corresponding monthly averaged daytime LST. The model was trained on 80% of the data and tested on the remaining 20%. The T max model was then applied to generate high‐resolution maps of decadal seasonal average maximum air temperatures for the first and second decades and analyze the impact of urban growth.

We utilized high spatial resolution maps of surface PM2.5 derived from a data fusion approach. Instead of developing models, we employed the surface PM2.5 data set of Van Donkelaar et al. (2021), which integrates MAIAC AOD, other NASA aerosol products, surface observations, and chemical transport model outputs. This validated data set offers monthly average PM2.5 maps (PMavg), enabling the estimation of changes in PM2.5 concentrations due to urban growth.

3.6 Social Vulnerability Analysis

We compared SVI spatial patterns between 2000 and 2020 to assess how urban growth affects the social vulnerability of the HMA. Census tracts were grouped into four categories based on SVI values (0.0–0.25, 0.25–0.50, 0.50–0.75, and 0.75–1.0). For each group, decadal average LST, NDVI, AOD, maximum air temperature, and surface PM2.5 were computed to evaluate the combined effects of UHI, UPI, and social vulnerability.

4 Results

4.1 Multi‐Decadal LULC Change Impacts in the Houston Metropolitan Area

Between 2001 and 2019, the HMA witnessed a substantial 1,345 km2 (22.6%) increase in urban land cover (Figure 2, Table 1), equivalent to adding two times the entire urban area of New York City. This growth primarily resulted from converting agricultural and woodland classes to urban and built‐up land cover. Medium density urban land cover showed the highest growth in both area and percentage, followed by high‐density urban land cover, while the open space urban category exhibited the least growth. The spatial pattern of urban growth expanded in all directions, with higher growth in the northern, northwestern, and southwestern directions. The HMA consistently experienced positive urban growth, with annual percent increases ranging from 1.8% to 5.1% between 2001 and 2019, reflecting a steady yet occasionally accelerated urbanization trend.

Figure 1 (a) The spatial extent of the study area for this paper, with the red box in the inset (b) showing the location of the study area within the state of Texas.

Figure 2 2001 (a) and 2019 (b) NLCD LULC across the HMA, and the new urbanization footprint (in red) from 2001 to 2019 (c). Reds in (a) and (b) indicate urbanized classifications, where DLI, DMI, and DHI are low, medium, and high intensity urban respectively. The complete list of abbreviation meanings can be found in Supporting Information S1.

Table 1 Mean Values of Selected MODIS Variables As Well As Surface T max and PM2.5 Across the HMA Urban Census Tracts, Isolated by SVI Bins

	2000–2009 (2010–2019)	2010s–2000s difference	
SVI bin	LST day	LST night	NDVI	AOD	T max	PM2.5	LST day	LST night	NDVI	AOD	T max	PM2.5	
0.0–0.25	300.48 (300.97)	289.83 (289.94)	0.49 (0.49)	0.15 (0.14)	304.71 (304.99)	11.39 (9.02)	1.023	0.456	−0.021	−0.011	0.407	−2.323	
0.25–0.50	300.12 (300.78)	289.40 (289.80)	0.52 (0.50)	0.14 (0.14)	304.37 (304.66)	11.03 (8.72)	0.816	0.404	−0.018	−0.005	0.341	−2.306	
0.50–0.75	300.39 (300.92)	289.73 (289.84)	0.49 (0.49)	0.15 (0.14)	304.61 (304.74)	11.25 (8.73)	0.750	0.370	−0.016	−0.010	0.318	−2.298	
0.75–1.0	300.84 (301.83)	289.84 (290.28)	0.46 (0.45)	0.15 (0.15)	304.83 (305.41)	11.30 (9.36)	0.837	0.344	−0.019	−0.002	0.327	−2.236	

Figure 3 shows decadal averages of daytime and nighttime LST as well as AOD. Spatial patterns of daytime and nighttime LST align with urban land cover growth in Figures 2c and 3a–3d, exhibiting higher values over urban areas. Decadal‐averaged daytime LST for urban majority census tracks surpassed non‐urban tracts during both 2000–2009 and 2010–2019 (Table S3 in Supporting Information S1). The contrast in decadal average daytime (nighttime) LST between urban and non‐urban tracts increased from ∼4.2 (1.8) K to 4.5 (1.9) K between the two decades. Urban tracks’ decadal‐averaged daytime (nighttime) LST increased by ∼0.86 (0.40) K from 2000 to 2009 to 2010–2019, with a mean difference of ∼0.30 (∼0.13). We used the decadal average daytime LST for urban census tracts during the 2000–2009 period as a threshold (∼303 K) to define the urban LST footprint. Using this threshold, the daytime urban LST footprint expanded from 527 km2 (2000–2009) to 1,433 km2 (2010–2019). Figure 4 reveals the spatial difference in decadal average AOD and daytime and nighttime LST between the 2010–2019 and 2000–2009 periods. Both the day and night LST difference patterns align with the pattern of differences in urban land cover between the two decades (Figure 2c), although less pronounced in the case of nighttime LST. Decadal average analysis of NDVI and AOD also reveals differences in mean values between urban and non‐urban regions (Table S3 in Supporting Information S1). Decadal average NDVI for non‐urban census tracts remains stable between 2000–2009 and 2010–2019, but urban census tracts show a decrease of ∼0.006. The spatial pattern of NDVI changes shows maximum differences (Figure S3 in Supporting Information S1) in areas of urban land cover growth (Figure 2c). Decadal average AOD also show higher mean values over urban census tracts compared to non‐urban tracts in both decades. However, the decrease in mean decadal AOD between the two decades is less for urban census tracts than for non‐urban tracts. While a mean difference in AOD of −0.004 is found in the newly urbanized regions, the spatial pattern of decadal average AOD do show small, localized increases in pockets over such areas (Figure 4c).

Figure 3 Decadal‐averaged spatial distributions of MODIS‐derived monthly mean Daytime LST (a, b), Nighttime LST (c, d), and AOD (e, f). The 2000–2009 decade is displayed in the top row, and the 2010–2019 decade is displayed in the bottom row. Note that AOD data begins in 2002, while LST Day and Night both extend back to 2000.

Figure 4 (a) Daytime LST, (b) Nighttime LST, and (c) AOD 2010s decadal average minus 2000s decadal average over the HMA.

Utilizing the seasonal Mann‐Kendall test, we assessed statistically significant monotonic trends in time series of monthly mean daytime and nighttime LST, NDVI, and AOD spatially averaged for urban census tracts (Figure S2 in Supporting Information S1). For the 2000–2009 period, statistically significant monotonic trends were observed only for NDVI and nighttime LST when analyzed separately for each decade (Table S4 in Supporting Information S1). However, when considering the entire multi‐decadal period, all variables displayed statistically significant monotonic trends. Daytime and nighttime LST exhibited increasing trends, with the magnitude of Sen slope indicating a higher daytime trend compared to nighttime. In contrast, both NDVI and AOD showed decreasing trends over the multi‐decadal study period.

4.2 Impact of Urbanization on Daytime Patterns of Near‐Surface Air Temperature and PM2.5

We used in situ observations to assess if trends in satellite‐derived LST and AOD correspond to changes in near‐surface air temperature and PM2.5 concentrations (Table S6, Figure S4 in Supporting Information S1). Therefore, we examined time series trends of maximum and minimum air temperatures and average PM2.5 concentration observations and how they vary spatially.

We found that most stations exhibited statistically significant, monotonically increasing trends in maximum and minimum surface air temperature (Table S6 in Supporting Information S1). Minimum temperature trends generally exceeded maximum temperature trends, with higher values observed in regions of urban growth. Urban stations also showed greater trends, with maximum temperature trends higher by ∼0.013 K year−1 and minimum temperature trends higher by ∼0.040 K year−1. The Galveston (KGLS) and Hooks Memorial Airport (KDWH) localities, coinciding with observed urban growth, displayed the largest trends in maximum and minimum temperatures, respectively. Houston Hobby Int’l Airport (KHOU), within the urban core with minimal land cover changes, exhibited lower Sen slope estimates.

All three PM2.5 monitoring stations that reported at least 10 years of continuous data (EPA stations 482010058, 482011039, and 482011035, where 48201 is the code for Harris County, Texas; Figure 1) exhibited statistically significant decreasing trends, indicating a general reduction in particulate pollution. Station 0058, surrounded by the least impervious surface, showed the lowest trend magnitude, while station 1035, within the highest impervious surface, displayed the highest trend. The spatial pattern of these trends reflects an overall reduction in particulate pollution influenced by Clean Air Act regulations, which has been shown in prior studies (Environmental Protection Agency, 2019). However, localized urban pollution sources weaken the effectiveness of these regulations, and this effectiveness also varies based on the degree of urbanization. This limited analysis, in line with spatial patterns in decadal average AOD (Figures 3e, 3f, and 4c), underscores the necessity for a more extensive surface PM2.5 network for broader validation.

We employed a multiple linear regression model for monthly average maximum air temperature (T max) to infer high‐resolution spatial patterns. The constructed T max model (Table S7 in Supporting Information S1) yielded an R2 which value of 0.844 and an RMSE of 1.62 K, demonstrating performance comparable to or better than similar models in prior studies (de Souza et al., 2022; Yuvaraj, 2020). Utilizing the T max model, we generated monthly maps of average maximum temperature for the HMA region. These maps were then used to compute decadal average seasonal and annual maximum temperature maps for both the 2000–2009 and 2010–2019 decades (Figure 5).

Figure 5 Decadal‐averaged spatial distributions of predicted T max (a, b), produced from multi‐variate ordinary least squares linear regression and reported in K, and surface PM2.5 (d, e). Decadal differences of T max (c) and PM2.5 (f) are also displayed.

Spatial maps of maximum temperature reveal contrasts between urban regions and their surroundings, intensifying during the second decade. The difference in decadal average maximum temperature between the two decades (Figure 5) highlights a robust warming signal along the ring of new urbanization from 2001 to 2019 (Figure 2c), aligning with earlier time series trend analysis. Although decadal differences can exceed 2 K in specific locations, the mean difference in decadal average maximum temperature for urban regions is approximately ∼0.25 K. This enhanced warming signal associated with urbanization is present throughout all seasons (Figure S6 in Supporting Information S1), with the most notable contrast occurring during the spring season.

We conducted a similar analysis using downscaled surface PM2.5, combining in situ observations with satellite data sets. Despite the absence of a notable trend in AOD (Tables S3 and S4 in Supporting Information S1), in situ sensors indicated some decreasing trends in PM2.5 (Figure S5 in Supporting Information S1). While the downscaled analysis was expected to mirror the trends observed by in situ sensors, our findings reveal contrasting implications. First, decadal average patterns of downscaled PM2.5 show higher values over urban regions compared to the surroundings. Mean PM2.5 concentrations in urban locations exceeded the primary standards set by the EPA (12 µgm−3) in the first decade but not in the second decade. Differences in decadal average surface PM2.5 indicate the least change in concentrations along the ring of new urbanization throughout the study period (Figures 2c and 5f). Notably, the decadal difference in PM2.5 falls below −3 µgm−3 at certain locations, but the mean decadal difference is −2.34 µgm−3. Second, this discrepancy between AOD and PM2.5 trends suggests that column‐total AOD, which includes contributions from atmospheric transport, has different controlling factors compared to in situ observations. Thus, PM2.5 trends may be more reflective of the ongoing urbanization process as they include less transport from other regions.

4.3 Socioeconomic Impacts of Rapid Urbanization

Finally, we assessed changes in social vulnerability within the HMA over the multi‐decadal study period and investigated how interactions between heat, pollution, and social vulnerabilities contribute to environmental inequalities. Census tracts were categorized based on SVI percentile rankings: 0.0–0.25 (least socially vulnerable), 0.25–0.50, 0.50–0.75, and 0.75–1.0 (most socially vulnerable). The spatial distribution of the SVI revealed a significant increase in social vulnerability within the urban regions of the HMA in the last two decades (Figure 6). The largest driver of this change in the main urban core of the HMA was found to be socioeconomic status, followed closely by household characteristics. The most socially vulnerable urban census tracts in the HMA increased by 629 tracts during this period.

Figure 6 2000 (a) and 2020 (b) SVI ranking (0‐1) for HMA census tracts.

Spatial patterns of decadal average day and nighttime LST, AOD, maximum surface air temperature, and surface PM2.5 concentrations reveal census tracts with the highest SVI are situated in areas where populations are most susceptible to heat and pollution exposure (Figures 3 and 5). The mean of decadal average day and nighttime LST, AOD, maximum air temperature, and PM2.5 indicates higher values for the most vulnerable census tracts during both decades (Table 1).

The change in LST and surface air temperature between the two decades is most significant for the least socially vulnerable census tracts, possibly due to increased urbanization in suburban regions of Houston (Figure 2c). Although NDVI also shows the highest decrease for the least vulnerable census tracts, the mean decadal average NDVI is lowest for census tracts with the most social vulnerability, suggesting reduced access to green spaces in these areas. Mean decadal averages of both AOD and surface PM2.5 indicate a decrease in particulate matter air quality across all census tract categories. However, the decrease between the two decades is least for the most socially vulnerable census tracts and higher for census tracts with the least social vulnerability. Although this analysis was conducted for specific SVI bins across the entire HMA, focusing on a newly urbanized census tract to another with identical vulnerability but no urbanization, the newly urbanized tract experiences the most significant increase in heat stress.

5 Uncertainty and Limitations

The following methodological limitations need to be considered in the interpretation of the study findings. In addition to errors inherent in the retrieval of MODIS land and aerosol products, retrieval is only valid under cloud‐free conditions. While long‐term averaging mitigates the impact of retrieval errors, cloud cover modulates the number of valid retrievals at any given geographical location. Thus, the number of retrievals used to construct decadal averages varies at each geographical location. Further, spatial variability of cloud cover itself can be affected by urban growth (Theeuwes et al., 2019; Vo et al., 2023), thereby impacting computed statistics of the variables utilized in the study (and thus spatial patterns in decadal averages and trend analysis). Given these limitations of the study, future work may include relatively high‐resolution numerical model experiments incorporating scenarios of realistic LULC changes over the multidecadal period.

6 Discussion and Conclusions

Our analysis shows that the HMA, one of the largest metropolitan areas in the United States, has experienced rapid urbanization over the past two decades, adding ∼1,345 km2 of developed land cover, which is nearly twice as large as the area of New York City. This urban growth is accompanied by both UHI and UPI growth as well as heightened social vulnerability. We found that both the day and nighttime UHI LST footprint also expanded with urban growth and the decadal average values showed increases even within regions that were urbanized before the study period, with corresponding increases in UHI as reflected by increases in maximum and minimum air temperatures (Table 1, Table S3 in Supporting Information S1).

While the UPI associated with the HMA experienced an overall reduction in particulate pollution, the magnitude of the decrease is less when compared to rural areas. Over the urban areas, decadal changes in AOD and surface PM2.5 both show higher spatial variability compared to the rural surrounding areas, and exhibit a more homogenous and drastic reduction during the second decade of the study. While our analysis shows that the Clean Air Act of 1990, with subsequent evaluations having occurred most recently in 2012, resulted in an overall decrease in particulate pollution within the UPI. However, the increased volume of vehicular traffic (Figure S7 in Supporting Information S1) and other sources associated with urban growth appear to locally counteract the overall reduction of emissions from individual sources. These findings align with similar studies (Lim et al., 2020; Miller et al., 2020; Southerland et al., 2022), specifically Miller et al. (2020), which found local pockets of air pollution within the HMA that are not detected by existing sparse‐density air quality networks. Our finding that the co‐occurrence of such features with localized maxima of urban heat and social vulnerability highlights the necessity for advanced data fusion methodologies that incorporate satellite data, numerical modeling, and low‐cost sensor inputs for urban heat and gaseous pollution mapping (Pochwała et al., 2020; Sulzer et al., 2022).

We find synergies arising from population growth, heightened social vulnerability, UHI expansion, and localized PM2.5 increases within the UPI exacerbate environmental inequalities in the HMA, with implications amid climate change concerns. A recent study links approximately 700 yearly deaths to nationwide heat stress (Vaidyanathan et al., 2004), and since Anderson and Bell (2011) report ∼2.5% mortality risk rises per 1°F (∼0.55°C) heatwave intensity, heat stress‐induced deaths may increase as heatwaves are projected to occur more frequently (Meehl & Tebaldi, 2004; Perkins et al., 2012).

Our regression analysis indicates that urban growth may amplify heatwave events in the HMA by an additional 0.55 K, elevating risks. Combined exposure to heat stress and particulate pollution can result in a 250% mortality risk increase from individual effects (Analitis et al., 2018; X. Hu et al., 2022; Stafoggia et al., 2023). Our analysis reveals urban growth contributes to localized particulate pollution. Numerous studies highlight higher risks for socially vulnerable populations from heat stress and particulate pollution (Khatana et al., 2022; O’Lenick et al., 2019). Thus, our findings of increased local social vulnerability with urban growth, alongside heightened heat stress and pollution, are concerning, especially with the projected increase in heatwave events, as was observed in the summer of 2023.

The United Nations (2022) projection that 87% of all Americans and 68% of the world population will live in urban areas by 2050 provides the perspective for considering the local and broader implications of our study. Further, the EPA Integrated Climate and Land Use Scenarios (ICLUS) data set suggests that the urban land cover within the HMA will continue expanding during the upcoming decades. Unless steps are taken to mitigate the UHI and UPI effects, their footprints will likely grow also. Therefore, if the observed relationships between urban social vulnerability, heat stress, and air quality continue, environmental inequalities will be exacerbated.

Our finding that environmental inequalities are worsened by the combined effects of urban heat islands (UHI), urban pollution islands (UPI), and social vulnerability is supported by prior studies (Y. Li et al., 2020; Nair et al., 2023; Sabrin et al., 2020; Tang et al., 2017) and is relevant to other U.S. cities and megacities worldwide. We also found that HMA did not experience increase in green spaces that mitigate UHI (Figure S3c in Supporting Information S1), highlighting the need of adopting such UHI mitigation strategies.

Our study demonstrates the utility of data fusion, incorporating satellite observations for high‐resolution mapping of the UHI and UPI to identify environmental inequalities and associated human health risks. This data fusion approach can help in the formulation of policies for mitigation and can easily be applied to other cities experiencing comparable urban growth.

Supporting Information

The supplemental information to this study, with figures and tables denoted with a leading “S” in the manuscript text when referenced, are provided in a separate document.

Conflict of Interest

The authors declare no conflicts of interest relevant to this study.

Supporting information

Supporting Information S1

Acknowledgments

This research was supported by NASA Grant 80MSFC22M0004. A. Blackford, T. Cowan, U. Nair, C. Phillips, A. Kaulfus, and B. Freitag have confirmed no financial, personal, or professional affiliations that could be perceived as having influenced the work reported in this study. No third‐party financial support was received to complete this work, and there are no relevant financial activities or intellectual property rights outside the submitted work. The authors have no additional relationships or activities that could be perceived to have influenced or potentially influenced the content of this manuscript.

Data Availability Statement

All data sets utilized for the analysis presented in this manuscript are publicly available. Software utilized for data analysis are Python, R, ArcGIS Pro, and Google Earth Engine (Gorelick et al., 2017). The land cover data was sourced from the National Land Cover Database (NLCD) for the years 2001 and 2019 (C. G. Homer et al., 2007; C. H. Homer et al., 2012; C. Homer et al., 2015). Land surface temperature, vegetation cover, and aerosol optical depth data were derived from MODIS products, specifically MOD11A2, MOD13Q1, and MCD19A2 (Didan et al., 2015; Lyapustin & Wang, 2018; Wan, 2014). Surface meteorological data was obtained from NOAA’s National Centers for Environmental Information (NCEI), and PM2.5 data was collected from EPA monitoring stations. Social vulnerability was assessed using the CDC’s Social Vulnerability Index (SVI) (Flanagan et al., 2020; Mah et al., 2023). Surface PM2.5 data was further analyzed using the data set provided by Van Donkelaar et al. (2021).
==== Refs
References

Analitis, A. , De’ Donato, F. , Scortichini, M. , Lanki, T. , Basagana, X. , Ballester, F. , et al. (2018). Synergistic effects of ambient temperature and air pollution on health in Europe: Results from the PHASE project. International Journal of Environmental Research and Public Health, 15 (9 ), 1–11. 10.3390/ijerph15091856
Anderson, G. B. , & Bell, M. L. (2011). Heat waves in the United States: Mortality risk during heat waves and effect modification by heat wave characteristics in 43 U.S. communities. Environmental Health Perspectives, 119 (2 ), 210–218. 10.1289/ehp.1002313 21084239
Bullard, R. D. (1983). Solid waste sites and the black Houston community. Sociological Inquiry, 53 (2–3 ), 273–288. 10.1111/j.1475-682X.1983.tb00037.x 11635985
de Souza, A. , dos Santos, C. M. , Ihaddadene, R. , Cavazzana, G. , Abreu, M. C. , de Oliveira‐Júnior, J. F. , et al. (2022). Analysis of extreme monthly and annual air temperatures variability using regression model in Mato Grosso do Sul, Brazil. Modeling Earth Systems and Environment, 8 (1 ), 647–663. 10.1007/s40808-021-01096-6
Didan, K. , Munoz, A. B. , Solano, R. , & Huete, A. (2015). MODIS vegetation index user’s guide (MOD13 series). (pp. 1–38). Vegetation Index and Phenology Lab, The University of Arizona. 10.5067/MODIS/MYD13Q1.006
Emerson, M. O. , Bratter, J. , Howell, J. , & Jeanty, W. (2012). Houston region grows more racially/ethnically diverse, with small declines in segregation: A joint report analyzing census data from 1990, 2000, and 2010, 21.
Environmental Protection Agency . (2019). Our nation’s air: A trends report. Retrieved from https://gispub.epa.gov/air/trendsreport/2019/#home
Flanagan, B. E. , Gregory, E. W. , Hallisey, E. J. , Heitgerd, J. L. , & Lewis, B. (2020). A social vulnerability index for disaster management. Journal of Homeland Security and Emergency Management, 8 (1 ). 10.2202/1547-7355.1792
Fry, R. , & Taylor, P. (2012). The rise of residential segregation by income (Vol. 202 , p. 26). Pew Research Center. Retrieved from http://www.pewsocialtrends.org/2012/08/01/the‐rise‐of‐residential‐segregation‐by‐income/
Gorelick, N. , Hancher, M. , Dixon, M. , Ilyushchenko, S. , Thau, D. , & Moore, R. (2017). Google Earth engine: Planetary‐scale geospatial analysis for everyone. Remote Sensing of Environment, 202 , 18–27. 10.1016/j.rse.2017.06.031
Homer, C. , Dewitz, J. , Yang, L. , Jin, S. , Danielson, P. , Xian, G. , et al. (2015). Completion of the 2011 national land cover database for the conterminous United States—Representing a decade of land cover change information. Photogrammetric Engineering & Remote Sensing, 81 (5 ), 345–354. 10.14358/PERS.81.5.345
Homer, C. G. , Dewitz, J. , Fry, J. , Coan, M. , Hossain, N. , Larson, C. , et al. (2007). Completion of the 2001 national land cover database for the conterminous United States. Photogrammetric Engineering & Remote Sensing, 73 (4 ), 337–341.
Homer, C. H. , Fry, J. A. , & Barnes, C. A. (2012). The national land cover database. US Geological Survey Fact Sheet, 3020 (4 ), 1–4.
Hu, L. , & Brunsell, N. A. (2013). The impact of temporal aggregation of land surface temperature data for surface urban heat island (SUHI) monitoring. Remote Sensing of Environment, 134 , 162–174. 10.1016/j.rse.2013.02.022
Hu, X. , Han, W. , Wang, Y. , Aunan, K. , Pan, X. , Huang, J. , & Li, G. (2022). Does air pollution modify temperature‐related mortality? A systematic review and meta‐analysis. Environmental Research, 210 , 112898. 10.1016/j.envres.2022.112898 35181304
Kendall, M. G. (1948). Rank correlation methods.
Khatana, S. A. M. , Werner, R. M. , & Groeneveld, P. W. (2022). Association of extreme heat with all‐cause mortality in the contiguous US, 2008‐2017. JAMA Network Open, 5 (5 ), E2212957. 10.1001/jamanetworkopen.2022.12957 35587347
Lane, H. M. , Morello‐Frosch, R. , Marshall, J. D. , & Apte, J. S. (2022). Historical redlining is associated with present‐day air pollution disparities in U.S. Cities. Environmental Science and Technology Letters, 9 (4 ), 345–350. 10.1021/acs.estlett.1c01012 35434171
Li, H. , Meier, F. , Lee, X. , Chakraborty, T. , Liu, J. , Schaap, M. , & Sodoudi, S. (2018). Interaction between urban heat island and urban pollution island during summer in Berlin. Science of the Total Environment, 636 , 818–828. 10.1016/j.scitotenv.2018.04.254 29727848
Li, Y. , Sun, Y. , Li, J. , & Gao, C. (2020). Socioeconomic drivers of urban heat island effect: Empirical evidence from major Chinese cities. Sustainable Cities and Society, 63 , 102425. 10.1016/j.scs.2020.102425
Lim, C. H. , Ryu, J. , Choi, Y. , Jeon, S. W. , & Lee, W. K. (2020). Understanding global PM2.5 concentrations and their drivers in recent decades (1998–2016). Environment International, 144 , 106011. 10.1016/j.envint.2020.106011 32795749
Logan, T. M. , Anderson, M. J. , Williams, T. G. , & Conrow, L. (2021). Measuring inequalities in urban systems: An approach for evaluating the distribution of amenities and burdens. Computers, Environment and Urban Systems, 86 , 101590. 10.1016/j.compenvurbsys.2020.101590
Lyapustin, A. , & Wang, Y. (2018). MCD19A2 MODIS/Terra+Aqua land aerosol optical depth daily L2G global 1km SIN grid V006 [Dataset]. NASA EOSDIS Land Processes Distributed Active Archive Center. 10.5067/MODIS/MCD19A2.006
Mah, J. C. , Penwarden, J. L. , Pott, H. , Theou, O. , & Andrew, M. K. (2023). Social vulnerability indices: A scoping review. BMC Public Health, 23 (1 ), 1–11. 10.1186/s12889-023-16097-6 36624437
Mann, H. B. (1945). Nonparametric tests against trend. Econometrica: Journal of the Econometric Society, 13 (3 ), 245. 10.2307/1907187
Meehl, G. A. , & Tebaldi, C. (2004). More intense, more frequent, and longer lasting heat waves in the 21st century. Science, 305 (5686 ), 994–997. 10.1126/science.1098704 15310900
Miller, D. J. , Actkinson, B. , Padilla, L. , Griffin, R. J. , Moore, K. , Lewis, P. G. T. , et al. (2020). Characterizing elevated urban air pollutant spatial patterns with mobile monitoring in Houston, Texas. Environmental Science and Technology, 54 (4 ), 2133–2142. 10.1021/acs.est.9b05523 31995368
Moser, A. , Uhl, E. , Rötzer, T. , Biber, P. , Dahlhausen, J. , Lefer, B. , & Pretzsch, H. (2017). Effects of climate and the urban heat island effect on urban tree growth in Houston. Open Journal of Forestry, 07 (04 ), 428–445. 10.4236/ojf.2017.74026
Nair, U. S. , Phillips, C. , Dasgupta, S. , & Mitra, C. (2023). Impact of growth of a medium‐sized Indian coastal city on urban climate: A case study using data fusion and analytics. Urban Climate, 49 , 101525. 10.1016/j.uclim.2023.101525
Nicholson, J. (2017). Segregated city: The geography of economic segregation in America’s metros. Journal of Real Estate Literature, 25 (2 ), 454.
Nowak, D. J. , & Walton, J. T. (2005). Projected urban growth (2000–2050) and its estimated impact on the US forest resource. Journal of Forestry, 103 (8 ), 383–389. 10.1093/jof/103.8.383
O’Dell, K. , Hornbrook, R. S. , Permar, W. , Levin, E. J. T. , Garofalo, L. A. , Apel, E. C. , et al. (2020). Hazardous air pollutants in fresh and aged western US wildfire smoke and implications for long‐term exposure. Environmental Science and Technology, 54 (19 ), 11838–11847. 10.1021/acs.est.0c04497 32857515
O’Lenick, C. R. , Wilhelmi, O. V. , Michael, R. , Hayden, M. H. , Baniassadi, A. , Wiedinmyer, C. , et al. (2019). Urban heat and air pollution: A framework for integrating population vulnerability and indoor exposure in health risk analyses. Science of the Total Environment, 660 , 715–723. 10.1016/j.scitotenv.2019.01.002 30743957
Perkins, S. E. , Alexander, L. V. , & Nairn, J. R. (2012). Increasing frequency, intensity and duration of observed global heatwaves and warm spells. Geophysical Research Letters, 39 (20 ), 1–5. 10.1029/2012GL053361
Pochwała, S. , Gardecki, A. , Lewandowski, P. , Somogyi, V. , & Anweiler, S. (2020). Developing of low‐cost air pollution sensor—Measurements with the unmanned aerial vehicles in Poland. Sensors, 20 (12 ), 1–17. 10.3390/s20123582
Sabrin, S. , Karimi, M. , & Nazari, R. (2020). Developing vulnerability index to quantify urban heat islands effects coupled with air pollution: A case study of Camden, NJ. ISPRS International Journal of Geo‐Information, 9 (6 ), 1–14. 10.3390/ijgi9060349
Sen, P. K. (1968). Estimates of the regression coefficient based on Kendall’s Tau. Journal of the American Statistical Association, 63 (324 ), 1379–1389. 10.1080/01621459.1968.10480934
Sexton, K. , & Linder, S. H. (2015). Houston’s novel strategy to control hazardous air pollutants: A case study in policy innovation and political stalemate. Environmental Health Insights, 9 (s1 ), 1–12. 10.4137/EHI.S15670
Sexton, K. , Linder, S. H. , Marko, D. , Bethel, H. , & Lupo, P. J. (2007). Comparative assessment of air pollution‐related health risks in Houston. Environmental Health Perspectives, 115 (10 ), 1388–1393. 10.1289/ehp.10043 17938725
Southerland, V. A. , Brauer, M. , Mohegh, A. , Hammer, M. S. , van Donkelaar, A. , Martin, R. V. , et al. (2022). Global urban temporal trends in fine particulate matter (PM2·5) and attributable health burdens: Estimates from global datasets. The Lancet Planetary Health, 6 (2 ), e139–e146. 10.1016/S2542-5196(21)00350-8 34998505
Stafoggia, M. , Michelozzi, P. , Schneider, A. , Armstrong, B. , Scortichini, M. , Rai, M. , et al. (2023). Joint effect of heat and air pollution on mortality in 620 cities of 36 countries. Environment International, 181 , 108258. 10.1016/j.envint.2023.108258 37837748
Streutker, D. R. (2002). A remote sensing study of the urban heat island of Houston, Texas. International Journal of Remote Sensing, 23 (13 ), 2595–2608. 10.12677/gst.2017.52004
Streutker, D. R. (2003). Satellite‐measured growth of the urban heat island of Houston, Texas. Remote Sensing of Environment, 85 (3 ), 282–289. 10.1016/S0034-4257(03)00007-5
Sulzer, M. , Christen, A. , & Matzarakis, A. (2022). A low‐cost sensor network for real‐time thermal stress monitoring and communication in occupational contexts. Sensors, 22 (5 ), 1828. 10.3390/s22051828 35270974
Tang, J. , Di, L. , Xiao, J. , Lu, D. , & Zhou, Y. (2017). Impacts of land use and socioeconomic patterns on urban heat Island. International Journal of Remote Sensing, 38 (11 ), 3445–3465. 10.1080/01431161.2017.1295485
Theeuwes, N. E. , Barlow, J. F. , Teuling, A. J. , Grimmond, C. S. B. , & Kotthaus, S. (2019). Persistent cloud cover over mega‐cities linked to surface heat release. Npj Climate and Atmospheric Science, 2 (1 ), 15. 10.1038/s41612-019-0072-x
Thiel, H. (1950). A rank‐invariant method of linear and polynomial regression analysis. Indagationes Mathematicae.
Ulpiani, G. (2021). On the linkage between urban heat island and urban pollution island: Three‐decade literature review towards a conceptual framework. Science of the Total Environment, 751 , 141727. 10.1016/j.scitotenv.2020.141727 32890803
United Nations . (2022). 2022 revision of world population prospects. Population Division of the Department of Economic and Social Affairs of the United Nations Secretariat. Retrieved from https://population.un.org/wpp
U.S. Census Bureau . (2022). Quick facts Houston City, Texas. Retrieved from http://www.census.gov/quickfacts/table/PST045215/4835000
Vaidyanathan, A. , Malilay, J. , Schramm, P. , & Saha, S. (2004). MMWR—Heat‐related deaths—United States, 2004–2018. Morbidity and Mortality Weekly Report, 69 (24 ), 729–734. Retrieved from https://www.cdc.gov/nchs/data/dvs/2e_volume1_2017.pdf
Van Donkelaar, A. , Hammer, M. S. , Bindle, L. , Brauer, M. , Brook, J. R. , Garay, M. J. , et al. (2021). Monthly global estimates of fine particulate matter and their uncertainty. Environmental Science and Technology, 55 (22 ), 15287–15300. 10.1021/acs.est.1c05309 34724610
Vo, T. T. , Hu, L. , Xue, L. , Li, Q. , & Chen, S. (2023). Urban effects on local cloud patterns. Proceedings of the National Academy of Sciences of the United States of America, 120 (21 ), 1–11. 10.1073/pnas.2216765120
Wan, Z. (2014). New refinements and validation of the collection‐6 MODIS land‐surface temperature/emissivity product. Remote Sensing of Environment, 140 , 36–45. 10.1016/j.rse.2013.08.027
Waren, W. (2013). Houston’s emerging exposure between African Americans and whites: Evidence of spatial assimilation or place stratification? Population, Space and Place, 19 (5 ), 633–643. 10.1002/psp.1753
Yuvaraj, R. M. (2020). Extents of predictors for land surface temperature using multiple regression model. The Scientific World Journal, 2020 , 1–10. 10.1155/2020/3958589
