Anomalous diffusion exists widely in polymer and biological systems. Pulsed-field gradient (PFG) anomalous diffusion is complicated, especially in the anisotropic case where limited research has beenreported. An general PFG signal attenuation expression, including the finite gradient pulse (FGPW) effect for free general anisotropic fractional diffusion { 0 < α, β ≤ 2 } based on the fractional derivative, has not been obtained, where α and β are time and space derivative orders. It is essential to derive a general PFG signal attenuation expression including the FGPW effect for PFG anisotropic anomalous diffusion research.In this paper, two recently developed modified-Bloch equations, the fractal differential modified-Bloch equation and the fractional integral modified-Bloch equation, were extended to obtain general PFG signal attenuation expressions for anisotropic anomalous diffusion. Various cases of PFG anisotropic anomalous diffusion were investigated, including coupled and uncoupled anisotropic anomalous diffusion. The continuous-time random walk (CTRW) simulation was also carried out to support the theoretical results.There is good agreement between the theory and the CTRW simulation. The obtained signal attenuation expressions and the three-dimensional fractional modified-Bloch equations are important for analyzing PFG anisotropic anomalous diffusion in NMR and MRI., the anomalous diffusion reduces to normal diffusion. Unlike normal diffusion, anomalous diffusion has a non-Gaussian propagator, and its mean β-th power of displacement is not linearly proportional to its diffusion time. Anomalous diffusion could be modeled by time-space fractional diffusion equations based on fractal derivative (see Appendix A) model and fractional derivative (see Appendix B) model [1,2,11,12,13,14,15,16,17,18,19]. Due to the non-Gaussian characteristic features of anomalous diffusion, it is difficult to investigate pulsed field gradient (PFG) anomalous diffusion. The PFG technique has been a powerful tool in studying normal diffusion [20,21,22,23,24], but many PFG normal diffusion theories may not be directly applicable to the investigation of anomalous diffusion without modification.There have been many efforts devoted to studying PFG anomalous diffusion [18,19,25,26,27,28,29,30,31,32,33,34,35,36], but most of them are related to isotropic anomalous diffusion. The theoretical research of PFG anisotropic anomalous diffusion is very limited [37,38,39]. G. Lin, General PFG signal attenuation expressions for anisotropic anomalous diffusion by……, 2017 2 (a) 90 G RF 180 δ PGSE (b) 90 G RF δ PGSTE 90 90 Fig. 1 Two typical PFG pulses: (a) PGSE pulse sequences, (b) PGSTE pulse sequence. The gradientpulse width is , and the diffusion delay is Δ. In PGSTE pulse sequence, the magnetization is transferred to the Z direction by the second radio frequency (RF) pulse to eliminate the T2 relaxation effect, which has the advantage over the PGSE pulse sequence.Anisotropic diffusion behaviors could exist in many materials [40,41,42,43,...