1The dependency of wave velocity in reaction-diffusion (RD) systems on the local front curvature determines not only the stability of wave propagation, but also the fundamental properties of other spatial configurations such as vortices. This Letter gives the first derivation of a covariant eikonal-curvature relation applicable to general RD systems with spatially varying anisotropic diffusion properties, such as cardiac tissue. The theoretical prediction that waves which seem planar can nevertheless possess a nonvanishing geometrical curvature induced by local anisotropy is confirmed by numerical simulations, which reveal deviations up to 20% from the nominal plane wave speed. [5]. The waves of excitation in cardiac tissue that initiate cardiac contraction are also being intensively studied, as abnormal wave propagation in the heart may lead to the formation of vortices and results in dangerous cardiac arrhythmias [6,7]. The key parameter that characterizes a propagating wave is its velocity; it turns out that the velocity of wave propagation in RD systems is affected by several factors, the curvature of the wave front being one of the most important. The dependency of the wave velocity on front curvature, i.e., the eikonal-curvature relation, was found in many systems [4,[8][9][10], and it was also shown to be essential to the stability of the wave front and properties of vortices in the particular RD system. The simplest eikonal-curvature relation is given by the following linear relationship [8,9]:with c 0 the speed of a traveling plane wave in the given medium and a medium-dependent constant close to the scalar diffusivity of the isotropic medium [8,11]. Although stability analysis of Eq. (1) ensures stable propagating waves for positive , stability in critical media where lies close to zero is not guaranteed [12]. Remarkably, even when > 0, the linearized equation (1) Fortunately, the mathematical treatment of local anisotropy is greatly simplified when the curved-space formalism introduced in [15] is adopted. This formalism will not only allow us to generalize Eq. (1) to homogeneous excitable media of an arbitrary anisotropy type, but it also simplifies calculations such that the quadratic curvature corrections to the wave speed can be calculated for the first time. Our mathematical derivation provides an alternative and a generalization to Keener's seminal proof for the isotropic case, which relies on a boundary layer approximation and is thus restricted to the limit of steep wave fronts [9]. Here, Kuramoto's approach [8] is elaborated in a Riemannian context and pursued to higher order in curvature.To start the derivation, the set of coupled reactiondiffusion equations (RDE) is written in terms of a column matrix u of state variables u ðmÞ , with the summation convention for repeated spatial indices, @ t uðr; tÞ ¼ @ i ðD ij ðrÞ@ j Puðr; tÞÞ þ Fðuðr; tÞÞ:Local anisotropy of the medium is embodied by the spatially varying diffusion tensor D ij ðrÞ. The constant projection matrix P selects only those sta...