29Estimating species divergence with gene flow has been crucial for characterizing the gray 30 zone of speciation, which is the period of time where lineages have diverged but have not yet 31 achieved strict reproductive isolation. However, estimates of divergence times and gene flow 32 often ignores spatial information, for example the formation and shape of hybrid zones. Using 33 population genomic data from the eastern ratsnake complex (Pantherophis obsoletus), we infer 34 phylogeographic groups, gene flow, changes in demography, the timing of divergence, and 35 hybrid zone widths. We examine the spatial context of diversification by linking migration and 36 timing of divergence to the size, shape, and types of hybridization (e.g., F1, backcrosses) in 37 hybrid zones. Rates of migration between lineages are associated with the width and shape of 38 hybrid zones. Timing of divergence is not related to migration rate across species pairs and is 39 therefore a poor proxy for inferring position in the gray zone. Artificial neural network 40 approaches are applied to understand how landscape features and past climate have influenced 41 population genetic structure among these lineages prior to hybridization. The Mississippi River 42 produced the deepest divergence in this complex, whereas Pleistocene climate and elevation 43 secondarily structured lineages. 44 45