Abstract
BACKGROUND: Anopheles mosquitoes pose notable public health concerns as competent vectors of malaria and other diseases. Although malaria is no longer endemic in the United States, recent locally acquired cases in states including Maryland highlight the need to better understand Anopheles dynamics in the region. This study aimed to identify geographic hotspots of Anopheles presence in Maryland and evaluate how land cover, microclimatic conditions, and socioeconomic vulnerability shape their spatial and temporal distribution. METHODS: Monthly Anopheles occurrence data (1999-2024) from Global Biodiversity Information Facility (GBIF) were aggregated at county and Census Block Group (CBG) scales. Counties were ranked by mean annual presence to identify hotspots. Associations with land cover, microclimatic, and socioeconomic conditions were assessed using Spearman's rank correlation (ρ; P < 0.05). At the CBG scale, significantly correlated variables (|ρ| ≥ 0.25) were used to fit an Extreme Gradient Boosting model to quantify the relative importance of environmental and socioeconomic predictors. Spatial dependence was addressed through blocked cross-validation, and model interpretability was evaluated with SHapley Additive exPlanations (SHAP) values. RESULTS: Prince George's and Anne Arundel Counties emerged as primary hotspots of presence with highest observer effort during the analysis period. Seasonal analysis revealed an annual cycle, with peak presence from May to September, coinciding with warmer conditions favorable to vector proliferation. SHAP analysis at the CBG-scale identified habitat availability as the most influential predictor (33.1% of total model impact for low impervious surface percentage) with woody wetland emerging as the most preferred habitat; followed by humid conditions (24.6%), and low elevation (18.2%). Notably, cooler and more humid microclimates within the warm season provide optimal habitat, reflecting fine-scale environmental controls on Anopheles distribution. Therefore, CBG-level analysis within Prince George's County revealed a negative correlation between Area Deprivation Index and Anopheles presence (ρ = -0.35), indicating fine-scale ecological drivers-such as woody wetland habitat, low impervious surface, and humid cooler microclimates-more prevalent in affluent suburban residential neighborhoods. CONCLUSIONS: This study demonstrates that fine-scale habitat characteristics and warm-season microclimates structure Anopheles mosquito presence in Maryland. These insights support more spatially targeted vector control and improved public health surveillance strategies.