Nucleon form factors at large momentum transfer are important for understanding the transition from nonperturbative to perturbative QCD and have been the focus of experiment and phenomenology. We calculate proton and neutron electromagnetic form factors GE,M (Q2) from first principles using nonperturbative methods of lattice QCD. We have accumulated large Monte Carlo statistics to study form factors up to momentum transfer Q2 ≲ 8 GeV2 with a range of lattice spacings as well as quark masses that approach the physical point. In this paper, results of initial analyses are presented and compared to experiment, and potential sources of systematic uncertainty are discussed.