Topological order is now being established as a central criterion for characterizing and classifying ground states of condensed matter systems and complements categorizations based on symmetries. Fractional quantum Hall systems and quantum spin liquids are receiving substantial interest because of their intriguing quantum correlations, their exotic excitations and prospects for protecting stored quantum information against errors. Here we show that the Hamiltonian of the central model of this class of systems, the Toric Code, can be directly implemented as an analog quantum simulator in lattices of superconducting circuits. The four-body interactions, which lie at its heart, are in our concept realized via Superconducting Quantum Interference Devices (SQUIDs) that are driven by a suitably oscillating flux bias. All physical qubits and coupling SQUIDs can be individually controlled with high precision. Topologically ordered states can be prepared via an adiabatic ramp of the stabilizer interactions. Strings of qubit operators, including the stabilizers and correlations along non-contractible loops, can be read out via a capacitive coupling to read-out resonators. Moreover, the available single qubit operations allow to create and propagate elementary excitations of the Toric Code and to verify their fractional statistics. The architecture we propose allows to implement a large variety of many-body interactions and thus provides a versatile analog quantum simulator for topological order and lattice gauge theories.Topological phases of quantum matter [1,2] are of scientific interest for their intriguing quantum correlation properties, exotic excitations with fractional and non-Abelian quantum statistics, and their prospects for physically protecting quantum information against errors [3].The model which takes the central role in the discussion of topological order is the Toric Code [4-6], which is an example of a Z 2 lattice gauge theory. It consists of a two-dimensional spin lattice with quasi local fourbody interactions between the spins and features 4 g degenerate, topologically ordered ground states for a lattice on a surface of genus g. The topological order of these 4 g ground states shows up in their correlation properties. The topologically ordered states are locally indistinguishable (the reduced density matrices of a single spin are identical for all of them) and only show differences for global properties (correlations along non-contractible loops). Moreover local perturbations cannot transform these states into one another and can therefore only excite higher energy states that are separated by a finite energy gap.In equilibrium, the Toric Code ground states are therefore protected against local perturbations that are small compared to the gap, which renders them ideal candidates for storing quantum information [4]. Yet for a selfcorrecting quantum memory, the mobility of excitations needs to be suppressed as well, which so far has only been shown to be achievable in four or more lattice dimensions [7,8]...