Pests and disease have become an increasingly common issue as globalized trade brings non-native species into unfamiliar systems. Emerald ash borer (Agrilus planipennis), is an Asiatic species of boring beetle currently devastating the native population of ash (Fraxinus) trees in the northern forests of the United States, with 85 million trees having already succumbed across much of the Midwest. We have developed a reaction-diffusion partial differential equation model to predict the spread of emerald ash borer over a heterogeneous 2-D landscape, with the initial ash tree distribution given by data from the Forest Inventory and Analysis. As expected, the model predictions show that emerald ash borer consumes ash which causes the local ash population to decline, while emerald ash borer spreads outward to other areas. Once the local ash population begins to decline emerald ash borer also declines due to the loss of available habitat. Our model's strength lies with its focus on the county scale and its linkage between emerald ash borer population growth and ash density. This enables one to make accurate predictions regarding emerald ash borer spread which allows one to consider various methods of control as well as to accurately study the economic effects of emerald ash borer spread.