We investigate the effect of an anisotropic substrate on the turbulent dynamics of a simulated two dimensional active nematic. This is introduced as an anisotropic friction and an effective anisotropic viscosity, with the orientation of the anisotropy being defined by the substrate. In this system we observe the emergence of global nematic order of topological defects that is controlled by the degree of anisotropy in the viscosity and the magnitude of the active stress. No global defect alignment is seen in passive liquid crystals with anisotropic viscosity or friction confirming that ordering is driven by the active stress. We then closely examine the active flow generated by a single defect to show that the kinetic energy of the flow is orientation dependent, resulting in a torque on the defect to align them with the anisotropy in the substrate.Active nematic liquid crystals are fluids consisting of self (or mutually) propelling rod shaped particles resulting in an anisotropic fluid with broken rotational symmetry that drives itself at the microscopic scale [1,2]. This interesting combination of broken rotational symmetry and out of equilibrium, active behaviour has lead to an explosion of interest from both experimental and theoretical physics [1][2][3][4][5][6]. There have been many successful experiments reproducing active nematics often utilising biological components, including microtubule kinesin suspensions [1, 2, 7], acto-myosin gels [8,9] and elongated cells [10,[12][13][14], but also from inert components such as vibrated monolayers of granular rods [15].These systems have been shown to display a rich phenomenology depending on many factors such as the degree to which the system is driven [4], the confining geometry [16], the density [13] and the boundary conditions [17]. By varying these factors it is possible to observe diverse spatiotemporal patterns including vortices [14,16], oscillating textures [8,18] and travelling bands [8,16]. When the driving force is sufficiently high, active nematics can spontaneously nucleate many topological defects, generating flows and interacting chaotically in a regime referred to as low reynolds number, active turbulence [1,4,10,11]. These defects have been shown to exert an elastic torque on each other [19] and experiments have indicated that long range nematic order of defects is possible in a state of active turbulence [20] though this hasn't been reproduced theoretically. It has been shown that the position and orientation of these defects can be influenced by the substrate on which the active nematic is placed. By changing the geometry of the substrate it is possible to reorient defects and sort them by charge [21], by changing the topology of the substrate it is possible to control the total number of defects and their trajectories [18]. A defect ordered active nematic has been recreated by placing a 2D active nematic on top of a passive liquid crystal that can be controlled by an external magnetic field [22]. When the passive liquid crystal layer is ordered ...