Let A ∈ Sym(n × n) be an elliptic 2-tensor. Consider the anisotropic fractional Schrödinger operator L s A +q, where L s A := (−∇ · (A(x)∇)) s , s ∈ (0, 1) and q ∈ L ∞ . We are concerned with the simultaneous recovery of q and possibly embedded soft or hard obstacles inside q by the exterior Dirichlet-to-Neumann (DtN) map outside a bounded domain Ω associated with L s A + q. It is shown that a single measurement can uniquely determine the embedded obstacle, independent of the surrounding potential q. If multiple measurements are allowed, then the surrounding potential q can also be uniquely recovered. These are surprising findings since in the local case, namely s = 1, both the obstacle recovery by a single measurement and the simultaneous recovery of the surrounding potential by multiple measurements are longstanding problems and still remain open in the literature. Our argument for the nonlocal inverse problem is mainly based on the strong uniqueness property and Runge approximation property for anisotropic fractional Schrödinger operators.