Three growing metros, Houston, Goiânia and Kelowna, and the climate they are growing into
The same tools in three countries: how much each metro grew, what the growth took, and what the climate scenarios indicate for heat, drought, fire and heavy rain. Read in scenarios, not in forecasts: the signal is broad, and it is the place that decides what is exposed.
Loading the numbers…
—inhabitants in the Houston metro area, 1 July 2025 measured
—inhabitants in the Goiânia metropolitan region, IBGE 2025 estimate measured
—inhabitants in the Kelowna metropolitan area, 1 July 2025 measured
—times Goiânia’s built-up area in 30 years, 25 km circle (1990–2020) measured
measured by an official table, a station or a satellite.
estimated a reading that depends on a choice (a grid, a model, a threshold).
missing a part that did not respond in this run.
1
Three metros, three countries, the same tools
Houston, in the United States, Goiânia, in Brazil, and Kelowna, in Canada, are metros that grow fast in each country’s official numbers, and each already has a climate story told by UmaTerra. In 2024–25, Houston was the metropolitan area that gained the most people in the United States. This story reads the three with the same satellite data and with UmaTerra’s tools: how much each grew, what the growth took, and what the climate scenarios indicate for heat, drought, fire and heavy rain. It reads in scenarios, not in forecasts.
2
What the official tables say
Each country defines its metropolitan area differently: in the United States a county-based statistical area; in Canada a census metropolitan area, built from municipalities linked by commuting; in Brazil a metropolitan region set by law, whose list of municipalities IBGE publishes. The three series below come from each country’s official tables and cover different periods; they are not the same ruler.
3
What the satellite sees
To the left of the bar, the 50 km area around the centre seen by Landsat 5 in the mid-1980s; to the right, by Sentinel-2 in 2025. The lines show each metro’s legal window and the fixed 10, 25 and 50 km circles, the story’s common ruler. Pick a city in the selector and drag the bar.
4
How much the built-up area grew
The Global Human Settlement Layer (GHSL) measures built-up surface in 100 m cells, from 1975 to 2030. This story uses only the observed epochs, up to 2020: the 2025 and 2030 ones are projections and are left out. The chart puts the three cities on the same base, 1990 equal to 100, to compare pace; the table gives the hectares.
5
Up or outwards?
Dividing population by built-up surface gives the density of the urban footprint. To have a continuous series we use GHSL’s modelled population; it is an estimate, and the table sets it against the official numbers where the years coincide. Density that falls while the city grows points to outward expansion; density that holds, to filling in.
6
Which ring took the growth
The map shows, by colour, the epoch in which each 100 m cell became built-up, up to 2020. The chart adds the new built-up surface between 1990 and 2020 in four rings of distance from the centre. It is the most direct reading of where the city grew: near the centre, in the middle, or on the far edge.
7
Land use, in each country, from 1985 to 2025
Each country has its own land-use map: MapBiomas in Brazil, the Annual NLCD in the United States and the Annual Crop Inventory in Canada. They are not the same product and do not add up. The table shows, in a 10 km circle, land use at the start and at the end of the series in two places of each city: an urban one and a rural one. In Canada the inventory is agricultural, starts in 2011 and changed its wild-land classes between versions; so only city, water and forest are read with confidence.
8
The canopy that went, and why we cannot say why
Hansen counts all tree-cover loss since 2001, from felling, fire, storm or pests. So it is not a deforestation map: next to the growing city it also records Kelowna’s 2023 fire and the felling of plantations. The map shows loss by period; the chart, hectares per year in the 25 km circle.
9
The scenario is broad; the place decides what is exposed
Climate models see in cells of about 25 km: a city’s whole circle gets the same signal, and in Kelowna the two places fall in the same cell. What changes from one place to the other is what is on the ground. The table crosses, at each place, the land use with what the models indicate for the cell and shows how much of the area is sensitive to each hazard. It is a way of saying that these broad numbers need local information to be read, and that none of them is a forecast for your street.
10
The property and the parcel, up close
Camposcope analyses a rural property from Brazil’s register (the CAR, self-declared) or a parcel from British Columbia’s register. Here it reads four examples: in Goiânia, one property on the urban edge and one rural; in Kelowna, one urban parcel and one rural, both owned by a local government, so that no one’s home is exposed. The United States is not in Camposcope yet: coverage is being expanded, and Houston is left out of this page.
11
Hot days since 1981, in what was measured
The ERA5-Land reanalysis gives the daily maximum temperature since 1950 in cells of about 9 km. For each year we count the days when the maximum passed 30 °C, which serves as a common ruler for the three cities, and those that passed the threshold each country uses. This is what has already happened, before any scenario.
12
Three countries, three thresholds
UmaTerra’s climate tool uses, for a hot day, 30 °C in Canada, 32.2 °C (90 °F) in the United States and 35 °C in Brazil: in Canada 35 °C almost never occurs, and the threshold has to mean something locally. The count of days by each threshold measures different things and cannot be compared between cities. Nor can it be compared between tools: the scenarios use a bias-corrected 25 km grid, the reanalysis a 9 km one, and the number of hot days from one does not hold in the other. What is compared is the change in degrees, and that is what the scenarios page uses.
13
The core against the ring, by satellite
MODIS measures land-surface temperature since 2000, in 1 km pixels, day and night. We compare the 5 km core around the centre with a 20–30 km ring, in each city’s warm season, without water pixels. It is a measured surface difference, not an air one, and the ring may have different terrain and vegetation: it does not attribute the heat to the city’s growth.
14
Exposure, not prediction
Nine CMIP6 climate models, downscaled by NASA, compare 2041–2070 with a 1985–2014 baseline under two emission scenarios: the intermediate (SSP2-4.5) and the high (SSP5-8.5). A hazard counts as present when the model median passes a threshold and at least two-thirds agree on its sign; as strong, at a higher threshold with at least 80%. The result describes what the models indicate for a cell, never a given year’s weather nor anyone’s street.
15
How much warmer, and how many more days
The change in mean temperature and in the hottest day of the year, in degrees, is the measure comparable across the three cities. The table also gives each country’s hot days, which are read only within their own city. The global frame is the IPCC’s: 2.0 °C above 1850–1900 under the intermediate scenario and 2.4 °C under the high one, in 2041–2060.
16
Less rain, more evaporation, more dry days
Drought is read through three numbers: dry days per year, annual rainfall, and the water balance, which is rain minus potential evaporation. A warmer year with the same rain ends up drier; so the balance can worsen even where rain does not change. The hazard appears as present or strong by the rule on the previous page.
17
Fire seen by satellite and fire weather in the scenarios
The map shows the years each 500 m pixel burned between 2001 and 2025, from MODIS. The chart adds what burned per year in the 25 km circle. For the future, the Canadian Fire Weather Index (FWI) counts days of high and extreme danger in the scenarios. Fire weather is a condition, not a fire, and heat and drought drive it in many regions.
18
The water that has stayed and what the models indicate
The map shows the water the satellite saw in at least a quarter of images since 1984 and, in Houston, the radar reading of Harvey’s 2017 flood. For the future, the models give the change in the heaviest rainy day and in heavy-rain days; a warmer atmosphere holds more water vapour, and extreme rains tend to intensify. Houston’s and Kelowna’s events have their own stories, linked at the end.
19
Which place today each city’s 2041–2070 climate resembles
To picture a future climate, the tool looks for the cell whose present climate is closest to what the models indicate for the city: mean temperature, annual rainfall and dry days. The reading compares three annual means, not extremes, seasons or soil; it says which place the climate resembles, not that the land will look like that place.
20
A strong El Niño, and what past seasons show
NOAA’s Climate Prediction Center tracks an El Niño under way, strong by the relative index (RONI). For each city, the table groups the October–March seasons since 1981 by ENSO phase and shows the mean rain and how many were below normal. There are few strong events, and a single one can carry the mean. The models also do not agree on how ENSO itself will change in a warmer climate.
21
Heat, drought, fire and heavy rain, in both scenarios
The matrix summarises, for each city’s cell, the level of each hazard (0 no signal, 1 present, 2 strong) in the two scenarios, with how many of the nine models agree on the sign of the change. It shows where the signal is robust and where it is weak; it does not say what will happen in a given year or where inside the cell.
22
Drawing places prone to each hazard
Inside each city’s 25 km circle, three strata: lowland (height above drainage up to 3 m), repeated fire (burned in at least two years between 2001 and 2025) and heat (the hottest tenth of the daytime surface temperature). From each stratum a draw with a fixed seed picks four points, and Naturametrics’ buffer service reads a 1 km circle around each: land cover, built-up area, temperature, years burned and height above drainage. An empty stratum stays empty.
23
Three cities, the same questions, different answers
The final table gathers what the story measured, city by city, without ranking them: growth, the shape of growth, the heat that has already come, the intermediate scenario’s signal and the heat island. The differences come from the place, from the definition of a metro and from the data; the story does not claim any city is better or worse prepared.
umaterra.com.br · Houston · Goiânia · Kelowna
Loading…measured
Loading the numbers…
What this story does not say
Not a forecast. The scenarios describe what nine models indicate for a 25 km cell in 2041–2070, not a given year’s weather nor what will happen in a neighbourhood. A city’s two places are told apart only by what is on the ground, and in Kelowna they fall in the same cell.
The three metros are not defined the same way. A county-based statistical area, a commuting-based census region and a legal metropolitan region have different sizes and logics; that is why the story’s common ruler is the fixed 10, 25 and 50 km circles.
GHSL built-up and population are modelled. Only the epochs up to 2020 are used; earlier ones come from lower-resolution imagery and built-up is a surface per cell, not a map of buildings. Population is an estimate and differs from the official count.
Hot-day thresholds change by country. The day count is by each country’s threshold and cannot be compared between cities nor between the scenarios’ grid and the reanalysis; the change in degrees is what is compared.
Surface temperature is not air temperature. MODIS measures the ground and the canopy top in 1 km pixels; the 20–30 km ring has its own terrain and vegetation (in Kelowna, hills and forest). The difference does not attribute the heat to the city’s growth.
Canopy loss is not deforestation, and each country has its own land-use map. Hansen counts fire, felling, storms and pests; MapBiomas, the NLCD and AAFC’s inventory are different products, and the last is agricultural and changed classes between versions: in Canada only city, water and forest are read with confidence.
El Niño has few strong events in the sample. Five strong El Niños since 1981: a single one can carry the mean. The models do not agree on how ENSO will change.
Camposcope and the draws are examples, not samples. They are four properties, two per city; the CAR is self-declared and Camposcope does not yet cover the United States. The draws have a documented seed and design, but four points per stratum do not represent the city; in Houston the up-to-3 m lowland covers most of the circle and does not discriminate, and the repeated-fire stratum is empty.
What is not here. People’s vulnerability, infrastructure, adaptation and costs: the story measures growth and climate exposure, and leaves the question of who is most exposed, and why, to local work.
Written with Anthropic’s Claude Sonnet 5.5, from the tables and papers cited and the numbers Earth Engine and the UmaTerra tools computed: no fact or number comes from the model. Versão em português.
Notes and references
US Census Bureau. Metropolitan and Micropolitan Statistical Areas Population Totals: 2020–2025 (Vintage 2025, file cbsa-est2025-alldata.csv), Houston-Pasadena-The Woodlands, TX (CBSA 26420). Boundaries: cartographic boundary file 2023. census.gov.
Statistics Canada. Table 17-10-0148-01, Population estimates, July 1, by census metropolitan area and census agglomeration, 2021 boundaries; Kelowna CMA (59915), with the nine municipalities and reserves of the statistical classification (SGC 2016) and the 2021 census polygons. statcan.gc.ca.
IBGE. 2010 and 2022 Demographic Censuses (SIDRA, tables 202 and 4714), 2024 and 2025 population estimates (table 6579) and Composition of the Metropolitan Regions 2025 (21 municipalities of the Goiânia RM); 2022 municipal mesh. The 2021 estimate, on a pre-census base, is not used. sidra.ibge.gov.br.
Pesaresi, M., & Politis, P. (2023). GHS-BUILT-S R2023A: GHS built-up surface grid, derived from Sentinel2 composite and Landsat, multitemporal (1975–2030). European Commission, JRC. doi:10.2905/9F06F36F-4B11-47EC-ABB0-4F8B7B1D72EA. Read in Earth Engine (JRC/GHSL/P2023A/GHS_BUILT_S), 100 m cells; epochs 1975–2020 observed, 2025 and 2030 projected and not used. The map counts a cell as built with ≥ 500 m² of built-up surface.
Schiavina, M., Freire, S., & MacManus, K. (2023). GHS-POP R2023A: GHS population grid multitemporal (1975–2030). European Commission, JRC. doi:10.2905/2FF68A52-5B5B-4A22-8F40-C41DA8332CFE. Population modelled from census data disaggregated by built-up area: an estimate, not a count.
Imagery: Landsat 5 TM Collection 2 Level 2 (USGS/NASA), median of 1984–1986 (Houston, Goiânia) and of June to September 1984–1991 (Kelowna, where snow scenes fill the rest of the year); Sentinel-2 harmonised SR (Copernicus/ESA), median of May to September 2025, clouds and shadows masked. Read in Earth Engine.
Hansen, M. C. et al. (2013). High-resolution global maps of 21st-century forest cover change. Science 342, 850–853. doi:10.1126/science.1244693. Version v1.13 (loss 2001–2025), read in Earth Engine; tree-cover loss from any cause.
Land use as muda-clima reads it: MapBiomas Collection 11 (Brazil, 1985–2025), USGS Annual NLCD (conterminous United States, 1985–2025) and AAFC’s Annual Crop Inventory (Canada, 2009–2025, national since 2011; an agricultural inventory whose wild-land classes changed between versions). muda-clima’s exposure groups.
Thrasher, B. et al. (2022). NASA Global Daily Downscaled Projections, CMIP6. Scientific Data 9, 262. doi:10.1038/s41597-022-01393-4; Eyring, V. et al. (2016). Overview of the Coupled Model Intercomparison Project Phase 6 (CMIP6) experimental design and organization. Geoscientific Model Development 9, 1937–1958. doi:10.5194/gmd-9-1937-2016. Nine models (EC-Earth3, MPI-ESM1-2-HR, CMCC-ESM2, MRI-ESM2-0, GFDL-ESM4, MIROC6, NorESM2-MM, IPSL-CM6A-LR, INM-CM5-0), 0.25° cells and the NEX-GDDP-FWI fire-weather index; muda-clima’s precomputed grid.
IPCC (2021). Summary for Policymakers. In: Climate Change 2021: The Physical Science Basis, Table SPM.1 (global warming in 2041–2060 above 1850–1900: 2.0 °C under SSP2-4.5 and 2.4 °C under SSP5-8.5, best estimates). Cambridge University Press. doi:10.1017/9781009157896.001.
Muñoz-Sabater, J. et al. (2021). ERA5-Land: a state-of-the-art global reanalysis dataset for land applications. Earth System Science Data 13, 4349–4383. doi:10.5194/essd-13-4349-2021. Daily maximum temperature (ECMWF/ERA5_LAND/DAILY_AGGR), ~9 km cell under each city’s centre, 1981–2025.
Wan, Z. (2014). New refinements and validation of the collection-6 MODIS land-surface temperature/emissivity product. Remote Sensing of Environment 140, 36–45. doi:10.1016/j.rse.2013.08.027. MOD11A2 v061, 1 km, warm-season means (Houston and Kelowna: June–August; Goiânia: August–October), water pixels (occurrence ≥ 50%) removed; Peng, S. et al. (2012). Surface urban heat island across 419 global big cities. Environmental Science & Technology 46, 696–703. doi:10.1021/es2030438; Oke, T. R. (1982). The energetic basis of the urban heat island. Quarterly Journal of the Royal Meteorological Society 108, 1–24. doi:10.1002/qj.49710845502.
Giglio, L. et al. (2018). The Collection 6 MODIS burned area mapping algorithm and product. Remote Sensing of Environment 217, 72–85. doi:10.1016/j.rse.2018.08.005. MCD64A1 v061, 500 m; a pixel counts as burned in a year if it burned in any month.
Pekel, J.-F. et al. (2016). High-resolution mapping of global surface water and its long-term changes. Nature 540, 418–422. doi:10.1038/nature20584. Water occurrence (JRC/GSW1_4). Yamazaki, D. et al. (2019). MERIT Hydro: a high-resolution global hydrography map based on latest topography datasets. Water Resources Research 55, 5053–5073. doi:10.1029/2019WR024873. Height above the nearest drainage (HAND).
Cai, W. et al. (2014). Increasing frequency of extreme El Niño events due to greenhouse warming. Nature Climate Change 4, 111–116. doi:10.1038/nclimate2100; Cai, W. et al. (2021). Changing El Niño–Southern Oscillation in a warming climate. Nature Reviews Earth & Environment 2, 628–644. doi:10.1038/s43017-021-00199-z.
NOAA Climate Prediction Center. ENSO Diagnostic Discussion and RONI table, read by muda-clima on 2 Nov 2026 (state: El Niño, RONI of +1.69 in Jul–Sep 2026). TerraClimate (rain, ~4 km) by ENSO phase, October to March, 1981–2025. cpc.ncep.noaa.gov.
Fischer, E. M., & Knutti, R. (2015). Anthropogenic contribution to global occurrence of heavy-precipitation and high-temperature extremes. Nature Climate Change 5, 560–564. doi:10.1038/nclimate2617; Jolly, W. M. et al. (2015). Climate-induced variations in global wildfire danger from 1979 to 2013. Nature Communications 6, 7537. doi:10.1038/ncomms8537.
van Oldenborgh, G. J. et al. (2017). Attribution of extreme rainfall from Hurricane Harvey, August 2017. Environmental Research Letters 12, 124009. doi:10.1088/1748-9326/aa9ef2. The radar reading of the flood (Sentinel-1, 30 Aug 2017) is the one in the Houston Harvey story.
Zanaga, D. et al. (2022). ESA WorldCover 10 m 2021 v200. doi:10.5281/zenodo.7254221. Global land cover, the same in the three countries, for the drawn plots.
The UmaTerra tools used here: muda-clima (scenarios, land use, climate analogue and ENSO), Camposcope (the CAR in Brazil, ParcelMap BC in Canada) and Naturametrics’ buffer service; there is no public description of the tools beyond the apps’ own pages. muda-clima’s hot-day thresholds: 30 °C Canada, 32.2 °C United States, 35 °C Brazil.
Definitions of a metropolitan area: United States, core-based statistical area (OMB; county-based); Canada, census metropolitan area (Statistics Canada; urban core plus municipalities with strong commuting links); Brazil, a metropolitan region set by state law and listed by IBGE. The three are not equivalent.