Planetary-scale answers, unlocked.
A Hands-On Guide for Working with Large-Scale Spatial Data. Learn more.
Authors
Data center site selection is the process of finding land for a facility of tens to hundreds of megawatts. On suitable land, the facility can be powered, connected, cooled, and permitted, and kept clear of natural hazards. Analysts screen candidate land against grid capacity, transmission, fiber, water, hazards, zoning, and latency to users. They then rank the land that passes on cost and time to build. This article tests one claim: power limits where a data center can go, and the rest of the screen finds suitable land near power.
Most of those criteria are geographic. The screen combines map layers in a geographic information system (GIS) and ranks land with multi-criteria decision analysis (MCDA), a method that turns several criteria into one score.
Hyperscalers (the largest cloud operators), colocation operators (who rent out data center space), and their advisors run these screens for new campuses. Utilities run them in reverse to find where they can serve a large new load. Investors use them to check portfolios for hazard exposure.
The stakes come from demand. Masanet and colleagues estimated that global data centers used about 205 TWh (terawatt hours) in 2018, around 1% of world electricity. The estimate appears in Recalibrating global data center energy-use estimates (Science, 2020). The Lawrence Berkeley National Laboratory 2024 United States Data Center Energy Usage Report by Shehabi and colleagues puts U.S. data center use at 176 TWh in 2023. That was about 4.4% of U.S. electricity. The report projects 325 to 580 TWh by 2028. The International Energy Agency's Energy and AI report projects global use rising from about 415 TWh in 2024 to about 945 TWh in 2030.
Each criterion maps to data layers and to either a hard filter (a constraint that removes land) or a score (a factor that ranks it):
Power capacity and interconnection queues. A hyperscale campus needs firm capacity, power the utility commits to deliver at all hours. The grid connection usually sets the schedule. Generators also wait in interconnection queues, the lists of projects waiting to connect to the grid. Lawrence Berkeley National Laboratory's Queued Up: 2024 Edition by Rand and colleagues counted the queues at the end of 2023. They held nearly 2,600 GW (gigawatts; a gigawatt is a million kilowatts) of generation and storage. The typical project built in 2023 took nearly five years from request to commercial operation.
Transmission lines and substations. Distance to a high-voltage line and a substation with spare capacity sets the cost of the connection. New lines take years to build.
Fiber. Data centers need diverse fiber routes to carriers and other data centers. Fiber maps are mostly commercial, so open-data screens use highway and rail corridors as a rough stand-in.
Water and cooling. Evaporative cooling consumes water, and climate sets how many hours a year outside air can cool the servers. In dry regions water becomes a hard constraint.
Land and zoning. A campus needs a large, flat, contiguous parcel zoned for the use, or a realistic path to rezoning. Slope comes from a digital elevation model, wetlands and forest from land cover, and the legal unit is the parcel.
Natural hazards. The Federal Emergency Management Agency (FEMA) National Flood Hazard Layer maps the 1% and 0.2% annual chance flood zones that most screens exclude (see flood risk data). Seismic hazard comes from the U.S. Geological Survey (USGS) 2023 National Seismic Hazard Model described by Petersen and colleagues (Earthquake Spectra, 2024). Wildfire and wind layers complete the hazard screen.
Latency to population. Latency is the delay for data to travel between users and servers. A common rough rule divides the straight-line distance times 2.1 by the speed of light in fiber. Bozkurt and colleagues cite it in Dissecting latency in the Internet's fiber infrastructure (2018). Workloads that need fast responses stay near users and internet exchanges. AI training can sit farther away.
Tax incentives and land cost. State sales and use tax exemptions on equipment, local property tax rates, and land prices enter as scores.
The data center site selection process turns a region into a ranked shortlist:
Each layer enters through a spatial join, which matches records by location, the core operation of geospatial analysis for siting. Distance to transmission is a distance join. Flood exposure is a point-in-polygon test. Slope comes from zonal statistics, a summary of elevation values inside each cell or parcel.
Weights from the Analytic Hierarchy Process. Thomas Saaty introduced the Analytic Hierarchy Process (AHP) in A scaling method for priorities in hierarchical structures (Journal of Mathematical Psychology, 1977). Decision makers compare criteria two at a time on a 1 to 9 scale. A matrix calculation, the principal eigenvector, turns the comparisons into weights. Saaty's rule of thumb is to revisit any matrix whose consistency ratio, a check for contradictory judgments, is above 0.10.
WherobotsDB runs a first-pass screen as one SQL query over open layers in the Havasu catalog. The worked example is an illustrative screen of eastern Loudoun County, Virginia. The study box runs from 77.75°W to 77.30°W and 38.85°N to 39.15°N. It covers Leesburg, Ashburn, Sterling, Dulles Airport, and Chantilly. The unit of analysis is the H3 resolution 8 hexagon, under a square kilometer and close to the size of a campus parcel. ST_H3CellIDs returns the 1,808 cells whose centers fall in the box.
Latency is left out on purpose: every cell sits inside the Washington metro area, so it barely varies between cells. Population from Census block groups and American Community Survey (ACS) estimates belongs in a regional screen that compares metros.
Overture power, building, and road layers feed the four factors, and the Overture water layer feeds the constraint. Each Overture table is stored in Havasu, Wherobots' spatial table format built on Apache Iceberg. Row counts cover each table's search box in the query:
ST_H3CellIDs
base_infrastructure
buildings_building
transportation_segment
base_water
The water layer stands in for FEMA flood zones, which a production screen loads from the National Flood Hazard Layer (NFHL) through the ArcGIS Feature Service reader.
The weights are illustrative. A real study sets them with stakeholders, for example through AHP, and tests the ranking's sensitivity to them.
The inventory query counts the power features that intersect the box, grouped by each feature's highest tagged voltage. It ran in 91 seconds on a Small SQL session:
-- Power lines and substations in the study box, by highest tagged voltage SELECT 'power' AS layer, class, CASE WHEN kv >= 500000 THEN '500 kV' WHEN kv >= 230000 THEN '230 kV' WHEN kv IS NULL THEN 'untagged' ELSE 'under 230 kV' END AS band, COUNT(*) AS n FROM (SELECT class, array_max(transform(split(source_tags['voltage'], ';'), v -> CAST(v AS INT))) AS kv FROM wherobots_open_data.overture_maps_foundation.base_infrastructure WHERE bbox.xmax > -77.75 AND bbox.xmin < -77.30 AND bbox.ymax > 38.85 AND bbox.ymin < 39.15 AND subtype = 'power' AND class IN ('power_line', 'substation')) GROUP BY 1, 2, 3 ORDER BY 2, 3;
Each criterion becomes a column per cell, and the score is the weighted sum:
-- Illustrative multi-criteria screen for data center sites on H3 resolution 8 cells WITH cells AS ( -- 1,808 cells over eastern Loudoun County and its neighbors, centers projected to UTM zone 18N meters SELECT /*+ REPARTITION(8) */ h3, ST_Transform(ST_Centroid(ST_H3ToGeom(array(h3))[0]), 'EPSG:4326', '+proj=utm +zone=18 +datum=WGS84 +units=m +no_defs') AS cm FROM (SELECT explode(ST_H3CellIDs(ST_PolygonFromEnvelope(-77.75, 38.85, -77.30, 39.15), 8, false)) AS h3) ), hv AS ( -- Power lines and substations tagged 230 kV or more (OpenStreetMap voltage tags, in volts) SELECT class, ST_Transform(geometry, 'EPSG:4326', '+proj=utm +zone=18 +datum=WGS84 +units=m +no_defs') AS gm FROM wherobots_open_data.overture_maps_foundation.base_infrastructure WHERE bbox.xmax > -77.82 AND bbox.xmin < -77.23 AND bbox.ymax > 38.80 AND bbox.ymin < 39.20 AND subtype = 'power' AND class IN ('power_line', 'substation') AND array_max(transform(split(source_tags['voltage'], ';'), v -> CAST(v AS INT))) >= 230000 ), roads AS ( SELECT ST_Transform(geometry, 'EPSG:4326', '+proj=utm +zone=18 +datum=WGS84 +units=m +no_defs') AS gm FROM wherobots_open_data.overture_maps_foundation.transportation_segment WHERE bbox.xmax > -77.82 AND bbox.xmin < -77.23 AND bbox.ymax > 38.80 AND bbox.ymin < 39.20 AND subtype = 'road' AND class IN ('motorway', 'trunk') ), wet AS ( -- Constraint: every cell a mapped river or lake touches (a flood-plain proxy that stands in for FEMA flood zones) SELECT DISTINCT explode(ST_H3CellIDs(geometry, 8, true)) AS h3 FROM wherobots_open_data.overture_maps_foundation.base_water WHERE bbox.xmax > -77.75 AND bbox.xmin < -77.30 AND bbox.ymax > 38.85 AND bbox.ymin < 39.15 AND subtype IN ('river', 'lake') ), built AS ( -- Building footprint area per cell, each building assigned to the cell holding its centroid SELECT ST_H3CellIDs(ST_Centroid(geometry), 8, false)[0] AS h3, SUM(ST_AreaSpheroid(geometry)) AS m2 FROM wherobots_open_data.overture_maps_foundation.buildings_building WHERE bbox.xmax > -77.76 AND bbox.xmin < -77.29 AND bbox.ymax > 38.84 AND bbox.ymin < 39.16 GROUP BY 1 ), -- Distance from each cell center to the nearest feature, with the small side broadcast d_sub AS (SELECT /*+ BROADCAST(h) */ c.h3, MIN(ST_Distance(c.cm, h.gm)) AS m FROM cells c JOIN hv h ON h.class = 'substation' AND ST_DWithin(c.cm, h.gm, 5000) GROUP BY c.h3), d_line AS (SELECT /*+ BROADCAST(h) */ c.h3, MIN(ST_Distance(c.cm, h.gm)) AS m FROM cells c JOIN hv h ON h.class = 'power_line' AND ST_DWithin(c.cm, h.gm, 3000) GROUP BY c.h3), d_road AS (SELECT /*+ BROADCAST(r) */ c.h3, MIN(ST_Distance(c.cm, r.gm)) AS m FROM cells c JOIN roads r ON ST_DWithin(c.cm, r.gm, 5000) GROUP BY c.h3), f AS ( SELECT c.h3, ds.m AS sub_m, dl.m AS line_m, dr.m AS road_m, COALESCE(bu.m2, 0) / ST_AreaSpheroid(ST_H3ToGeom(array(c.h3))[0]) AS built_share, w.h3 IS NOT NULL AS wet FROM cells c LEFT JOIN d_sub ds ON ds.h3 = c.h3 LEFT JOIN d_line dl ON dl.h3 = c.h3 LEFT JOIN d_road dr ON dr.h3 = c.h3 LEFT JOIN built bu ON bu.h3 = c.h3 LEFT JOIN wet w ON w.h3 = c.h3 ), s AS ( -- Factor scores: 1, 0.5, or 0 SELECT *, CASE WHEN sub_m <= 2000 THEN 1 WHEN sub_m <= 5000 THEN 0.5 ELSE 0 END AS s_sub, CASE WHEN line_m <= 1000 THEN 1 WHEN line_m <= 3000 THEN 0.5 ELSE 0 END AS s_line, CASE WHEN built_share < 0.02 THEN 1 WHEN built_share < 0.05 THEN 0.5 ELSE 0 END AS s_land, CASE WHEN road_m <= 2000 THEN 1 WHEN road_m <= 5000 THEN 0.5 ELSE 0 END AS s_road FROM f ) -- Illustrative weights: substation 0.30, 230 kV+ line 0.25, open land 0.30, highway 0.15 SELECT lower(hex(h3)) AS h3, ROUND(sub_m) AS sub_m, ROUND(line_m) AS line_m, ROUND(road_m) AS road_m, ROUND(built_share, 4) AS built_share, wet, s_sub, s_line, s_land, s_road, CASE WHEN wet THEN NULL ELSE ROUND(0.30 * s_sub + 0.25 * s_line + 0.30 * s_land + 0.15 * s_road, 3) END AS score FROM s ORDER BY score DESC NULLS LAST, sub_m;
The distance joins copy the small layers (substations, lines, highways) to every worker. EXPLAIN FORMATTED shows three such broadcast index joins and no Cartesian product, the slow all-pairs comparison. Buildings and water join on the H3 cell ID, so no polygon is intersected with another. The one query built the grid, read four tables, and ran eight joins: three distance joins and five joins on the H3 cell ID. It scored all 1,808 cells in 13.9 seconds on a Small SQL session.
EXPLAIN FORMATTED
To count cells per factor score, and to test the weights, the summary query reuses the main query. It keeps every named subquery, or common table expression (CTE), through s and adds a comma after it. It then replaces the final SELECT with the scored CTE and SELECT below. It ran in 12.5 seconds:
s
SELECT
scored
scored AS ( SELECT *, CASE WHEN wet THEN NULL ELSE 0.30 * s_sub + 0.25 * s_line + 0.30 * s_land + 0.15 * s_road END AS score, CASE WHEN wet THEN NULL ELSE 0.25 * (s_sub + s_line + s_land + s_road) END AS score_equal FROM s ) SELECT crit, lvl, COUNT(*) AS cells FROM (SELECT stack(8, 'substation', CAST(s_sub AS DOUBLE), 'line_230kv', CAST(s_line AS DOUBLE), 'open_land', CAST(s_land AS DOUBLE), 'highway', CAST(s_road AS DOUBLE), 'near_water_excluded', CAST(IF(wet, 1, 0) AS DOUBLE), 'score', ROUND(CAST(score AS DOUBLE), 3), 'score_equal', ROUND(CAST(score_equal AS DOUBLE), 3), 'both_085', CAST(IF(ROUND(score, 3) >= 0.85 AND ROUND(score_equal, 3) >= 0.85, 1, 0) AS DOUBLE)) AS (crit, lvl) FROM scored) GROUP BY crit, lvl ORDER BY crit, lvl;
Each criterion removes a different part of the map:
The river and lake constraint gives 395 cells a null score, which leaves 1,413 scored. 41 cells score 1.0, full marks on every factor. The top cell, 882aaabb19fffff near Arcola, sits 61 m from a 230 kV substation, 117 m from a 230 kV line, and 253 m from a highway. Buildings cover 0.57% of it.
The open-land factor penalizes the existing cluster. Cells in the substation ring around Ashburn and Sterling score mid-range because buildings already cover them. The Dulles runways count as open land. The shortlist goes to parcel, zoning, and ownership checks next.
The ranking holds up when the weights change. With equal weights of 0.25 on every factor, 112 cells score 0.85 or more under both weightings, and the same 41 cells keep full marks. The step scores of 1, 0.5, and 0 cause ties, which the fuzzy scoring of Jiang and Eastman (cited below) reduces.
A production screen replaces each stand-in with the authoritative layer. Parcels come from a File Geodatabase, and slope from GLO-30, the Copernicus global surface model with 30 m cells through RS_ZonalStats. The Site Selection use case shows the pattern, and an AI agent can run it through the Wherobots Model Context Protocol (MCP) Server. These joins add facts about physical places to geospatial AI models trained on text, documents, databases, and the internet.
The production version covered 52 counties in the Far West and West weather zones of the Electric Reliability Council of Texas (ERCOT), which runs most of the Texas grid. It ran as one Wherobots job run in 648 seconds on a Medium runtime. The unit is the parcel: 907,845 Regrid records reduce to 871,593 unique polygons. Setbacks remove land near water, pipelines, airports, schools, hospitals, and settled areas. A slope screen on GLO-30 removes land steeper than 5%. That leaves 153,910 parcels with at least 4 ha, and 30,690,000 buildable acres.
The job read the parcels, the Overture tables, and the GLO-30 raster in WherobotsDB: raster and vector in one engine; most stacks pick one.
Distance to power uses an estimated line route, the straight-line distance times 1.32, the routing allowance in an independent balanced model built on Wherobots for Australia. The transmission cost estimation guide from the Midcontinent Independent System Operator (MISO) estimates line length the same way, adding 30% to the straight-line distance between substations. Within a 5 km estimated route, 5,570,879 buildable acres lie near substations of 138 kV or more. 1,192,666 acres, 3.9% of the buildable total, lie near 345 kV.
The balanced model puts 0.40 on distance to 138 kV or more and 0.25 on distance from sensitive receptors (schools, hospitals, and homes). It puts 0.20 on distance to a wastewater feature, and 0.15 on assemblage size (adjoining parcels with one owner). A power-first weighting puts 0.45 on 345 kV and 0.15 on 138 kV, and the other three factors share 0.40 in the same ratio. The two top 50 lists share 9 parcels. The balanced list spans 8 counties and the power-first list 21. Part of that shift is built in, since a heavier 345 kV weight lifts parcels near 345 kV. The Loudoun top set held under equal weights; the West Texas top 50 did not.
The wastewater layer needs a source check. The study box holds 160 Overture water features classed as wastewater, all unnamed OpenStreetMap polygons. In Reeves and Winkler counties they may be oilfield ponds. All 50 balanced top parcels sit within 10 km of one, and 32 sit in those two counties.
Planned power changes the ranking too. A timed score treats the end substations of planned lines and generators with signed agreements as power sources from their planned dates. At the end of 2027, it changes 9 of the power-first top 50. The dates are plans.
Data center siting borrows from fifty years of land suitability research.
Ian McHarg's Design with Nature (1969) made the stacked map overlay the planner's tool. Collins, Steiner, and Rushman trace its path to GIS in Land-use suitability analysis in the United States (Environmental Management, 2001).
Steve Carver joined multi-criteria evaluation to GIS for siting in Integrating multi-criteria evaluation with geographical information systems (International Journal of Geographical Information Systems, 1991). Eastman, Jin, Kyem, and Toledano added constraint maps, factor maps, and weighted linear combination to raster GIS in Raster procedures for multi-criteria/multi-objective decisions (Photogrammetric Engineering and Remote Sensing, 1995). The same paper allocated land among competing objectives.
Jacek Malczewski wrote On the use of weighted linear combination method in GIS (Transactions in GIS, 2000). The paper found that analysts often applied the weighted overlay without a full understanding of its assumptions. Two of them shape the Loudoun screen. Every factor sits on the same 0 to 1 scale. Each weight sets how much of one factor can make up for a shortfall in another. His survey of GIS-based multicriteria decision analysis (International Journal of Geographical Information Science, 2006) classified more than 300 refereed articles published from 1990 to 2004.
Jiang and Eastman replaced hard thresholds with fuzzy membership functions in Application of fuzzy measures in multi-criteria evaluation in GIS (2000). A fuzzy membership function scores land on a smooth curve from 0 to 1, which removes the sudden score drop at each threshold. Ligmann-Zielinska and Jankowski mapped where suitability depends on the weights in Spatially-explicit integrated uncertainty and sensitivity analysis of criteria weights in multicriteria land suitability evaluation (Environmental Modelling and Software, 2014).
Covas, Silva, and Dias applied MCDA to sustainable data center location in Multicriteria decision analysis for sustainable data centers location (International Transactions in Operational Research, 2013).
No public file gives the number site selection needs most: spare capacity at each substation by year. Open data covers lines, substations, and planned in-service dates. Spare capacity comes from the utility.
Siddik, Shehabi, and Marston mapped the footprint in The environmental footprint of data centers in the United States (Environmental Research Letters, 2021). One-fifth of the direct water footprint of U.S. data center servers comes from moderately to highly water-stressed watersheds. Nearly half of the servers draw power fully or partly from plants in water-stressed regions. Evaporative cooling uses less electricity but consumes water, so the cheapest site to power can be the costliest to cool.
A county screen fits on a laptop. The West Texas screen joined 871,593 parcels to substations, receptors, and an elevation model. Jia Yu, Jinxuan Wu, and Mohamed Sarwat introduced spatial joins across many machines on Apache Spark with GeoSpark (ACM SIGSPATIAL 2015). Yu, Zhang, and Sarwat described its spatial partitioning and indexing in Spatial data management in Apache Spark: the GeoSpark perspective and beyond (GeoInformatica, 2019). GeoSpark became Apache Sedona, the engine under WherobotsDB. It splits those joins across the machines of a cluster.
Both screens show that data center site selection starts with power. Constraints, weights, and source layers shape the shortlist.
Screen land against power, building, road, and hazard layers in WherobotsDB with a Wherobots free trial at cloud.wherobots.com.
Data center sites are chosen in two passes over candidate land. A first screen removes land that fails a hard constraint: inside a flood zone, on steep ground, in a protected area, or too far from high-voltage transmission. A second pass scores the remaining cells or parcels on power, fiber, water, latency, incentives, and cost. A weighted sum of those scores ranks the land, a method called GIS multi-criteria decision analysis. Due diligence and utility discussions then narrow the ranked shortlist to a site.
The main data center site selection criteria are power capacity and the grid interconnection timeline, distance to transmission lines and substations, and fiber connectivity and route diversity. The other main criteria are water for cooling, natural hazard exposure (flood, seismic, wildfire, severe wind), land size, zoning and slope, latency to users, tax incentives, and land cost. Power sets the schedule in the examples this article cites.
A large campus needs firm capacity, power the utility commits to deliver at all hours. The nearest substation may not have spare capacity, and new transmission takes years to permit and build. In 2022 Dominion Energy told data center customers that power for some new facilities in eastern Loudoun County, Virginia, would be delayed for years because transmission had not kept pace.
Data center site selection uses parcel boundaries and zoning, transmission lines and substations with voltage, interconnection queue records, and fiber route maps. It also uses FEMA flood zones and other hazard layers, elevation models for slope, water availability and watershed stress, building footprints, census population, and incentive zone boundaries. Open sources such as Overture Maps and the U.S. Census cover power lines, substations, roads, buildings, and population.
A first-pass geospatial screen takes days once the data is assembled, and a single SQL query can score every cell in a county. Due diligence on the shortlist, including utility discussions, environmental review, and land negotiation, takes longer. Grid connection often takes longest. The typical U.S. generation project built in 2023 took nearly five years from interconnection request to commercial operation, according to Lawrence Berkeley National Laboratory.
The distance is set locally. County zoning ordinances set setbacks (the minimum distance from a building to a property line or a neighbor), noise limits, and screening for data center uses. Backup generators fall under state air permits. Site selection screens often treat nearby homes as a constraint, because noise and community opposition can stop a rezoning.
The five key components of a data center are power (utility supply, backup generators, and uninterruptible power), cooling, network connectivity, physical space for racks and equipment, and security. Site selection finds land where the first four can be supplied at the scale and date a project needs.
Multi-criteria decision analysis (MCDA) combines several criteria into one ranking. In GIS-based MCDA, hard constraints remove land first. Each remaining factor is converted to a common 0 to 1 scale and multiplied by a weight. The sum is a suitability score per cell or parcel, a method called weighted linear combination or weighted overlay. Malczewski’s 2006 survey in the International Journal of Geographical Information Science classified more than 300 articles that use the approach.
The Analytic Hierarchy Process, introduced by Thomas Saaty in 1977 and set out in his 1980 book, derives criteria weights from pairwise comparisons on a 1 to 9 scale. A standard matrix calculation, the principal eigenvector of the comparison matrix, turns the comparisons into weights. A consistency ratio above 0.10 signals judgments that contradict each other and should be revisited. AHP is widely used to weight GIS site selection criteria.
U.S. data centers used about 176 TWh in 2023, about 4.4% of U.S. electricity, according to the Lawrence Berkeley National Laboratory 2024 report. The same report says use could reach 325 to 580 TWh by 2028, or 6.7% to 12% of U.S. electricity. Globally, the International Energy Agency (IEA) projects data center use rising from about 415 TWh in 2024 to about 945 TWh in 2030.
What is a Kelvin Wave? Speed, El Niño and 2026 Data
A Kelvin wave is a long, low ocean wave that carries warm water east along the equator in an El Niño. In the 30 days to 10 November 2026, the sea at San Diego stood a record 41.2 cm above NOAA's tide prediction.
What Is an ETL Pipeline? Steps, Examples and Spatial ETL
An ETL (extract, transform, load) pipeline moves data from sources into storage. Learn the steps with examples. Then see a spatial mistake that passed the usual checks: rounding coordinates emptied 27% of Denver's building outlines.
What is an MCP Server? Spatial Context for AI Agents
An MCP server gives AI agents tools to call. With no tools, Claude Code gave no count of San Francisco's buildings. Via Wherobots MCP, it counted 162,538.
share this article
Awesome that you’d like to share our articles. Where would you like to share it to: